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Stimulated  scattering  was  performed  with  infrared  lasers  at  1064  nm  and  with  a  ring  dye  laser  at  532  nm. 
For  the  measurements  at  1064  nm,  a  Nd:YAG  laser  was  used  as  the  pump  laser  and  a  tunable  diode  laser  was  used 
as  the  probe  laser.  With  this  apparatus  stimulated  Brillouin  and  stimulated  Rayleigh  scattering  were  observed  in 
hexane  and  other  liquids.  These  were  the  first  high  resolution  stimulated  scattering  performed  without  a  ring  dye 
laser  and  the  first  stimulated  scattering  gain  spectroscopy  in  the  infrared  to  our  knowledge.  Both  electrostrictive 
and  absorptive  versions  of  Brillouin  and  Rayleigh  scattering  were  obtained.  The  transition  from  electrostrictive  to 
absorptive  Rayleigh  scattering  was  observed  for  the  first  time  to  our  knowledge.  The  influence  of  scattered  pump 
light  was  determined  not  to  be  a  problem  for  our  measurements.  A  supercritical  cell  was  designed  and  constructed 
capable  of  operation  up  to  600  K  and  2000  psi.  A  second  generation  cell  was  built  to  reduce  the  effects  of 
convection  at  the  windows.  The  first  stimulated  scattering  measurements  in  a  supercritical  fluid  were  performed. 
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OBJECTIVES 


The  objectives  of  this  research  are  to  develop  stimulated  scattering  as  a  diagnostic  for 
supercritical  fluids,  and  use  this  technique  to  improve  our  understanding  of  fluids  in  the 
supercritical  state. 

The  study  of  supercritical  fluids  and  flows  requires  new  diagnostic  techniques.  Currently 
available  techniques  such  as  laser-induced  fluorescence  (LIF)  and  coherent  anti-Stokes  Raman 
scattering  (CARS)  are  complicated  by  increased  molecular  interactions,  leading  to  stronger 
quenching,  larger  absorption  and  refractive  index,  and  incomplete  understanding  of  the  influence 
of  local  conditions  on  spectroscopic  parameters  such  as  linewidths,  nonresonant  background 
contributions,  and  quenching  rates.  This  research  shows  that  stimulated  scattering  techniques  hold 
great  promise  for  studying  supercritical  fluids. 


STATUS  OF  EFFORT 

Research  was  performed  on  stimulated  scattering  as  a  diagnostic  for  supercritical  fluids. 
New  diagnostics  are  needed  in  the  supercritical  regime  because  low  pressure  diagnostics  do  not 
work  well.  Using  stimulated  scattering  measurements,  thermal,  compressional,  and  compositional 
properties  of  supercritical  fuels  may  be  determined  in  situ.  These  techniques  should  improve  our 
knowledge  of  fluid  properties  in  the  supercritical  state. 

Stimulated  scattering  was  performed  with  infrared  lasers  at  1064  nm  and  with  a  ring  dye 
laser  at  532  nm.  Because  of  strong  absorption  by  thermally-stressed  jet  fuel  in  the  visible, 
measurements  at  1064  nm  are  preferred.  For  the  measurements  at  1064  nm,  a  Nd:YAG  laser  was 
used  as  the  pump  laser  and  a  tunable  diode  laser  was  used  as  the  probe  laser.  Modifications 
performed  to  allow  sensitive  measurements  in  the  infrared  include  substitution  of  an  InGaAs 
detector  for  a  silicon  detector  and  substitution  of  an  optical  fiber  for  the  pinhole  spatial  filter.  With 
this  apparatus  stimulated  Brillouin  and  stimulated  Rayleigh  scattering  were  observed  in  hexane  and 
other  liquids.  These  were  the  first  high  resolution  stimulated  scattering  performed  without  a  ring 
dye  laser  and  the  first  stimulated  scattering  gain  spectroscopy  in  the  infrared  to  our  knowledge. 
Both  electrostrictive  and  absorptive  versions  of  Brillouin  and  Rayleigh  scattering  were  obtained. 
The  transition  from  electrostrictive  to  absorptive  Rayleigh  scattering  was  observed  for  the  first  time 
to  our  knowledge.  The  influence  of  scattered  pump  light  was  determined  not  to  be  a  problem  for 
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our  measurements.  Scattered  pump  light  had  been  a  significant  concern  because  our  infrared 
measurements  were  performed  with  10  times  less  probe  beam  power  and  over  four  times  smaller 
gain  coefficient  than  previous  measurements  at  532  nm.  A  supercritical  cell  was  designed  and 
constructed  capable  of  operation  up  to  600  K  and  2000  psi.  Titanium  bolts  and  glass  windows 
were  used  to  avoid  leaks  associated  with  differential  expansion  at  high  temperatures.  A  second- 
generation  cell  was  built  to  reduce  the  effects  of  convection  at  the  windows.  The  first  stimulated 
scattering  measurements  in  a  supercritical  fluid  were  performed. 


INTRODUCTION 

Supercritical  fuels  and  mixing  are  becoming  very  important  for  Air  Force  projects  such  as 
advanced  air  breathing  propulsion  systems  and  advanced  launch  boosters.  The  relative  lack  of 
experience  with  and  understanding  of  the  fluid  properties,  fluid  dynamics  and  mixing,  and 
chemical  kinetics  under  such  conditions  requires  the  development  and  application  of  diagnostic 
techniques  that  are  well  suited  to  characterize  the  new  features  and  processes  that  may  occur. 

These  techniques  are  especially  important  because  diagnostic  approaches  for  low  pressures  cannot 
be  easily  adapted  to  supercritical  conditions.  This  research  on  stimulated  scattering  is  very  well 
suited  to  performing  diagnostics  on  supercritical  fluids. 

Supercritical  conditions  are  important  in  two  significant  areas.  For  air  breathing 
propulsion,  increasing  demands  are  being  placed  on  fuel  as  a  coolant  for  the  engine,  environmental 
control  systems,  electronics,  and  at  high  Mach  numbers,  even  the  airframe.  These  cooling  needs 
will  heat  the  fuel  into  the  supercritical  phase.  Relatively  little  is  known  about  the  properties  of  fuels 
in  the  supercritical  regime  or  the  effects  of  injecting  supercritical  fuel  into  subcritical  turbines  or 
ramjets.  Diagnostics  are  required  to  understand  the  state  and  properties  of  the  fuel  as  it  is  heated 
and  to  help  understand  the  effect  of  these  supercritical  fuels  on  engine  performance. 

For  high  pressure  rocket  propulsion,  fuel  or  oxidant  may  transform  from  a  subcritical  to  a 
supercritical  state  during  injection  into  the  combustion  chamber.  For  example,  in  liquid  oxygen- 
hydrogen  rockets,  liquid  oxygen  is  injected  at  supercritical  pressure  and  subcritical  temperature  and 
is  heated  to  supercritical  temperature  in  the  combustion  chamber.  Improved  diagnostics  are 
required  to  study  the  significance  of  this  transcritical  mixing  on  rocket  engine  performance. 

High  temperatures  create  many  unusual  problems  for  fuel  systems.  Thermal  cracking 
changes  fuel  properties  and  can  produce  gums  and  solids  that  freeze  valves  and  foul  fuel  nozzles 
and  heat  exchangers.  In  the  supercritical  regime,  fuels  have  liquid-like  densities  and  gas-like 
diffusivities  and  viscosities.  The  changing  heat  load  during  the  flight  of  advanced  aircraft  will  lead 
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to  widely  varying  fuel  properties  during  a  mission;  fuel  systems  must  perform  well  under  all  these 
conditions.  In  the  region  of  the  critical  point,  destructive  flow  and  pressure  instabilities  can  arise. 
These  various  effects  demand  thorough  investigation  of  proposed  fuels  and  fuel  systems  under 
anticipated  conditions. 

We  performed  research  into  new  diagnostic  techniques  that  can  aid  the  design  and 
diagnostic  evaluation  of  new  supercritical  fuel  systems  for  next  generation  Air  Force  propulsion 
systems.  These  techniques  may  be  used  for  obtaining  a  better  understanding  of  supercritical  fuels 
behavior  and  mixing  and  for  evaluating  fuel  delivery  systems  and  supercritical  fuel  mixing.  This 
information  is  critical  in  ensuring  reliable  performance  of  next  generation  aircraft. 

Optical  Diagnostics  for  Supercritical  Fluids 

Traditional  optical  diagnostic  techniques  can  be  broadly  grouped  into  two  classes: 

(1)  hydrodynamic  diagnostics,  including  schlieren,  interferometry,  and  tracer  particles,  and 

(2)  molecular  diagnostics,  including  laser-induced  fluorescence  (LIF)  and  coherent  anti-Stokes 
Raman  scattering  (CARS).  Techniques  from  these  two  classes  will  continue  to  be  important  in 
characterizing  supercritical  mixing  and  combustion.  Hydrodynamic  diagnostics  will  be  useful  for 
determining  qualitative  behavior,  gross  structure,  and  mapping  of  flow  fields.  Indeed,  this  type  of 
diagnostic  has  been  used  for  several  studies  of  supercritical  mixing.1"4 

Molecular  diagnostics  provide  a  fine-scale  picture  of  the  details  of  molecular  interactions, 
radical  concentrations,  local  temperature,  and  velocities.  Quantitative  interpretation  of  molecular 
diagnostics  remains  a  relative  weakness,  especially  in  the  supercritical  regime.  Extending 
techniques  such  as  LIF  and  CARS  to  high  pressure,  multiphase  environments  is  complicated  by 
increased  molecular  interactions,  leading  to  stronger  quenching,  larger  absorption  and  refractive 
index,  broader  linewidths,  and  incomplete  understanding  of  the  influence  of  local  conditions  on 
spectroscopic  parameters  such  as  linewidths,  nonresonant  background  contributions,  and 
quenching  rates. 

Our  research  involves  a  third  class  of  diagnostics  techniques  that  can  effectively  investigate 
the  unusual  conditions  expected  in  supercritical  fluids.  These  diagnostics  are  based  on  stimulated 
Rayleigh,  Brillouin,  and  Raman  scattering  measurements.  We  call  this  class  collective  diagnostics, 
because  they  involve  collective  material  motions.  These  material  motions  are  in  turn  fundamentally 
linked  to  useful  and  significant  underlying  physical  parameters. 

Near  the  critical  point,  the  gaseous  and  liquid  forms  of  a  fluid  become  indistinguishable. 
Such  a  fluid  has  unusual  properties,  such  as  liquid-like  densities  and  gas-like  diffusivities  and 
viscosities,  with  negligible  surface  tension,  no  heat  of  vaporization,  infinite  isothermal  compress- 
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ibility,  and  enhanced  chemical  reactivity.  The  large  compressibility  and  thermodynamic  equiva¬ 
lence  of  the  liquid  and  gaseous  forms  at  the  same  pressure  mean  that  a  near-critical  fluid  is 
continually  condensing  and  evaporating,  which  really  mean  aggregating  and  dissociating,  because 
there  is  no  actual  separation  of  phases. 

This  brief  description  illustrates  that  near  critical  and  supercritical  fluids  will  experience  a 
range  of  collective  effects  that  may  greatly  influence  the  mechanisms  of  mixing  and  combustion. 
For  example,  the  concept  of  a  droplet  is  no  longer  well  defined,  yet  the  formation  and  breakup  of 
droplets  remain  important  to  understanding  both  how  a  "liquid-like"  fuel  mixes  with  a  "gas-like" 
oxidizer  and  how  chemical  reactions  are  modified.  The  spatial  scale  of  these  collective  effects 
should  be  about  1-1000  pm,  much  larger  than  molecules,  and  difficult  to  visualize  and  quantify 
with  path-integrated  measurements  such  as  schlieren  or  interferometry.  Thus  there  is  a  need  for 
new  diagnostic  techniques,  which  is  filled  by  stimulated  scattering  techniques. 

Stimulated  Scattering 

Rayleigh,  Brillouin,  and  Raman  scattering  occur  commonly  as  spontaneous  scattering. 
These  scattering  processes  arise  from  natural  oscillation  modes  of  materials  and  can  be  used  to 
determine  the  physical  parameters  responsible  for  those  oscillations.  When  these  collective  modes 
are  excited  with  a  powerful  laser,  the  mode  oscillations  can  be  driven  so  hard  that  they  grow 
exponentially.  This  phenomenon  is  called  stimulated  scattering.  The  dominant  advantage  of 
stimulated  scattering  is  that  the  scattered  signal  can  be  made  arbitrarily  large;  otherwise,  these 
processes  produce  extremely  weak  signals.  By  using  a  probe  to  measure  the  induced 
amplification,  we  can  obtain  very  good  quantitative  results.  This  technique  is  distinct  from  the 
stimulated  scattering  that  builds  up  from  noise,  in  which  case  quantification  is  very  difficult. 

The  large  signals  from  stimulated  scattering  are  particularly  helpful  for  investigating 
Rayleigh  and  Brillouin  scattering,  where  the  weak  signals  available  from  spontaneous  scattering 
are  difficult  to  discriminate  from  background  excitation  light.  Other  advantages  of  stimulated 
scattering  include  excellent  temporal  resolution,  and  improved  spectral  resolution  and  signal-to- 
noise  ratio.  Furthermore,  the  use  of  two  laser  beams  allows  spatial  registration  and  point 
measurement  of  local  conditions. 

With  a  single  detection  system,  all  three  processes — Rayleigh,  Brillouin,  and  Raman — can 
be  measured.  These  processes  together  provide  measurements  of  a  wide  range  of  material 
properties.  Rayleigh  scattering  provides  information  on  thermal  properties,  compositional 
fluctuations,  and  density;  Brillouin  scattering  on  compressional  or  elastic  properties  and  density; 
and  Raman  scattering  on  chemical  and  compositional  properties,  density,  and  temperature.  While 
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spontaneous  Brillouin5  and  Raman4-6-7  scattering  have  been  applied  to  supercritical  fluids,  the  use 
of  stimulated  scattering  for  supercritical  fluids  is  new. 

The  basic  arrangement  for  stimulated  scattering  measurements  is  shown  in  Figure  1.  The 
pump  laser  sets  up  an  electric  polarization  oscillating  at  the  characteristic  frequency  of  a  scattering 
mode  of  the  material.8  For  a  strong  input  laser  pulse,  this  polarization  acts  as  a  driving  force, 
leading  to  amplification  of  both  the  material  oscillation  and  a  scattered  optical  wave.  The  optical 
amplification  is  detected  as  a  gain  or  loss  on  the  probe  beam.  Spatial  resolution  is  determined  by 
the  overlap  volume  of  the  pump  and  probe  beams. 


Figure  1 .  Basic  experimental  arrangement  for  stimulated  scattering  measurements. 

When  using  spontaneous  scattering,  it  can  be  quite  difficult  to  discriminate  the  Rayleigh 
and  Brillouin  scattered  signals  against  the  background  excitation  light,  even  under  ideal  detection 
conditions.  However,  the  large  signals  from  stimulated  scattering  processes  result  in  no  problems 
in  discriminating  against  background  light.  This  can  be  seen  in  Figures  4,  5,  9,  and  10  in  the 
Accomplishments  /  New  Findings  section.  Even  though  these  measurements  are  performed  within 
a  few  GHz  (a  fraction  of  1  cm  -1)  of  the  wavelength  of  the  pulsed  excitation  laser,  there  is  no 
contribution  of  scattered  background  light  to  the  measurement. 

Other  advantages  of  stimulated  scattering  over  spontaneous  scattering  include  exceptional 
temporal  resolution,  improved  spectral  resolution  (limited  by  the  laser  linewidths  rather  than  a 
spectrometer  or  interferometer),  and  very  high  signal-to-noise  ratios.  Furthermore,  the  use  of  two 
laser  beams  allows  spatial  registration  and  point  measurement  of  local  conditions. 

Rayleigh  scattering  results  from  refractive  index  variations  due  to  thermal  waves  or 
diffusive  density  fluctuations  and  compositional  fluctuations.  Brillouin  scattering  results  from 
refractive  index  variations  due  to  sound  waves  or  traveling  density  or  pressure  fluctuations. 

In  a  single-component  fluid,  the  positions,  strengths,  line  shapes,  and  linewidths  in  the 
combined  Rayleigh  and  Brillouin  scattering  spectrum  provide  information  about  the  temperature, 
density,  isothermal  compressibility  %T,  the  heat  capacity  ratio  Cp/cv,  the  thermal  diffusivity  DT,  the 
acoustic  attenuation  coefficient  T,  the  adiabatic  sound  velocity  cs,  and  the  longitudinal  kinematic 
viscosity  or  effective  mass  diffusion  coefficient  Dv.9  Values  of  these  parameters  obtained  from 
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stimulated  scattering  in  hexane  are  given  in  Table  1  of  the  Accomplishments  /  New  Findings 
section.  An  additional  oscillatory  mode10-11  provides  information  on  the  exchange  of  energy 
between  the  internal  vibrational  modes  and  translational  modes.  Near  the  critical  point,  several  of 
these  parameters  approach  exceptional  values,  e.g.,  cs  — >  0  and  %T  — >  °°.  Critical  opalescence  is 
an  extreme  case  of  scattering  off  density  fluctuations.  One  of  the  first  practical  uses  of 
spontaneous  Brillouin  scattering  in  gases  was  the  investigation  of  the  speed  of  sound  near  the 
critical  point.5 

In  two-component  fluids,  additional  information  can  be  obtained  from  Rayleigh  and 

Brillouin  scattering  measurements.12  Even  in  a  nominally  well-mixed  sample,  fluctuations  will 

exist  in  relative  concentration  or  mole  fraction.  These  fluctuations  contribute  to  scattering  that 

appears  as  part  of  the  central  Rayleigh  line,  with  a  linewidth  depending  on  the  binary  diffusion 

coefficients  Dy.  In  addition,  the  differences  in  molecular  weight  between  the  two  components 

means  that  there  are  actually  three  speeds  of  sound  in  a  binary  mixture:  the  hydrodynamic  or 

1/2 

adiabatic  speed  of  sound  of  the  mixture,  cs  -  [kTAfn^  +  m2%2)]  >  and  the  speeds  of  sound  of 

the  individual  components,  ~  [kT/mJ  and  ~  [kT/m2]  .  Acoustic  waves  traveling  at  the 

component  speeds  are  rapidly  damped,  but  show  up  as  well-defined  features  at  characteristic 
frequencies  in  the  Brillouin  spectrum  of  the  mixture.  The  linewidths  of  these  features  could  be  used 
to  determine  the  binary  viscosity  coefficients,  T|jj.  In  principle,  such  Brillouin  scattering 
measurements  in  mixtures  could  be  used  to  determine  both  the  local  chemical  composition  and  the 
rate  of  chemical  mixing. 

Raman  scattering  from  molecules  is  a  well  established  technique  in  combustion  diagnostics. 
Spontaneous  Raman  scattering  is  one  of  the  few  optical  diagnostic  techniques  that  has  been  used  to 
date  on  supercritical  fluids.4’6’7’13  Stimulated  Raman  scattering  can  be  performed  with  essentially 
the  same  apparatus  as  the  stimulated  Rayleigh  and  Brillouin  measurements  and  provides  stronger 
signals,  higher  spectral  resolution,  and  better  signal-to-noise  ratios  than  spontaneous  Raman 
scattering.  Information  provided  by  stimulated  Raman  scattering  includes  chemical  composition, 
density,  and  temperature.  Interpretation  of  stimulated  Raman  scattering  measurements  is  simpler 
than  for  CARS,  because  the  signal  is  directly  proportional  to  the  imaginary  part  of  the  third-order 
nonlinear  susceptibility  (or  to  the  Raman  scattering  cross  section)  rather  than  to  the  squared 
magnitude  of  the  sum  of  real  and  imaginary  parts  of  the  third-order  susceptibility,  as  is  the  case  for 
CARS.14 

The  stimulated  scattering  techniques  are  widely  applicable.  We  have  used  stimulated 
Brillouin  and  Rayleigh  scattering  to  perform  measurements  in  glasses  and  crystals8’14’15  without 
the  requirement  of  detailed  spectroscopic  information  on  these  materials.  Our  supercritical 
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experiments  are  aided  by  long  term  research  at  SRI  on  supercritical  fluids,  including  studies  of 
supercritical  water  oxidation16  and  supercritical  fluid  extraction.17 

The  laser  requirements  for  our  stimulated  Rayleigh  and  Brillouin  measurements  are  quite 
modest.  The  pump  beam  energy  is  only  ~  40  |xJ,  and  the  probe  beam  is  provided  by  a  diode  laser. 
Thus,  stimulated  scattering  could  be  used  even  as  an  on  board  monitor  of  fuel  state  and  properties. 

Parameters  currently  of  interest  to  the  Air  Force  for  supercritical  fuels  include 
compressibility,  sonic  velocity,  density  and  phase,  and  chemical  composition.18  All  these 
parameters  can  be  provided  directly  by  stimulated  scattering  measurements.  Thus  stimulated 
scattering  will  play  a  very  important  role  in  understanding  fuel  performance  and  in  diagnosing  fuel 
system  performance,  critical  steps  in  assuring  the  performance  of  next  generation  aircraft.  We  are 
one  of  only  a  handful  of  groups  that  have  performed  stimulated  Rayleigh  and  Brillouin 
experiments,  and  we  are  the  only  group  applying  these  techniques  to  supercritical  fluids. 


ACCOMPLISHMENTS  /  NEW  FINDINGS 

We  applied  stimulated  Rayleigh  and  Brillouin  scattering  to  property  measurements  in  pure 
supercritical  fluids.  The  physical  properties  that  can  be  measured  using  this  technique  are  the 
compressibility  (and  hence  speed  of  sound  and  temperature),  acoustic  damping  times  or  kinematic 
viscosities,  thermal  diffusivity,  and  density.  We  performed  stimulated  scattering  measurements  at 
both  532  nm  and  1064  nm.  The  measurements  at  1064  nm  are  the  first  pump-probe  stimulated 
scattering  measurements  in  the  infrared  to  our  knowledge.  We  made  experimental  modifications  to 
reduce  scattered  light  and  improve  the  signal-to-noise  ratios.  First-  and  second-generation 
supercritical  cells  were  designed,  built,  and  tested.  Measurements  with  high  spectral  resolution 
and  signal-to-noise  ratio  were  obtained  in  supercritical  fluids.  Our  work  complements  well  work 
on  optical  diagnostics  of  supercritical  fuels  being  performed  at  Wright-Patterson  using  spontaneous 
Raman  and  fluorescence  measurements.19 

Stimulated  Scattering  Measurements 

The  basic  experimental  arrangement  for  stimulated  scattering  is  shown  in  Figure  1  above. 
As  a  pump  laser  for  our  stimulated  scattering  measurements,  we  use  an  injection-seeded  Nd.YAG 
laser.  This  laser  was  previously  modified  to  allow  production  of  linewidths  down  to  20  MHz  and 
below.20  This  narrow  linewidth  is  possible  through  generation  of  long  pulse  lengths  produced 
through  lowered  oscillator  gain  and  pulse  lengthening  in  the  amplifier.  Narrow  linewidths  for  the 
pump  laser  are  important  for  resolving  Brillouin  features,  which  can  be  as  narrow  as  13  MHz.8 
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The  optical  properties  of  jet  fuels  in  the  supercritical  state  are  important  for  this  work.  In 
recent  experiments,  fuels  taken  to  supercritical  conditions  are  found  to  be  black  on  returning  to 
ambient  conditions.18  Little  is  known  about  the  optical  properties  of  the  fuels  while  they  are  still  at 
supercritical  conditions.  We  measured  the  transmission  of  thermally  stressed  JP7  fuel  obtained 
from  Tim  Edwards  of  AFRL,  Wright-Patterson.  The  results  are  shown  in  Figure  2.  Although  the 
thermally  stressed  fuel  is  highly  absorbing  in  the  visible,  it  is  transparent  in  the  near  infrared, 
where  we  performed  our  measurements. 

Because  stimulated  scattering  measurements  are  nonresonant,  they  can  be  performed  at  any 
wavelength  for  which  the  medium  is  transparent.  Previous  stimulated  Brillouin  measurements  in 
solids  were  performed  at  532  nm,8>15’21  using  the  second  harmonic  of  an  injection-seeded 
Nd:YAG  laser  as  a  pump  laser.  The  probe  beam  was  produced  with  a  cw  dye  laser  pumped  by  an 
argon  ion  laser.  Recent  advances  in  diode  laser  technology  allow  use  of  a  tunable  diode  laser  as 
the  probe  laser.  In  this  case,  the  fundamental  of  the  Nd:YAG  laser  is  used  as  the  pump  laser  and 
measurements  are  performed  at  1064  nm.  We  chose  to  use  1064  nm  because  of  the  lower 
absorption  in  thermally  stressed  fuel,  the  lower  operating  costs  compared  with  the  argon  laser,  and 
the  reduced  time  required  for  keeping  the  ring  dye  laser  in  good  operation.  Operation  at  532  nm 
has  the  advantage  of  higher  gain  coefficients  and  higher  probe  beam  powers. 


Figure  2.  Transmission  through  1  cm  of  JP7  fuel,  thermally  stressed  to  1200°  F. 
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We  purchased  a  tunable  diode  laser  from  Environmental  Optical  Sensors  Inc.  (EOSI,  since 
acquired  by  Newport  Corporation),  which  had  previously  produced  lasers  at  the  required 
wavelength.  However,  EOSI  had  considerable  difficulty  producing  a  laser  for  us.  The  first  two 
lasers  EOSI  sent  ceased  functioning  prematurely,  and  considerable  delays  were  encountered  before 
receiving  a  third  laser. 

To  maintain  progress  on  our  research  program  while  waiting  for  the  diode  laser,  we 
performed  measurements  at  532  nm,  using  the  apparatus  shown  in  Figure  3.  We  borrowed  an 
argon  ion  laser  and  ring  dye  laser  from  another  researcher  in  our  laboratory.  A  wavemeter  (New 
Focus)  was  used  to  monitor  the  probe  laser  wavelength.  Both  lasers  were  focused  through  the 
sample  using  effective  focal  lengths  of  about  2  m.  To  prevent  saturation  of  the  detector,  the 
chopper  wheel  reduces  the  dye  laser  power  except  during  a  short  period  corresponding  to  the  laser 
pulse.  The  spatial  filter  (lens  and  pinhole  aperture)  reduces  the  amount  of  scattered  pump  laser 
light  incident  on  the  detector. 


Dye 

Laser 


Cell 

Figure  3.  Apparatus  for  stimulated  scattering  at  532  nm. 


With  this  laser  and  the  second  harmonic  of  the  Nd:YAG  laser,  we  obtained  stimulated 
Brillouin  spectra  such  as  that  shown  for  methanol  in  Figure  4.  There  are  two  peaks  in  the 
spectrum.  The  positive  peak  at  the  left  is  a  gain  peak,  which  occurs  when  energy  is  transferred 
from  the  pump  laser  to  the  probe  (dye)  laser.  The  negative  peak  at  the  right  is  a  loss  peak,  which 
occurs  when  energy  is  transferred  from  the  probe  laser  to  the  pump  laser.  Information  that  can  be 
extracted  from  these  peaks  includes  the  overall  density  (from  the  peak  height),  the  acoustic  velocity 
or  compressibility  (from  the  peak  separation),  and  the  acoustic  lifetime  or  kinematic  viscosity  (from 
the  peak  widths). 
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Several  experimental  improvements  were  required  to  obtain  a  stimulated  Rayleigh  scattering 


signal  at  532  nm.  The  dye  laser 
has  large  intensity  fluctuations 
both  near  1  MHz  and  near  80 
MHz.  These  were  reduced  by 
using  low-  and  high-pass  filters 
that  transmitted  most  of  the 
frequency  components  in  the 
pump  pulse.  Amplitude 
fluctuations  associated  with  jitter 
of  the  pump  pulse  were  reduced 
by  triggering  the  boxcar  used  for 
data  acquisition  from  the  laser 


light  rather  than  from  the  Q-switch 


Frequency  (GHz) 


trigger  pulse.  The  data  acquisition  „  . 

Figure  4.  Stimulated  Brillouin  scattering  spectrum  for  methanol 

electronics  were  configured  to 


perform  single-shot  measurements 
to  allow  subtraction  of  the 
scattered  pump  laser  light  on  a 
shot-by-shot  basis.  Following  the 
signal-to-noise  improvements,  we 
were  able  to  observe  stimulated 
Rayleigh  scattering  in  methanol. 
This  peak  is  shown  in  Figure  5 
along  with 

a  fit  to  the  appropriate  theoretical 
lineshape  (the  imaginary  part  of  a 
complex  Lorentzian  function). 

From  our  measurements  at 


Frequency  (GHz) 

Figure  5.  Stimulated  Rayleigh  spectrum  for  methanol. 


532  nm,  we  obtained  enough 

information  to  begin  design  of  a  cell  for  supercritical  experiments.  To  reduce  the  influence  of 
scattered  pump  light,  we  used  a  cell  that  is  long  enough  to  prevent  overlap  of  the  pump  and  probe 
beams  on  the  windows,  and  has  tilted  windows. 
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Supercritical  Cell 

We  designed  and  built  a  cell  for  experiments  at  high  pressure  and  temperature,  up  to  500  K 
and  2000  psi.  The  cell  design  and  high  pressure  system  are  shown  in  Figures  6  and  7.  In 
designing  this  cell,  we  drew  on  information  gained  from  conversations  with  Tim  Edwards  and 
Chris  Bunker  of  AFRL  Wright-Patterson  in  their  work  on  high  pressure/high  temperature  systems. 
The  cell  volume  was  kept  small  to  minimize  the  risk  from  the  failure  of  the  graphite  foil  used  to  seal 
the  windows.  At  supercritical  conditions,  any  hydrocarbons  leaking  from  the  cell  causes  a  fire 
because  the  temperatures  exceed  the  autoignition  temperature.  For  our  cell  design,  we  used  a 
heater  that  totally  encloses  the  cell  to  maintain  uniform  conditions  throughout  the  cell.  The  initial 
design  used  steel  bolts  and  fused  silica  windows.  However,  the  differential  expansion  of  the  bolts 
and  the  windows  resulted  in  leaks  through  the  graphite  foil  seals  at  high  temperatures.  We 
substituted  titanium  bolts  and  glass  windows,  which  have  very  similar  thermal  expansion 
coefficients,  and  the  cell  was  tested  to  500  K  and  2000  psi  simultaneously  without  leaks. 


Thermocouples 


Inlet 


Exit  Windows 

Figure  6.  Design  for  high  pressure/high  temperature  cell. 


Backup  Temperature  Controller 


Figure  7.  Overall  high  pressure  system. 
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Infrared  Stimulated  Scattering  Measurements 

During  construction  of  the  supercritical  cell,  we  received  a  functioning  tunable  diode  laser, 
which  provides  8  mW  at  1064  nm.  We  substituted  this  laser  for  the  ring  dye  laser  and 
reconfigured  the  experiment  to  operate  at  1064  nm  as  shown  in  Figure  8.  Although  the 
experiments  we  had  performed  at  532  nm  greatly  aided  our  work  at  1064  nm,  several  changes 
were  required.  The  8-mW  power  was  less  than  the  approximately  100  mW  we  had  using  the  ring 
dye  laser.  In  addition,  the  silicon  detector  we  used  to  detect  the  probe  signal  had  a  poor  sensitivity 
at  1064  nm. 


Acousto-optic  Optical 


Figure  8.  Apparatus  for  stimulated  Rayleigh  and  Brillouin  scattering  at  1064  nm. 

To  obtain  improved  detection  sensitivity  for  the  1064im  diode  laser,  we  substituted  an 
InGaAs  detector  for  the  silicon  detector  used  at  532  nm.  InGaAs  has  a  larger  capacitance  per  unit 
area  than  silicon.  To  obtain  comparable  time  response,  we  had  to  substitute  a  0.1-mm-diameter 
InGaAs  detector  for  the  2.54-mm-diameter  silicon  detector.  The  silicon  detector  had  been  large 
enough  to  detect  all  the  probe  beam  light  when  placed  directly  behind  the  pinhole  used  to  exclude 
scattered  pump  light,  but  the  InGaAs  detector  was  too  small  to  collect  all  the  probe  light  in  this 
manner  and  too  large  to  act  as  a  spatial  filter  in  place  of  the  pinhole.  To  remedy  this  situation,  we 
replaced  the  pinhole  with  an  optical  fiber,  which  could  then  be  coupled  directly  to  the  InGaAs 
detector.  This  made  alignment  easier  than  it  had  been  using  the  pinhole  and  silicon  detector. 
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Because  of  the  lower  laser  power,  it  was  not  necessary  to  use  a  chopper  wheel  for  measurements  at 
1064  nm. 

The  diode  laser  was  more  sensitive  to  scattered  pump  light  than  the  dye  laser,  and 
amplitude  oscillations  were  produced  even  with  an  optical  isolator  placed  in  front  of  the  diode 
laser.  To  reduce  these  effects,  we  added  a  second  optical  isolator,  added  a  spatial  filter  to  the  probe 
beam,  and  added  apertures  to  the  probe  and  pump  beams  on  either  side  of  the  fluid  cell. 

With  these  modifications,  we  obtained  stimulated  spectra  using  the  tunable  diode  laser.  An 
example  of  a  stimulated  Rayleigh-Brillouin  spectrum  in  hexane  is  shown  in  Figure  9.  This 
spectrum  has  three  peaks — a  double-peaked  Rayleigh  peak  centered  at  the  frequency  of  the  incident 
laser  beam  and  two  Brillouin  peaks  shifted  above  and  below  the  incident  laser  frequency  by  an 
amount  proportional  to  the  speed  of  sound.  To  our  knowledge,  these  are  the  first  stimulated 
scattering  spectra  taken  without  using  a  single-mode  dye  laser.  The  most  significant  improvement 
is  that  the  diode  laser  has  much  less  amplitude  noise  than  the  dye  laser.  Properties  extracted  from 
these  measurements  are  summarized  in  Table  1. 

A  special  feature  of  stimulated  scattering  measurements  at  1064  nm  is  the  large  stimulated 
Rayleigh  peak,  which  is  due  to  the  presence  of  stimulated  thermal  scattering.  Electrostrictive  and 
absorptive  (thermal)  effects  are  two  primary  mechanisms  for  the  coupling  between  the  electric  field 
of  light  and  the  materials  fluctuations  responsible  for  the  scattering — acoustic  waves  for  Brillouin 
scattering  and  non-propagating  entropy  waves  for  Rayleigh  scattering.22  The  electrostrictive  effect 
(the  force  produced  by  the  electric  field  through  the  linear  polarizability)  produces  stimulated 
scattering  in  the  absence  of  absorption.  With  absorption,  thermal  forces  produce  new  forms  of 
stimulated  scattering. 

Liquids  containing  hydrogen  atoms  typically  display  absorption  due  to  overtone  or 
combination  bands  of  the  OH  or  CH  stretch  vibrations  at  1064  nm.  This  absorption  causes 
significant  changes  in  the  stimulated  Rayleigh  and  stimulated  Brillouin  spectra  as  compared  with 
measurements  performed  in  the  visible  spectral  region  (such  as  in  Figure  4),  where  many  liquids 
have  no  absorption.  For  pure  electrostrictive  scattering,  stimulated  Rayleigh  scattering  is  typically 
on  the  order  of  100  times  smaller  than  stimulated  Brillouin  scattering,  and  extremely  difficult  to 
observe.  However,  with  thermal  effects,  the  stimulated  Rayleigh  scattering  can  become  larger  than 
the  stimulated  Brillouin  scattering.  This  allows  the  stimulated  Rayleigh  peak  to  be  easily  observed 
at  1064  nm.  We  have  also  observed  electrostrictive  stimulated  Rayleigh  scattering  at  1064  nm,  and 
the  transition  from  electrostrictive  to  thermal  Rayleigh  scattering,  in  which  case  the  two  types  of 
scattering  can  entirely  cancel. 
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Figure  9.  Stimulated  Rayleigh  and  Brillouin  spectrum  for  hexane. 


Table  1.  Properties  measured  in  hexane  using  stimulated  Rayleigh  and  Brillouin  scattering. 


Brillouin  shift 

2.91  GHz 

Brillouin  width 

151MHz 

Rayleigh  width 

7.9  MHz 

Acoustic  damping  time 

1.06  ns 

Acoustic  velocity 

1125  m/s 

Acoustic  attenuation 
coefficient 

8.4  x  105  nr1 

Effective  mass  diffusion 
coefficient 

7.1  x  10-6  m2/s 

Isentropic  compressibility 

1.2  x  10-9  m2/N 

Thermal  diffusivity 

1.2  x  10-7  m2/s 

We  had  been  concerned  that  scattered  pump  light  would  be  a  significant  problem  for 
measurements  with  the  diode  laser.  The  diode  laser  had  about  10  times  less  power  than  the  dye 
laser,  and  the  gain  coefficient  was  over  4  times  smaller  at  1064  nm  than  at  532  nm.  However,  we 
found  that  the  scattered  light  was  almost  insignificant  for  our  measurements  at  1064  nm,  primarily 
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due  to  the  steps  we  had  taken  to  reduce  scattered  light  for  532  nm,  the  improved  spatial  filtering 
provided  by  the  fiber,  the  lower  noise  of  the  diode  laser  (which  allowed  measurement  of  smaller 
stimulated  gains,  or  smaller  pump  energies),  and  the  generally  decreasing  trend  in  the  amount  of 
scatter  at  longer  wavelengths.  This  is  very  good  news  because  we  expected  scattered  pump  light  to 
be  the  primary  obstacle  to  performing  measurements  in  the  infrared. 

When  operating  the  cell  at  high  temperatures,  we  found  that  convection  near  the  windows 
caused  beam  wander.  This  made  measurements  difficult,  particularly  in  the  region  of  the  critical 
point.  To  avoid  these  effects,  we  produced  a  second-generation  cell.  Having  found  that  scattered 
light  was  less  of  a  problem  at  1  pm,  we  made  a  shorter  cell  with  the  windows  normal  to  the  cell 
rather  than  tilted  as  before.  This  cell  had  a  total  internal  length  of  2.5  cm.  A  small  tilt  of  the  entire 
cell  was  sufficient  to  reduce  scattered  light  to  required  levels.  The  shorter  cell  could  be  more 
evenly  heated  than  the  long  cell,  reducing  convection.  To  further  reduce  convection,  aluminum 
spacers  were  placed  between  the  ends  of  the  cell  and  the  end  of  the  heater,  and  insulating  end  plates 
were  placed  on  the  ends  of  the  oven. 

Figure  10  shows  an  example  i  1  1  r  T  ~T“" 

of  stimulated  scattering  performed  at  II 

200°C  and  630  psi.  Our  _ J\ Signal  L _ 

measurements  are  the  first  stimulated  l|  \  ( 

Rayleigh  and  stimulated  Brillouin  II  1[ 

spectra  taken  at  supercritical  « _ J  _ L _ 

conditions  to  our  knowledge.  Note  m  .  1  \j 

that  we  are  still  able  to  obtain  3  \  y 

measurements  at  a  high  signal-to-  _ J  \ _ _ _ _ 

noise  ratio  at  these  conditions.  Also  \ 

shown  in  Figure  10  are  fits  of  || 

theoretical  lineshapes  to  the  data  and  ^Thermal  Lineshapes  A _ ^ _ 

the  electrostrictive  and  thermal  • 

portions  of  those  fits.  Note  that  the  _15  .1>0  _0  5  0  0  0.5  i.o  1.5 

signal  data  show  an  asymmetry  in  the  Frequency  (GHz) 

„  .  ,  .  ,  .  Figure  10.  Stimulated  Rayleigh  and  Brillouin  spectrum  for 

Brillouin  (outer)  peaks.  Thts  ^  at  2C0.C  and  630  psj 

asymmetry  is  produced  by  the 

imaginary  Lorentzian  lineshape  of  the  thermal  stimulated  Brillouin  peak.  Stimulated  thermal 
Brillouin  and  Rayleigh  scattering  is  caused  by  the  moderate  absorption  of  the  Nd:  YAG 
fundamental  wavelength  in  organic  liquids.  Without  the  thermal  effects,  the  stimulated  Rayleigh 
peak  would  be  extremely  small  and  difficult  to  measure. 


Thermal  Lineshapes 


Frequency  (GHz) 

Figure  10.  Stimulated  Rayleigh  and  Brillouin  spectrum  for 
hexane  at  200°C  and  630  psi. 
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We  have  shown  that  stimulated  Brillouin  and  Rayleigh  scattering  is  a  robust  and  precise 
technique  for  performing  measurements  in  supercritical  fluids.  This  work  forms  an  excellent 
foundation  for  future  studies  of  supercritical  fluids. 
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ABSTRACT 

We  have  demonstrated  stimulated  Rayleigh  /  Brillouin  scattering  at  a  wavelength  of  1  pm 
using  an  injection  seeded  Nd:YAG  laser  as  a  pump  laser  and  a  tunable  diode  laser  as  a  probe 
laser.  Spectra  with  good  signal-to-noise  ratio  are  obtained  in  spite  of  the  low  probe  beam  power 
and  small  gain  coefficient  in  the  infrared.  Stimulated  Rayleigh  scattering  is  readily  observable  in 
most  liquids  due  to  absorption  by  the  OH  and  CH  overtone  bands.  The  absorption  also  causes  an 
asymmetry  in  the  stimulated  Brillouin  peak.  A  Rayleigh  linewidth  of  8  MHz  is  measured  with 
this  approach. 
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Stimulated  scattering  may  be  used  to  perform  noninvasive  in  situ  measurements  of 
material  physical  properties.  Stimulated  Rayleigh  and  stimulated  Brillouin  scattering  provide 
information  on  the  thermal  and  elastic  properties,  respectively,  of  the  material  under  study. 
Linewidths  for  stimulated  Rayleigh  and  stimulated  Brillouin  scattering  are  typically  on  the  order 
of  10  and  100  of  MHz,  respectively.  Thus,  narrowband  lasers  are  required  for  Rayleigh  and 
Brillouin  scattering  measurements. 

These  measurements  are  non-resonant,  and  may  be  performed  at  any  wavelength  for 
which  the  medium  under  study  is  transparent  or  nearly  transparent.  Previous  experiments  on 
stimulated  Rayleigh  and  Brillouin  scattering  have  used  a  ring  dye  laser  as  either  a  pump  or  probe 
laser.1'3  As  an  alternative  to  ring  dye  lasers,  we  have  used  an  external  cavity  laser  diode 
operating  at  1064  nm.  Although  this  laser  provides  over  an  order  of  magnitude  less  power  than  a 
ring  dye  laser  and  the  gain  coefficient  is  more  than  4  times  smaller  at  1064  nm  than  at  532  nm, 
we  find  very  good  results  with  the  laser  diode.  In  addition,  the  absorption  present  in  most  liquids 
at  1064  nm  makes  stimulated  Rayleigh  scattering  readily  observable. 

Liquids  containing  hydrogen  atoms  typically  display  absorption  due  to  overtone  or 
combination  bands  of  the  OH  or  CH  stretch  vibrations  at  our  measurement  wavelength  of  1064 
nm.  This  absorption  causes  significant  changes  in  the  stimulated  Rayleigh  and  stimulated 
Brillouin  spectra  as  compared  with  measurements  performed  in  the  visible  spectral  region,  where 
many  liquids  have  little  or  no  absorption. 

Brillouin  and  Rayleigh  scattering  involve  scattering  of  light  from  density  fluctuations 
caused  by  acoustic  waves  and  non-propagating  entropy  waves,  respectively.  For  stimulated 
Rayleigh  and  Brillouin  scattering,  there  are  two  primary  mechanisms  for  coupling  the  electric 
field  of  the  light  and  the  material  fluctuations-electrostriction  and  absorption.4  The 
electrostrictive  effect  (the  force  produced  by  the  electric  field  through  the  linear  polarizability) 
produces  stimulated  scattering  in  the  absence  of  absorption.  With  absorption,  thermal  forces  also 
can  cause  material  fluctuations.  For  pure  electrostrictive  scattering,  stimulated  Rayleigh 
scattering  is  typically  on  the  order  of  100  times  smaller  than  stimulated  Brillouin  scattering,  and 
extremely  difficult  to  observe.  However,  with  thermal  effects,  the  stimulated  Rayleigh  scattering 
can  become  larger  than  the  stimulated  Brillouin  scattering,  and  is  easily  observed. 

Our  experimental  apparatus  for  stimulated  scattering  at  1  pm  is  shown  in  Fig.  1.  The 
pump  laser  is  an  injection-seeded  Nd:YAG  laser.  This  laser  is  an  attractive  pump  laser  for 
stimulated  scattering  measurements  as  it  can  provide  high  peak  powers  to  enhance  the  nonlinear 
response  and  near-Fourier  transform  limited  linewidths  to  provide  high  spectral  resolution.  By 
reducing  the  oscillator  pump  energy,  Fourier  transform-limited  spectral  linewidths  as  narrow  as 
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20  MHz  have  been  obtained  with  this  laser.5  The  probe  laser  is  a  diode  laser,  tuned  through 
feedback  from  an  external  grating  (EOSI  model  2010).  This  laser  has  not  been  actively 
stabilized,  so  we  rely  on  passive  stability  for  our  measurements.  The  elongated  beam  profile  of 
the  diode  laser  is  reduced  and  circularized  using  an  anamorphic  prism  pair  (Melles  Griot 
06  GPA  003)  producing  a  beam  of  roughly  1  mm  in  diameter.  An  acousto-optic  modulator 
directs  the  probe  beam  to  a  Fizeau  wavelength  meter  (New  Focus  model  7711),  except  for  a  15 
ms  window  centered  at  the  time  of  pump  pulse,  when  the  modulator  is  turned  off  and  the  probe 
beam  passes  to  the  sample.  Both  the  pump  and  probe  beams  are  spatially  filtered  to  remove 
higher  order  transverse  mode  structure.  The  pump  and  probe  beams  are  overlapped  in  a  cell 
containing  the  liquid  in  a  counterpropagating  geometry  with  a  crossing  angle  of  approximately 
25  mrad.  The  gain  signal  is  detected  with  a  InGaAs  photodiode  (Fermionics  FD100FC)  and  an 
transimpedance  amplifier. 

The  Nd:YAG  laser  has  approximately  107  times  more  power  than  the  diode  laser  during  a 
typical  measurement.  To  prevent  scattered  Nd:YAG  light  from  influencing  the  diode  laser,  two 
optical  isolators  are  used  on  the  diode  laser  beam.  To  prevent  scattered  pump  laser  light  from 
overwhelming  the  gain  signal,  the  detected  light  is  spatially  filtered  using  an  iris,  a  lens  and  a 
62.5  pm  core  optical  fiber  coupled  to  the  InGaAs  detector. 

The  wavelength  meter  is  calibrated  between  400  and  1000  nm  only.  For  these 
measurements  we  require  accurate  relative  wavelength  changes,  not  absolute  wavelength  values, 
we  tested  the  accuracy  of  the  relative  wavelength  readings  by  comparing  the  wavelength 
readings  with  the  heterodyne  frequency  obtained  by  beating  the  diode  laser  and  seed  laser  beams. 
The  two  measurements  agreed  to  within  1.5%. 

A  stimulated  Rayleigh/Brillouin  spectrum  for  hexane  obtained  with  this  system  is  shown 
as  the  top  curve  in  Fig.  2.  The  abscissa  is  the  difference  in  frequency  between  the  pump  and 
probe  lasers.  The  outer  two  peaks  are  Brillouin  peaks;  the  central  pair  of  peaks  is  due  to 
stimulated  Rayleigh  scattering.  The  positive-going  Brillouin  peak  to  the  left  is  a  gain  peak, 
corresponding  to  transfer  of  power  from  the  Nd:  YAG  laser  to  the  probe  laser.  The  negative¬ 
going  peak  to  the  right  is  a  loss  peak,  wherein  power  is  transferred  from  the  probe  laser  to  the 
pump  laser. 

We  performed  fits  to  the  measured  lineshapes  to  determine  linewidths,  lineshifts,  and 
peak  heights.  The  electrostrictive  Brillouin  lineshape  is  described  by  the  real  part  of  the  complex 
Lorentzian  profile;  the  thermal  Brillouin  and  Rayleigh  lineshapes  are  described  by  the  imaginary 
parts  of  a  complex  Lorentzian  profile.4  The  measured  lineshapes  are  described  by  the 
convolution  of  the  Gaussian  spectral  lineshape  of  the  Nd:YAG  pump  laser  with  these  Lorentzian 
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profiles.  These  convolved  lineshape  may  be  expressed  as  the  real  and  imaginary  parts  of  the 
complex  error  function  for  the  Gaussian  for  the  Brillouin  and  Rayleigh  peaks,  respectively  (the 
former  is  a  Voigt  profile).  Fits  were  performed  using  the  algorithm  of  Humlicek,6,7  and  the  Igor 
graphics  program  (Wavemetrics).  We  use  the  Fourier-transform  limited  linewidth  calculated 
from  pulsewidth  as  the  spectral  linewidth. 

The  full  fit  to  the  Rayleigh/Brillouin  spectrum  is  shown  as  the  curve  labeled  "fit"  in  Fig.  2 
Very  good  agreement  is  found  between  the  measured  and  fit  lineshapes  (see  also  Figs.  3  and  4). 
The  electrostrictive  and  thermal  contributions  to  the  overall  fit  are  shown  separately  at  the 
bottom  of  Fig.  2.  Note  that  the  imaginary  Lorentzian  shape  of  the  thermal  Brillouin  scattering 
creates  an  asymmetry  for  the  Brillouin  peaks.  This  is  seen  more  clearly  in  Fig.  3,  which  shows 
the  Brillouin  gain  peak  from  Fig.  2.  The  spectroscopic  and  thermodynamic  properties  for 
n-hexane  determined  from  numerous  such  fits  are  given  in  Table  1. 

A  scan  over  the  stimulated  Rayleigh  peak  is  shown  in  Fig.  4.  As  the  wavemeter  does  not 
provide  sufficient  spectral  resolution  over  such  a  short  frequency  range,  the  frequency  values  for 
this  measurement  were  determined  from  the  heterodyne  signal  obtained  by  combining  a  portion 
the  seed  and  probe  lasers  on  a  photodiode.  This  peak  is  quite  narrow.  The  average  linewidth 
determined  from  multiple  fits  such  as  shown  in  Fig.  4  gives  value  of  7.9±3.2  MHz,  which  is  in 
good  agreement  with  the  value  expected  from  the  thermal  diffusivity  for  n-hexane.  These 
measurements  were  performed  with  a  pulsewidth  for  the  Nd:YAG  laser  of  approximately  30  ns, 
which  corresponds  to  a  full  width  half  maximum  spectral  linewidth  of  about  15  MHz.  Thus  the 
laser  linewidth  is  about  twice  the  measured  linewidth  of  the  Rayleigh  peak.  The  good  agreement 
with  the  expected  linewidth  confirms  that  the  Nd:YAG  laser  has  a  transform-limited  linewidth  of 
15  MHz. 

For  comparison,  a  stimulated  scattering  spectrum  of  methanol  taken  at  532  nm  is  shown 
in  Fig.  5.  This  spectrum  was  acquired  using  the  second  harmonic  of  the  same  Nd:YAG  laser  and 
a  cw  ring  dye  laser  (Coherent  CR-699)  as  a  probe  laser.  The  spectra  taken  at  1064  (Fig.  2)  and 
532  nm  (Fig.  5)  are  of  similar  quality,  except  that  the  infrared  measurement  has  the  strong 
stimulated  thermal  Rayleigh  peak. 

The  ability  to  perform  stimulated  scattering  in  the  infrared  opens  new  measurement 
possibilities  for  materials  that  are  only  transparent  in  the  infrared,  allowing  measurements  within 
the  bulk  material  rather  than  only  on  the  surface.  Measurement  of  spontaneous  Brillouin  or 
Rayleigh  scattering  in  the  infrared  are  greatly  hampered  by  smaller  scattering  cross  sections  and 
poorer  detector  performance.  However,  with  stimulated  scattering  the  signal  powers  can  be 
nearly  as  large  as  the  laser  power,  overcoming  these  limitations.  This  allows  bulk  measurement 
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in  materials  that  are  not  transparent  in  the  visible  such  as  most  semiconductors  and  liquids 
containing  multi-ring  aromatic  compounds. 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research  under  contract  F- 
49620-97 -C-0002 . 
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Table  1. 

Properties  measured  in  hexane  using  stimulated  Rayleigh  /Brillouin  scattering. 


Brillouin  shift 

2.91±0.06  GHz 

Brillouin  width 

151±8  MHz 

Rayleigh  width 

7.9±3.2  MHz 

Acoustic  damping  time 

1.06  ns 

Acoustic  velocity 

1125  m/s 

Acoustic  attenuation 
coefficient 

8.4  x  105  m-1 

Effective  mass  diffusion 
coefficient 

7.1  x  10'6  m2/s 

Isentropic  compressibility 

1.2  x  10'9  m2/N 

Thermal  diffusivity 

1.2  x  10-7  m2/s 
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Fig.  1.  Experimental  apparatus  for  stimulated  Brillouin  and  stimulated  Rayleigh  scattering. 
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Fig.  2.  Stimulated  Rayleigh/Brillouin  spectrum  for  n-hexane  (top),  fit  to  the  data  (middle)  and 
electrostrictive  and  thermal  components  to  the  fit  (bottom).  The  dashed  lines  indicate  zero  gain  level  for 

each  curve. 
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Gain  Signal  (arbitrary  units) 


Fig.  3.  Data  (points)  and  fit  (solid  line)  for  stimulated  Brillouin  gain  peak  for  n-hexane.  Stimulated  thermal 

Brillouin  scattering  gives  an  asymmetry  to  the  peak. 
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Fig.  5.  Stimulated  Brillouin  spectrum  for  methanol  measured  at  532  nm. 
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APPENDIX  B 


ELECTROSTRICTIVE  AND  THERMAL  STIMULATED  RAYLEIGH 
SPECTROSCOPY  IN  LIQUIDS,  MANUSCRIPT  IN  PREPARATION. 


Rayleigh  scattering  involves  the  scattering  of  light  by  non-propagating  isobaric  entropy 
fluctuations.  Density  variations  associated  with  these  fluctuations  cause  a  modulation  of  the 
refractive  index,  and  the  spatial  variations  of  the  refractive  index  result  in  a  scattering  of  light. 

When  the  incident  light  is  sufficiently  intense,  stimulated  Rayleigh  scattering  occurs.  In  this  case, 
coupling  of  the  incident  (pump)  light  and  the  fluctuations  can  produce  amplification  of  the  scattered 
light,  and  coupling  between  the  incident  and  scattered  light  fields  can  amplify  the  fluctuations. 

In  this  Letter  we  describe  the  use  of  tunable  Rayleigh  gain  spectroscopy  in  electrostrictive 
and  predominantly  absorptive  liquids.  This  is  the  first  observation  of  the  transition  from 
electrostrictive  to  thermal  stimulated  Rayleigh  scattering  to  our  knowledge.  This  transition  is  * 
produced  through  varying  the  medium  absorption.  In  agreement  with  theory,1  for  a  certain  amount 
of  absorption  the  electrostrictive  and  the  absorptive  part  of  the  gain  factor  cancel  (see  Eq.  1). 

Two  different  mechanisms  can  produce  the  coupling  between  light  and  medium  for 
simulated  Rayleigh  scattering — electrostriction  and  light  absorption.  In  any  medium,  the  electric 
field  of  the  electromagnetic  light  wave  can  create  a  dipole  moment  through  the  medium 
polarizability.  When  the  electric  field  is  not  uniform,  as  is  always  the  case  for  an  electromagnetic 
wave,  the  energy  associated  with  the  induced  dipole  moment  and  the  electric  field  will  vary 
spatially.  The  atoms  or  molecules  in  the  field  will  experience  a  force  associated  with  reaching  the 
lowest  energy  state.  This  force  is  called  electrostriction.  In  absorptive  media,  temperature  changes 
due  to  the  absorption  of  light  result  in  additional  density  fluctuations,  which  give  rise  to  stimulated 
thermal  Rayleigh  scattering. 


Electrostrictive  and  thermal  effects  tend  to  produce  opposite  influences  on  density. 
Electrostrictive  effects  generally  cause  the  density  to  increase  in  the  regions  with  higher  electric 
field.  Thermal  effects  cause  a  lower  density  in  the  region  of  the  higher  electric  field  as  the  heated 
medium  expands.  The  electrostrictive  effect  is  relatively  weak,  and  it  is  quite  difficult  to  observe 
pure  electrostrictive  stimulated  Rayleigh  scattering.  Stimulated  thermal  Rayleigh  scattering  is  more 
readily  observed  because  the  absorptive  gain  factors  are  large  compared  to  the  electrostrictive  ones. 

The  amplification  or  gain  for  the  scattered  light  produced  by  a  strong  pump  light  beam  can 
be  monitored  using  a  second  probe  light  beam.  The  spectral  variation  in  gain  for  the  probe  beam  is 
given  by  the  imaginary  part  of  a  complex  Lorentzian  function:1 


gRL  =  [geRL(max)  -  gaRL(max)] 


4(o/Trl 

1  +  (2o)/rRL)2  ' 


(1) 


where  Trl  is  the  spontaneous  Rayleigh  linewidth,  geRL(max)  and  gaRL(max)  are  the  maximum 
values  for  the  electrostrictive  and  the  absorptive  gain  factors,  respectively,  and  to=a)p-(Os  is  the 
difference  between  the  frequencies  top  of  the  pump  beam  and  tos  of  the  probe  (or  Stokes)  beam. 
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We  use  a  pump  and  probe  experiment  to  measure  stimulated  Rayleigh  scattering.  The  pump 
and  probe  beams  are  crossed  in  the  sample  cell  containing  the  liquid.  The  cw  probe  beam  is  used  to 
measure  the  gain  produced  by  the  intense  pulsed  pump  light.  The  maximum  gain  signal  is  obtained 
for  counterpropagating  pumping. 

The  experimental  apparatus  used  to  measure  the  frequency  dependence  of  the  gain  factor  is 
illustrated  in  Fig.  1.  Ray  trajectories  are  represented  by  dashed  lines.  We  use  the  light  from  an 
injection-seeded,  home-built,  single-mode  Nd:YAG  laser2  at  1064  nmto  provide  the  pump 
radiation.  For  our  experiment,  it  is  operated  with  a  pulse  duration  of  about  29  ns  (FWHM)  and  a 
repetition  rate  of  10  Hz.  The  spectral  bandwidth  of  the  laser  is  Fourier  transform  limited,  providing 
a  linewidth  of  about  15  MHz  for  a  29  ns  pulse.  The  probe  light  is  provided  by  a  tunable  single¬ 
mode  diode  laser  with  a  linewidth  smaller  than  300  kHz  in  50  ms.  An  acousto-optic  modulator 
directs  the  probe  beam  to  a  wavemeter  in  the  intervals  between  the  pump  pulses  for  determination 
of  diode  laser  wavelength. 

The  spatial  filters  serve  to  suppress  higher  order  transverse  spatial  modes.  Furthermore, 
the  spatial  filter  in  the  probe  trajectory,  together  with  the  two  optical  isolators,  prevents  pump  light 
from  reaching  the  resonator  of  the  tunable  laser,  which  would  otherwise  cause  frequency  and 
amplitude  fluctuations  on  the  probe  beam.  Inside  the  cell,  the  counterpropagating  (180°  crossing 
angle)  pump  and  probe  beams  are  overlapped  at  their  focal  points  to  obtain  a  large  gain  signal. 
Detection  is  performed  using  an  InGaAs  photodiode  with  a  transimpedance  amplifier.  The  output 
signal  is  split  with  a  bias  tee  into  a  high  frequency  part,  which  is  the  gain  signal  for  the  ~29  ns 
pulse,  and  a  low  frequency  part,  which  monitors  the  probe  beam  power.  The  energy  of  the 
Nd:YAG  beam  is  measured  with  a  pyroelectric  energy  detector  connected  to  an  energy  meter.  The 
measured  gain  signals  are  divided  by  the  probe  and  pump  beam  powers  to  compensate  for 
fluctuations  in  the  laser  powers.  A  fast  silicon  detector  is  used  to  measure  the  pulse  duration  of  the 
pump  beam. 

For  Rayleigh  gain  spectroscopy  in  liquids,  the  resolution  of  the  wavemeter  is  only  slightly 
smaller  than  the  scanning  interval  of  about  100  MHz.  To  allow  high-resolution  measurements  of 
stimulated  Rayleigh  scattering,  we  directly  measure  the  frequency  difference  to  between  the  pump 
and  the  signal  beams.  Glass  plates  are  employed  to  separate  a  few  percent  of  the  beams  from  the 
diode  laser  and  seed  laser  used  to  injection  seed  the  pump  laser.  These  beams  are  combined  on  a 
beamsplitter  and  focused  on  a  detector.  The  detector  measures  the  frequency  difference  between 
the  two  lasers;  the  sum  frequency  is  far  beyond  the  detection  bandwidth.  This  heterodyne  or  beat 
signal  provides  a  direct  measure  of  the  frequency  difference  between  the  lasers,  including  a  well 
defined  zero  point  of  the  frequency  shift.  Moreover,  frequency  fluctuations  of  the  seed  laser  are 
now  taken  into  account. 
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We  use  freon  113  (1,1,2-trichlorotrifluoroethane,  anhydrous)  to  observe  stimulated 
electrostrictive  Rayleigh  scattering.  Most  liquids  containing  hydrogen  atoms  cannot  be  used 
because  they  have  absorption  at  1064  nm  due  to  overtone  and  combination  vibrational  absorption 
from  O-H  and  C-H  bonds. 

The  pure  electrostrictive  stimulated  Rayleigh  peak  is  small.  We  used  pump  beam  energies 
of  about  5  mJ  in  a  30  ns  pulse  to  observe  the  electrostrictive  signal  in  freon  113.  This  corresponds 
to  a  power  of  about  10^  times  larger  than  the  ~1  mW  of  probe  beam  power  available  at  the  sample 
cell.  The  gain  signals  were  approximately  1%  of  the  probe  beam  power.  Thus  even  a  small 
fraction  of  scattered  pump  light  can  completely  overwhelm  the  electrostrictive  gain  signal.  This 
effect  is  exacerbated  by  the  fact  that  the  pump  and  probe  beam  have  nearly  identical  frequencies  in 
the  region  of  the  Rayleigh  peak.  Thus,  a  heterodyne  signal  between  the  probe  beam  and  scattered 
pump  light  is  also  present  at  the  detector  used  to  measure  the  gain  signal.  This  heterodyne  process 
provides  an  effective  amplification  of  the  scattered  pump  light  as  the  detected  heterodyne  signal  is 
proportional  to  the  square  root  of  the  product  of  the  pump  and  probe  powers.  Thus  the  1%  gain 
signal  is  only  as  big  as  the  heterodyne  signal  from  a  scattered  pump  light  power  of  KT4  of  the 
probe  beam  power,  or  10' 1 1  of  the  pump  beam  power. 

The  influence  of  scattered  pump  light  was  reduced  by  using  very  clean  windows  and  freon, 
by  tilting  the  sample  cell  to  direct  reflections  from  the  windows  away  from  the  detector,  blocking  all 
stray  reflections  and  other  scattered  light  with  black  screens,  making  the  pump  and  probe  beams  of 
approximately  the  same  diameter,  spatially  filtering  the  pump  and  probe  beams  to  remove  higher 
order  transverse  spatial  modes,  and  spatially  filtering  the  probe  beam  before  the  detector  using  an 
iris,  lens,  and  multimode  fiber.  However,  we  could  not  reduce  the  scattered  light  below  the  level  * 

from  the  spontaneous  Rayleigh  scattering  of  the  freon  in  the  volume  determined  by  the  overlap 
between  the  pump  and  probe  beams.  This  spontaneous  Rayleigh  scattering  is  the  dominant  noise 
source  for  our  measurements.  We  reduce  the  influence  of  the  heterodyne  signal  from  spontaneous 
Rayleigh  scattering  by  averaging  the  measured  gain  over  many  laser  shots  (typically  60  shots  for 
the  stimulated  Rayleigh  measurements). 

We  determined  the  height  of  the  Rayleigh  peak  by  fitting  an  appropriate  lineshape  to  the 
measured  gain  signal.  The  spectral  dependence  of  the  gain  signal  is  described  by  the  convolution 
of  the  stimulated  Rayleigh  lineshape  {imaginary  part  of  the  complex  Lorentzian,  Eq.  (1)}  and  the 
Gaussian  spectral  profile  of  the  pulse  laser.  This  convolution  may  be  expressed  in  terms  of  the 
imaginary  part  of  the  complex  error  function.  We  used  the  approximation  of  Humlicek3,4  for  the 
complex  error  function  to  perform  the  fitting.  The  linewidth  of  the  Fourier-transform  limited 
Nd:YAG  pulse  was  determined  by  measuring  the  temporal  duration  of  the  pulse.  The  contribution 
of  the  probe  laser  linewidth  is  negligible. 
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The  solid  curve  at  the  top  of  Fig.  2  shows  the  measured  electrostrictive  Rayleigh  gain  signal 
in  freon  113.  The  width  of  the  Rayleigh  peak  is  broadened  considerably  by  the  spectral  linewidth 
of  the  pulse  laser.  The  dashed  curve  shows  the  result  of  the  curve  fitting.  Freon  113  does  not 
absorb  at  the  Nd:  YAG  wavelength  of  1064  nm  and  shows  a  relatively  high  electrostrictive  gain.  To 
produce  some  absorption  in  the  freon,  we  use  ethanol,  which  mixes  with  freon  and  has  moderate 
absorption  at  1064  nm.  After  adding  0.10%  (volume  percentage)  of  ethanol,  the  electrostrictive 
and  the  absorptive  gain  have  the  same  value  and  cancel  each  other  out  (see  Eq.  1).  The  measured 
curve  is  shown  in  the  middle  of  Fig.  2.  The  Rayleigh  gain  signal  (solid  curve)  shows  no  Rayleigh 
peak,  and  the  curve  fitting  (dashed  curve)  produces  a  flat  line.  The  small  variations  in  the  measured 
gain  signal  are  due  to  the  heterodyne  signal  from  scattered  pump  light.  After  doubling  the  amount  r 
of  ethanol  in  freon  1 13,  we  get  a  gain  signal  that  has  about  the  same  height  as  in  pure  freon,  but  is 
inverted.  The  measured  gain  signal  for  0.19%  of  ethanol  is  shown  at  the  bottom  of  Fig.  2  (solid 
curve),  the  dashed  curve  is  the  fitted  profile. 

We  have  measured  the  absorption  coefficient  for  different  volume  percentages  of  ethanol  in 
freon  (see  Fig.  3).  The  absorption  coefficient  displays  a  linear  dependence  on  the  percentage  of 
ethanol.  This  measurement  shows  that  a  very  small  absorption  coefficient  (0.00012  cm-1)  is 
sufficient  to  cause  the  cancellation  of  the  electrostrictive  Rayleigh  scattering  at  0.1%  ethanol. 

For  small  amounts  of  ethanol  in  freon  the  physical  properties  (apart  from  the  absorption)  are 
changed  only  little.  Theory1  predicts  a  linear  growth  of  the  maximum  absorptive  gain  factor 
gaRL(max)  with  the  absorption  coefficient,  whereas  the  maximum  electrostrictive  gain  factor 
geRL(max)  does  not  depend  on  the  absorption.  Measurements  of  the  gain  signal  were  performed  for 
various  amounts  of  ethanol  in  freon,  and  the  height  of  the  Rayleigh  peak  was  determined.  The 
results  are  shown  in  Fig.  4  for  small  percentages  of  ethanol.  We  display  peak  height  rather  than 
peak  area  as  the  area  of  the  Rayleigh  peak  is  infinite  {Eq.  (1)  behaves  as  1/to  for  large  to}.  As 
expected,  the  height  displays  a  linear  dependence  on  the  absorption  coefficient  {or  geRL.(max)  - 
gaRL(max)}.  The  height  has  a  negative  value  for  gV/max^O  in  pure  freon,  becomes  zero  for 
gaRL(max)=geRL(rnax),  corresponding  to  0.1%  of  ethanol,  and  has  positive  values  for 

gV(max)>geRL(max). 

When  the  probe  laser  is  scanned  over  a  larger  frequency  range,  stimulated  Brillouin 
scattering  is  observed.  The  shapes  of  the  stimulated  electrostrictive  and  thermal  Brillouin  scattering 
peaks  are  described  by  the  real  and  imaginary  parts  of  the  complex  Lorentzian  profile, 
respectively.1  Only  stimulated  electrostrictive  Brillouin  scattering  is  observed  in  pure  freon  1 13 
(Fig.  5).  Gain  and  loss  peaks  are  observed  at  the  left  and  right  of  Fig.  5,  respectively.  The  gain 
and  loss  peaks  corresponds  to  transfer  of  energy  from  the  pump  laser  to  the  probe  laser  and  vice 
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versa.  The  stimulated  electrostrictive  Rayleigh  scattering,  which  is  typically  about  100  times 
smaller  than  the  stimulated  electrostrictive  Brillouin  scattering,  is  too  small  to  be  seen  in  Fig.  5. 

The  modest  absorption  coefficient  of  ethanol  at  1064  nm  (0.12  cm-1)  is  sufficient  to  cause 
very  large  stimulated  thermal  Rayleigh  scattering  in  ethanol  (Fig.  6).  The  stimulated  thermal 
Rayleigh  peak  is  much  larger  than  the  stimulated  electrostrictive  Brillouin  scattering  in  ethanol. 
Stimulated  thermal  Brillouin  scattering  causes  asymmetry  of  the  Brillouin  peaks  in  this  case. 

We  have  reported  an  experimental  investigation  of  stimulated  Rayleigh  scattering  in  liquids 
by  tunable  Rayleigh  gain  spectroscopy.  We  have  observed  the  electrostrictive  Rayleigh  peak  in 
freon  1 13,  which  is  the  first  observation  of  stimulated  electrostrictive  Rayleigh  scattering  in  a  liquid 
to  our  knowledge.  By  adding  small  amounts  of  ethanol,  we  have  increased  the  absorption 
coefficient  to  the  point  where  electrostrictive  and  absorptive  gain  cancel  out.  For  a  further  increase 
of  the  absorption,  the  peak  inverts  from  its  original  form.  In  agreement  with  theory,  a  linear 
relationship  between  the  height  of  the  peak  and  the  absorption  coefficient  has  been  observed. 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research  under  contract  F- 
49620-97 -C-0002. 
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Fig.  2.  Measured  stimulated  Rayleigh  gain  signal  (solid  curve)  and  fitted  profile  (dashed  curve)  in  pure 
freon  113  (top),  and  freon  with  0.1%  ethanol  (middle)  and  0.19%  ethanol  (bottom)  as  a  function  of  the 
frequency  difference  between  the  pump  and  probe  lasers. 
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Fig.  6.  Measured  stimulated  Brillouin  gain  signal  in  ethanol  as  a  function  of  the  frequency  difference 
between  the  pump  and  probe  lasers.  Absorption  produces  a  large  stimulated  thermal  Rayleigh  peak  and 
stimulated  thermal  Brillouin  scattering,  which  causes  asymmetry  of  the  stimulated  Brillouin  peaks. 
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ABSTRACT 

Two-photon-resonant  difference  frequency  generation  using  an  ArF  excimer  laser  provides 
widely  tunable  VUV  radiation  at  high  pulse  energies.  Two-photon  resonances  in  H2,  Kr,  and  HD 
are  within  the  tuning  range  of  the  ArF  laser.  With  this  technique,  we  have  directly  measured  over 
65  pJ  at  133  nm.  H2  has  a  significantly  smaller  phase  mismatch  that  Kr,  leading  to  more  efficient 
VUV  generation,  particularly  at  shorter  VUV  wavelengths.  However,  mixing  in  H2  also  produces 
additional  VUV  lines  from  two-photon  excited  amplified  spontaneous  emission  (ASE).  We  have 
observed  new  ASE  lines  produced  in  this  manner.  H2  is  ineffective  at  generation  of  Lyman-a 
radiation  due  to  the  production  of  H  atoms.  With  a  phase  matched  mixture  of  Kr  and  Ar,  we  have 
directly  measured  7  pi  at  Lyman-a.  A  physical  basis  for  the  asymmetric  tuning  profile  in  this  gas 
mixture  is  presented.  With  light  from  this  VUV  source,  we  have  performed  1+1  REMPI  in  Xe  at 
147  nm  and  two-photon-excited  fluorescence  in  neon  at  133  nm. 
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INTRODUCTION 


High  power  VUV  radiation  is  required  for  a  variety  of  applications  such  as  multiphoton 
spectroscopy  and  diagnostics,1-4  plasma  diagnostics,5'11  and  photochemistry.12'15  Approaches 
to  production  of  high  power  VUV  radiation  include  resonantly-enhanced  four-wave  mixing  in 
gases  and  metal  vapors,16  211  high  order  Raman  shifting,2,5,11’ 11,2 1-25  nonlinear  mixing  in  crystals 
(down  to  -170  nm), 24,25  and  direct  laser  sources.26,27  Of  the  broadly  tunable  sources,  Raman 
shifting  offers  simplicity  through  the  use  of  a  single  laser  source,  but  pronounced  amplitude 
variations  are  produced  for  the  large  number  of  anti-Stoked  shifts  required  to  reach  deep  into  the 
VUV.  The  resonantly-enhanced  mixing  schemes  can  provide  high  powers  deep  into  the  VUV,  but 
typically  require  two  lasers.  Two-photon  resonant  sum-  and  difference-frequency  mixing  may  be 
performed  using  resonances  in  Kr, 28,29  Hg,1,30  Xe31,32  and  other  gases33'35  using  frequency- 
doubled  dye  lasers  to  provide  the  two-photon-resonant  photons. 

Although  the  ArF  excimer  laser  has  a  relatively  narrow  tuning  range  of  approximately  1 
nm,  there  are  by  coincidence  a  number  of  two-photon  resonances  available  in  Kr,  H2,  and  HD 
within  this  range.3,36  The  ArF  laser  provides  a  number  of  advantages  for  two-photon  resonant 
difference-frequency  mixing.  The  generated  VUV  power  nominally  scales  as  the  square  of  the 
power  for  the  two-photon  resonant  wavelength,  for  which  the  ArF  can  provide  very  high  energies 
(up  to  100  mJ).  The  short  wavelength  of  the  two-photon-resonance  allows  tuning  to  short  VUV 
wavelengths  (the  nominal  tuning  range  is  1 10  to  180  nm)  and  reduces  the  phase  mismatch. 

The  application  of  ArF  radiation  to  tunable  VUV  generation  through  two-photon-resonant 
difference  frequency  mixing  was  proposed  as  early  as  1984  by  Hilbig  et  al.  (in  Kr),37  reported  in 
1991  by  Strauss  and  Funk38  and  Faris  and  Dyer,39  and  extended  to  femtosecond  VUV  generation 
in  1996  by  Kittelmann  et  al.40  The  ArF  excimer  laser  has  also  been  used  for  VUV  and  XUV 
generation  through  third  harmonic  generation,41  sum  frequency  mixing,42,43  two-photon  pumped 
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amplified  spontaneous  emission,44  Raman  shifting,3,2 1,45,46  and  other  nonlinear  processes.47'49 
All  these  processes,  with  the  exception  of  the  Raman  shifting,  use  two-photon  resonances  to 
provide  enhancement. 

This  paper  describes  the  use  of  an  ArF  laser  for  two-photon-resonant  difference  frequency 
mixing.  Mixing  is  performed  in  Kr  using  the  4p56p[3/2,2]  4p6  ]S  two-photon  resonance  at 

103363.4  cm'1  and  in  H2  using  the  E,F  %+  (v'=  6,J’=1)  X  1Ig+  (v"=0,J"=l)  two-photon 
resonance  at  103487.1  cm*1.  Factors  affecting  conversion  efficiency  and  tuning  range  are 
discussed  in  Section  2.  Experimental  results  for  VUV  generation  are  described  in  Section  3. 
Application  of  this  VUV  radiation  to  multi-photon  spectroscopy  is  presented  in  Section  4. 


2.  VUV  CONVERSION  EFFICIENCY 


For  the  four-wave  sum  and  difference-frequency  mixing  processes  described  by  0)4  = 

2(Di±C02,  the  generated  power  P4  scales  according  to 

P4  °c  N2  [X(3)]2  [Pj]2  p2  F  (1) 


where  N  is  the  number  density  of  atoms  or  molecules,  is  the  third  order  nonlinear 

susceptibility  for  this  mixing  process,  Pi  and  P2  are  the  incident  powers,  and  F  is  a  phasematching 

integral.  For  mixing  using  collinear  Gaussian  beams  with  a  common  focus  at  z=0,  F  is  given 
by50-52 


F(bi,b,Ak,za,Zb)  =  4  b 
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where  -za  and  Zb  are  the  positions  of  the  entrance  and  exit  windows  of  the  mixing  cell  and  b=b]/b2 
is  the  ratio  of  the  input  confocal  parameters.  The  phase  mismatch  between  the  k  vectors  for  the 
driving  waves  (kj  and  k2)  and  the  generated  wave  (k4)  is  Ak  =  k4  -  k'  where  k'=2ki±k2.  The 
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upper  (lower)  sign  refers  to  sum  (difference)  frequency  mixing  here  and  in  Eqs.  (3)  and  (5)  below. 
The  dimensionless  integration  variables  are  given  by  u  =  2z/bj  and  p  =  Vk'/b]  r,  where  z  and  r  are 


the  coordinates  parallel  and  perpendicular  to  the  laser  propagation  direction,  respectively.  When 
absorption  and  nonlinear  refractive  index  effects28  are  negligible  and  Ak  «  k*  the  factors  f(u)  and 
a(u)  are  given  by51'53 


where 


(l+iu)(l±ibu) 

g(u) 
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a(u)  =  f(u)  g(u)  (1+iu) 


(3) 

(4) 
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(5) 


The  integral  F  accounts  for  the  interaction  between  the  driving  polarization  produced  by  the 
pump  beams  and  the  generated  wave,  integrated  over  the  cell  length.  The  factors  a(u)  and  f(u) 
account  for  the  effects  of  focusing  on  the  frequency  conversion,  and  the  phase  mismatch  Ak 
accounts  for  dispersion. 


A .  TUNING  PROFILES 


Simple  expressions  for  F  may  be  obtained  for  two  special  cases,  (I)  sum  frequency 
generation  for  b  =  1  and  (II)  difference  frequency  generation  for  k2  -»  0.  In  both  cases,  Eq.(3) 
simplifies  to  f(u)  =  1  +  i2zt/b  i .  In  addition,  the  generated  beam  is  Gaussian  with  confocal 
parameter  bi.3®  Because  f  is  independent  of  u,  the  integrals  over  u  and  r  can  be  performed 
separately.  In  the  tight  focusing  limit  bi,  b2  «  za,  Zb,  the  phasematching  integral  simplifies  to: 


(I)  Sum  frequency  generation,  b  =  1  (bi  =  b2  =  b) 

(n2  (bAk)2  exp(bAk) 

Fi 


for  Ak<0 
for  Ak>0 


(6) 


(II)  Difference  frequency  generation,  k2  — >  0 
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Eq.  (7)  indicates  that  for  difference  frequency  generation,  F  ~  exp(bjAk)  for  Ak  <  0  and 
F  ~  exp(-b2Ak)  for  Ak  >  0.  For  b=l,  Eq.  (7)  simplifies  further  to  F  =  7t2exp(-blAkl). 


The  source  of  the  asymmetry  in  Eq.  (7)  for  b  ^  1  can  be  understood  by  examining  the 
phase  and  amplitude  variations  in  the  phasematching  integral.  These  phase  and  amplitude 
variations  can  be  visualized  with  the  vectorial  approach  used  by  Ward  and  New,54  which  they  refer 
to  as  vibration  curves.  The  contribution  of  each  differential  length  increment,  du,  for  the  integrand 
over  u  in  Eq.  (2)  is  plotted  as  an  infinitesimal  vector  in  the  complex  plane.  By  plotting  each  vector 
head  to  tail,  the  magnitude  of  the  integral  over  u  is  given  by  the  distance  from  the  tail  of  the  first 
vector  to  the  head  of  the  last  vector.  For  the  limiting  cases  of  Eqs.  (6)  and  (7),  F  is  proportional  to 
the  square  of  this  distance  because  the  integrals  over  u  and  p  are  separable. 


A  good  orientation  to  the  vibration  curve  is  provided  by  the  example  of  third  harmonic 
generation  as  shown  in  Fig.  1(a).  Three  curves  are  plotted  for  values  bAk  =  0  and  ±2.  To  provide 
a  measure  of  physical  distance,  open  circles  are  plotted  at  intervals  of  Au  =  1  (Az=b/2)  throughout 
Fig.  1 .  Each  set  of  three  curves  in  Fig.  1  crosses  at  the  point  u=0  (z=0).  For  bAk  =  0  the 
vibration  curve  is  a  circle,  starting  at  the  bottom,  and  returning  to  the  same  point;  third  harmonic 
generation  is  zero  for  bAk  =  0.  The  source  of  the  phase  variation  producing  this  circle  is  the  Guoy 
effect,55  which  is  an  additional  phase  shift  of  K  that  occurs  when  an  optical  beam  passes  through  a 
focus.  For  third  harmonic  generation,  the  driving  polarization  has  a  phase  shift  of  3k  (k  for  each 
input  beam),  while  the  generated  radiation  has  a  shift  of  only  K.  The  difference  in  phase  between 
the  driving  polarization  and  the  generated  radiation  produces  a  net  phase  shift  of  2jt,  equal  to  the 
angle  around  the  circle.  In  effect,  destructive  interference  between  the  driving  polarization  and  the 
radiation  generated  at  different  positions  in  the  cell  cancel  the  generated  beam  energy.  When  there 
is  phase  mismatch,  dispersion  provides  an  additional  phase  change.  For  Ak  <  0,  the  phase 
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variation  due  to  dispersion  can  counteract  the  phase  variation  due  to  the  Guoy  effect,  leading  to 
non-zero  third  harmonic  generation.  For  Ak  >  0,  phase  changes  from  dispersion  and  the  Guoy 
effect  have  the  same  sign,  leading  to  a  tighter  circle  (spiral),  and  again  no  generated  radiation.  This 
produces  the  well-known  result  that  third  harmonic  generation  only  occurs  for  Ak  <  0  (see  inset 
graph  of  F  versus  biAk  in  Fig.  1(a),  and  Eq  (6)). 

The  case  of  difference  frequency  generation  ((04  =  20)1-0)2)  and  equal  confocal  parameters 
(b=l)  is  shown  in  Fig.  1(b).  For  clarity,  we  present  the  special  case  k2  — >  0.  For  difference- 
frequency  generation,  the  Guoy  phase  shifts  for  002  and  a>i  have  opposite  signs.  Thus  the  driving 
polarization  and  the  generated  wave  both  have  a  phase  shift  of  it,  and  the  vibration  curve  is  a 
straight  line  for  Ak  — »  0  (for  k2  *  0,  there  will  be  small  curvature  to  this  line).  Positive  or 
negative  values  for  Ak  produce  equal  integrated  magnitudes,  so  the  tuning  profile  is  symmetric  (see 
also  inset  graph  in  Fig.  2(b)). 

When  the  confocal  parameters  are  not  equal  (b  ^  1),  the  Guoy  effect  for  the  driving  waves 
occurs  over  different  distances.  In  this  case,  the  vibration  curve  for  bAk  =  0  has  significant 
curvature  as  shown  in  Fig.  1(c),  again  plotted  for  k2  — >  0.  Now  positive  and  negative  values  of 
Ak  do  not  produce  equal  integrated  values.  A  strong  asymmetry  is  produced  in  the  tuning  profile, 
as  is  seen  in  the  inset  graph.  For  positive  Ak,  the  phase  mismatch  has  the  same  sign  as  the  Guoy 
effect  for  the  CO2  beam,  and  b2  governs  the  phasematching  tuning  width.  For  negative  Ak,  the 
phase  mismatch  and  the  Guoy  effect  for  coj  have  the  same  sign,  and  bi  governs  the  tuning  width. 
The  symmetry  between  bi  and  b2  for  Ak<0  and  Ak>0  in  Eq.  (7)  occurs  because  the  Guoy  effects 
of  one  of  the  03]  waves  and  the  generated  wave  cancel  for  k2  — >  0  (they  both  have  confocal 
parameters  bi).  When  k2  &  0,  there  is  an  additional  asymmetry  that  depends  on  the  ratio  k2/k’  in 
Eq.  (5).  However,  for  reasonable  values  of  k2/k'  and  b,  the  net  asymmetry  is  dominated  by  the 
size  of  b  rather  than  k2/k'. 

In  general,  difference  frequency  generation  using  equal  confocal  parameters  is  preferable. 
When  the  confocal  parameters  are  not  equal,  the  peak  difference  frequency  conversion  efficiency  is 
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reduced  by  a  factor  of  approximately  (l+b)(l  +  l/b)/4,  which  is  always  greater  than  one. 

However,  it  may  be  desirable  to  have  an  asymmetric  tuning  profile  in  certain  cases,  for  example  to 
enhance  or  reduce  difference  frequency  mixing  in  the  region  where  sum  frequency  mixing  occurs 
near  bAk=-2. 

B .  PHASE  MISMATCH 

We  have  calculated  the  phase  mismatch  for  Kr,  Ar,  and  H2  using  published  index  of 
refraction  data.  We  include  Ar  because  we  have  used  a  mixture  of  Kr  and  Ar  to  phase  match  VUV 
generation  in  the  region  of  Lyman-a.  For  Ar  we  use  the  Sellmeier  expression  given  by  Leonard56 
which  is  based  on  measurements  performed  down  to  1 10  nm  by  Chashchina  et  al.57  This 
expression  gives  much  better  agreement  with  other  refractive  index  measurements  in  the  range  of 
1 10-120  nm58'60  than  the  Sellmeier  expression  for  Ar  of  Bideau-Mehu  et  al.,61  which  is  based  on 
measurements  down  to  140  nm.  For  H2,  we  use  the  Sellmeier  expression  given  by  Chashchina 
and  Shreider,62  based  on  measurements  down  to  1 15  nm.  This  expression  agrees  with  the  more 
recent  measurements  between  1 10  and  1 15  nm  of  Gordon  et  al.60  to  within  experimental  error. 

For  Kr  we  use  the  analytic  expression  of  Bulanin  and  Kislyakov,63  which  is  based  on  measured 
oscillator  strengths  for  Kr  with  adjustment  of  the  oscillator  strengths  of  the  first  two  resonance 
transitions  to  optimally  match  measured  refractive  index  and  static  polarizability  data.  They  state 
the  expression  should  be  valid  up  to  the  second  resonance  transition.  Their  refractive  index 
expression  agrees  with  the  earlier  calculations  of  Mahon  et  al.64  for  (n-1)  to  within  2-3%  near 
Lyman-a,  which  lies  between  the  first  and  second  resonance  transitions.  Experimental 
measurements  of  the  refractive  index  for  Kr  near  Lyman-a  have  a  very  large  spread58*65’66  and 
don't  provide  a  good  test  of  the  accuracy  of  these  two  calculations  in  this  region. 

In  Fig.  2,  we  plot  the  phase  mismatch  per  atom  or  molecule  C=Ak/N  for  two-photon- 
resonant  four-wave  mixing  with  the  ArF  excimer  laser  in  Kr,  Ar,  and  H2,  where  N  is  the  number 
density  of  each  gas.  Note  that  the  phase  mismatch  for  H2  is  typically  much  less  than  the  phase 
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mismatch  for  Kr.  The  phase  mismatch  for  Kr  becomes  negative  at  123.58  nm,  and  increases  quite 
rapidly  until  the  next  resonance  at  1 16.49  nm,  passing  through  zero  near  120  nm. 

C.  SATURATION  EFFECTS 

The  description  of  conversion  efficiency  given  above  does  not  account  for  absorption  of  the 
VUV  wave  or  saturation  effects28,44,47,49,67"70  that  limit  the  conversion  at  high  powers.  Factors 
contributing  to  saturation  include  the  Kerr  effect,  and  for  resonantly-enhanced  processes,  power 
broadening.  Stark  shifts,  two-photon  absorption,  multiphoton  ionization,  stimulated  hyper-Raman 
scattering,  amplified  spontaneous  emission,  parametric  four-wave  mixing,  self-defocusing,  and 
related  nonlinear  effects.  These  processes  lead  to  changes  in  refractive  index,  population 
redistributions,  and  depletion  of  the  pump  waves  that  vary  both  along  and  transverse  to  the  pump 
beams.  Thus,  deviations  from  the  simple  descriptions  for  conversion  efficiency  are  to  be  expected. 


3.  VUV  GENERATION 

Our  experimental  arrangement  for  VUV  generation  is  shown  in  Fig.  3.  The  ArF  excimer 
laser  (Lambda  Physik  EMG-150  MSC)  has  separate  oscillator  and  amplifier  discharge  chambers. 
To  improve  the  beam  quality  for  nonlinear  frequency  conversion,  the  laser  has  been  modified  to 
include  a  spatial  filter  before  the  amplifier  and  to  triple  pass  the  amplifier.3  The  laser  beam 
divergence  is  about  5  times  over  the  diffraction  limit.  Energies  of  up  to  100  mJ  are  obtained, 
although  more  typical  operational  energies  are  about  20  mJ.  The  second  laser  is  a  Nd:YAG- 
pumped  dye  laser  (Quanta-Ray  DCR II  and  PDL).  To  produce  longer  VUV  wavelengths,  the 
second  harmonic  of  the  dye  laser  is  used.  A  control  circuit  (Lambda  Physik  model  EMG  97)  is 
used  to  minimize  drift  of  the  excimer  laser  pulse  timing.  The  beam  path  of  the  ArF  laser  is  purged 
with  argon  to  minimize  the  effects  of  oxygen  Schumann-Runge  absorption.3  The  ArF  and  dye 
beams  are  combined  using  a  high-power  45°  ArF  dielectric  mirror  and  focused  with  a  single  50-cm 
or  100-cm  lens  into  the  mixing  cell.  A  telescope  is  used  to  adjust  the  dye  laser  spot  size  and 
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compensate  for  the  chromatic  aberration  of  the  focusing  lens.  A  half  waveplate  is  used  to  produce 
horizontal  polarization  for  the  dye  laser  to  match  the  horizontal  polarization  of  the  ArF  laser.  When 
the  dye  laser  output  is  frequency  doubled,  the  waveplate  is  not  used  since  the  frequency  doubling 
process  provides  a  horizontally  polarized  output  beam. 

To  monitor  whether  the  ArF  wavelength  is  resonant  with  the  desired  two-photon  transition, 
ions  produced  through  2+1  resonantly-enhanced  multi-photon  ionization  (REMPI)  are  monitored 
using  an  electrode  biased  at  +1 .5  to  +3  V  relative  to  the  cell.  The  electrode  is  placed  before  the 
laser  focus  to  avoid  spectral  resolution  loss  due  to  power  broadening  at  the  focus. 

At  the  output  of  the  mixing  cell,  a  single  off-center  MgF2  lens  is  used  to  either  refocus  or 
collimate  the  VUV  beam,  to  dispersively  separate  the  VUV  and  pump  beams,  and  as  a  window  to 
the  cell.  This  single-optic  arrangement  minimizes  the  amount  of  absorption  loss  for  the  VUV 
beam.  This  lens  is  mounted  on  a  vacuum  x-y-z  adjustment  stage  to  allow  adjustment  of  the  VUV 
beam  position. 

To  prevent  a  potentially  rapid  decrease  in  the  VUV  transmission  caused  by  trace  amounts  of 
oil,  liquid  nitrogen  and  an  acetone/dry  ice  bath  are  used  to  isolate  the  VUV  beamline  and  mixing 
cell,  respectively,  from  the  vacuum  pumps.  Acetone/dry  ice  is  used  on  the  mixing  cell  to  avoid  a 
significant  boiloff  of  liquid  nitrogen  when  the  gas  in  the  cell  is  replaced.  High  purity  (99.997  or 
better)  Kr,  H2  and  Ar  were  used  in  the  mixing  cell. 

During  previously  reported  work,4  we  found  decay  of  the  MgF2  optics  transmission  after 
relatively  short  VUV  irradiation  (approximately  15  min).  We  built  a  VUV  generation  cell  with  no 
exit  window  cell  to  study  this  effect.  This  cell  consisted  of  a  30-cm  long,  4-mm  inside  diameter 
pyrex  tube,  sealed  with  a  window  at  the  entrance  and  open  at  the  exit.  The  cell  was  filled  at  the 
window  end  using  a  pulsed  valve.  2+1  REMPI  (resonantly  enhanced  multiphoton  ionization)  in 
H2  was  used  to  monitor  the  cell  pressure.  This  cell  functioned  well  for  VUV  generation,  and  was 
used  to  study  the  transmission  of  optics  inserted  into  the  generated  VUV  beam  using  a  vacuum 
feedthrough  fitting.  When  we  tested  recently  purchased  VUV  grade  MgF2  optics  with  the 
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windowless  cell,  we  did  not  find  any  rapid  transmission  decay.  Nor  have  we  found  rapid  optics 
decay  problems  using  newer  VUV-grade  optics  in  our  more  recent  VUV  generation  experiments. 
This  may  be  indicative  of  improvements  in  the  quality  of  VUV  grade  optical  material  since  the 
previous  studies. 

We  have,  however,  observed  a  decay  of  the  response  for  the  pyroelectric  energy  probes 
used  to  measure  absolute  VUV  energies.71  On  exposure  to  high  power  Lyman-a  radiation,  the 
probe  sensitivity  drops  for  the  VUV  as  well  as  for  the  pump  beams.  The  sensitivity  drops  over  2 
to  3  minutes  during  VUV  exposure,  and  recovers  when  the  VUV  is  removed,  for  example  when  a 
delay  is  introduced  between  the  timing  of  the  two  pump  laser  pulses.  Detection  of  the  VUV  with  a 
second  detector  (observation  of  scattered  VUV  light  with  a  Lyman-a  bandpass  filter  and  a  KBr 
photomultiplier  tube)  confirmed  that  the  drop  in  pyroelectric  probe  response  was  not  due  to  loss  of 
VUV  power.  The  drop  in  response  was  observed  for  both  a  Laser  Precision  model  RjP735  probe 
in  vacuum  and  for  a  Molectron  J3-09  probe  operated  under  an  Ar  purge.  The  sensitivity  of  the 
Laser  Precision  probe  recovered  over  periods  of  tens  of  minutes  to  hours,  while  the  Molectron 
probe  recovered  in  less  than  ten  minutes.  Under  certain  circumstances,  the  pulsed  electrical  output 
from  the  Molectron  probe  inverted  from  positive  to  negative.  For  the  Laser  Precision  probe,  the 
recovery  time  increased  with  successive  irradiation  episodes,  and  the  signal  did  not  recover  by 
illuminating  a  different  portion  of  the  detector. 

The  nature  of  this  sensitivity  loss  for  these  pyroelectric  probes  is  not  completely 
understood.  Pyroelectric  energy  meters  have  been  used  for  measurements  of  high  power  VUV  at 
wavelengths  of  130  nm  and  above.4’30  These  reports  do  not  mention  decay  in  the  response  of  the 
energy  meters,  indicating  that  the  sensitivity  decline  of  pyroelectric  energy  meters  may  not  occur  at 
130  nm  and  above.  Other  VUV  energy  or  power  detectors  exist  such  as  Schottky  photodiodes 
used  for  synchrotron  radiation72  and  ionization  cells.73  However,  these  detectors  may  saturate  for 
the  kW  power  levels  described  here. 
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We  have  used  both  Kr  and  H2  for  VUV  generation,  but  we  have  not  attempted  mixing  in 
HD.  For  mixing  in  Kr,  two-photon  resonances  in  the  4p5  6p[5/2,2],  4p5  6p[3/2,2],  and 
4p5  6p[  1/2,0]  states  may  be  accessed.  However,  the  4p^  6p[5/2,2]  and  4p-’  6p[  1/2,0]  lie  in  the 
wings  of  the  ArF  tuning  range  where  the  output  power  is  reduced;  therefore  we  have  used  only  the 
4p5  6p[3/2,2]  transition  for  mixing  in  Kr.  Seven  two-photon  resonances  for  H2  lie  within  the 
tuning  range  of  the  ArF  laser  corresponding  to  Q-branch  transitions  to  v=6  and  7  in  the  E-F  stated 
We  use  the  Q(l)  (6,0)  transition  most  commonly  for  VUV  generation  because  of  the  large 
Boltzmann  fraction  in  J=l,  and  because  the  Q(l)  (7,0)  line  lies  in  the  wing  of  the  tuning  curve. 
(Note  that  the  vibrational  level  v=6  in  the  E,F  state  is  sometimes  labeled  vg  =  2,  which 
corresponds  to  the  vibrational  number  of  the  inner  (E)  well  of  the  double  minimum  E,F  state.) 

Both  the  Kr  6p[3/2,2]  and  H2  Q(l)  (6,0)  transitions  overlap  with  Schumann-Runge  bands 
in  oxygen,  and  the  ArF  path  must  be  purged  with  argon.  The  absorption  is  particularly  strong  for 
the  Kr  6p[3/2,2],  which  requires  a  more  aggressive  Ar  purge  then  when  tuning  to  the  H2  Q(l) 

(6,0)  transition. 

The  selection  of  Kr  or  H2  as  the  nonlinear  mixing  medium  depends  on  the  application.  As 
discussed  below,  the  nonlinear  susceptibilities  for  the  Kr  6p[3/2,2]  and  H2  Q(l)  (6,0)  transitions 
are  similar  in  magnitude  near  Lyman-a.  When  mixing  in  H2,  VUV  amplified  spontaneous 
emission  (ASE)  is  generated,  as  is  discussed  in  Section  3A.  In  the  region  where  the  ASE  is  the 
strongest,  approximately  135-160  nm,  mixing  in  Kr  is  generally  preferred  to  avoid  the  additional 
VUV  from  ASE.  Near  135  nm  and  below,  the  phase  mismatch  for  Kr  becomes  quite  large  due  to 
the  resonance  lines  at  123.58  and  1 16.49  nm  (see  Fig.  2).  In  this  region  mixing  in  H2  provides 
higher  powers  than  for  Kr.  Using  a  pyroelectric  energy  meter,  we  have  directly  measured  65  at 
133  nm  when  mixing  in  H2.  Mixing  in  H2  does  not  work  well  in  the  region  of  Lyman-a,  as 
discussed  in  Section  3B.  In  the  regions  below  the  Kr  resonances,  the  phase  mismatch  is  negative 
in  Kr,  which  allows  the  possibility  of  phase  matching  by  adding  a  gas  with  positive  phase 
mismatch  as  described  in  Section  3B. 
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A. 


ASE 


When  the  ArF  is  tuned  to  two-photon  resonances  in  the  E,F  state  for  H244,74  or  HD,75 
sufficient  population  may  be  transferred  from  the  ground  state  to  cause  ASE  in  the  VUV.  This 
ASE  can  occur  due  to  population  transfer  from  the  E,F  state  to  the  B  state  through  near  infrared 
ASE  followed  by  ASE  from  the  B  to  X  state  (Lyman  band  emission).  Although  ASE  from  the  C 
state  to  the  X  state  (Wemer  bands)  also  can  occur  through  electron  collisional  transfer  from  the  E,F 
state  to  the  C  state,44  we  have  only  observed  Lyman  band  emission.  The  Lyman  band  ASE  occurs 
predominantly  at  wavelengths  from  about  135  nm  to  160  nm.  This  ASE  may  be  a  nuisance  when 
performing  difference  frequency  generation.  Although  the  optimum  pressure  for  generation  of 
ASE  tends  to  be  higher  than  the  optimum  pressure  for  difference  frequency  generation  4  it  is  not 
possible  to  perform  difference  frequency  generation  in  H2  without  significant  ASE  generation. 
However,  the  VUV  ASE  can  also  serve  a  beneficial  role  for  checking  VUV  system  performance. 
The  VUV  ASE  may  be  used  to  test  the  transmission  of  optics  and  mixing  gases  in  the  VUV 
independently  from  difference  frequency  mixing,  which  requires  careful  temporal  and  spatial 
overlaps  of  ArF  and  dye  beams.  The  ASE  may  also  be  used  as  a  source  of  discrete  VUV  lines. 

A  spectrum  of  ASE  lines  taken  with  a  VUV  spectrometer  (SPEX  model  1500SP)  is  shown 
in  Fig.  4.  This  measurement  was  taken  with  a  100  cm  lens  and  a  H2  pressure  of  approximately 
200  mbar  when  the  excimer  laser  was  tuned  to  the  E,F  (v'=6,  J'=l)  level.  Also  shown  are 
expected  line  positions  for  selected  B  — >  X  (Lyman  band)  transitions  from  Dabrowski.76  The 
spectrum  in  Fig.  4  has  not  been  corrected  for  the  spectral  response  of  the  spectrometer  and 
photomultiplier  tube.  Although  the  E,F  — >  B  ASE  can  populate  either  J=2  or  J=0  in  the  B  state, 
we  have  only  observed  VUV  ASE  from  J=2  in  the  form  of  P(3)  and  R(l)  emission;  P(l)  lines  are 
not  observed.  B  — >  X  transitions  from  v’=0  to  v"=4-8  and  v'=l  to  v"=5-9  are  observed.  For  each 
pair  of  lines,  the  P(3)  emission  is  stronger  than  the  R(l)  emission.  Only  the  P(3)  line  is  observed 
for  the  (1,5)  and  (1,9)  bands.  A  peak  near  162.8  nm  appears  to  be  the  (2,10)  P(3)  line.  At  this 
spectral  resolution,  it  is  not  possible  to  resolve  the  (2,9)  P(3),  (2,8)  P(3),  and  (2,8)  R(l)  lines 
from  the  (0,8)  R(l),  (0,7)  P(3),  and  (0,7)  R(l)  lines,  respectively.  The  unassigned  line  at  153.8 


12 


nm  is  due  to  a  stray  reflection  of  the  ArF  laser  in  the  spectrometer;  we  have  not  identified  the  line 
near  161.4  nm. 

In  addition  to  the  lines  shown  in  Fig.  4,  we  have  also  observed  ASE  emission  on  the  (1,2) 
Lyman  band  as  shown  in  Fig.  5.  ASE  on  this  band  has  not  been  reported  previously,  to  our 
knowledge.  The  general  pattern  of  observed  Lyman  band  ASE,  such  as  reduced  intensity  between 
the  (1,6)  and  (1,2)  lines  and  the  limited  range  for  emission  from  v'=2  agrees  with  the  variation  in 
the  Franck-Condon  factors  (overlap  of  the  vibrational  wavefunctions  in  the  B  and  X  states)  for 
these  bands  given  by  Spindler.77 

B .  LYMAN-oc  GENERATION 

Lyman-a,  corresponding  to  the  shortest  wavelength  connecting  the  ground  state  of  H 
atoms  at  121.567  nm  is  an  important  wavelength  for  spectroscopy,  photochemistry,  and  plasma 
diagnostics.  Nonlinear  production  of  Lyman-a  has  been  pursued  by  a  number  of  groups  using 
frequency  tripling;^’^®’^, 64,65, 78-87  Raman  shifting;9'77  electromagnetically  induced 
transparency;88  and  resonantly-enhanced  four-wave  mixing  in  Mg,89  Hg,7’90’91  and 
jq..  13,14,29,71  \y e  have  investigated  the  performance  of  our  VUV  source  at  this  wavelength. 

VUV  generation  with  H2  produces  about  two  orders  of  magnitude  more  energy  than  with 
Kr  in  the  region  of  Lyman-a.  The  VUV  energy  generated  near  Lyman-a  in  65  mbar  of  pure  H2  is 
shown  as  a  function  of  VUV  wavelength  as  the  solid  line  in  Fig.  6.  A  dramatic  dip  in  VUV 
generations  occurs  at  Lyman-a  due  to  the  production  of  H  atoms.  There  is  a  strong  asymmetry  to 
this  dip,  with  reduced  efficiency  in  the  long  wavelength  side.  Others  also  have  reported 
asymmetric  dips  near  Lyman-a  for  VUV  generation  in  H2.9,97,9^  The  asymmetry  is  due  to 
contributions  of  H  atoms  to  the  phase  mismatch.  The  index  of  refraction  of  the  H  atoms  varies 
rapidly  near  the  strongly-allowed  Lyman-a  transition,  which  produces  a  profound  effect  on  the 
phase  matching.  On  the  long  wavelength  side  of  Lyman-a,  the  phase-mismatch  increases  and 
conversion  efficiency  drops.  On  the  short  wavelength  side,  H  atoms  make  a  negative  contribution 


13 


to  Ak  and  improve  the  conversion  efficiency  through  reduced  phase  mismatch.  This  interpretation 
is  supported  by  measurements  in  negatively  dispersive  gases:  when  generating  VUV  in  Kr  with  H2 
added,  the  asymmetry  is  reversed  (see  Fig.  7)  because  Ak  for  Kr  is  negative  in  this  spectral  region. 

We  have  modeled  the  shape  of  the  tuning  profile  in  Fig.  6  by  calculating  values  for  the 
phase  matching  function,  F,  as  a  function  of  VUV  wavelength.  This  calculation  uses  the  phase 
mismatch  per  molecule  for  H2  from  Fig.  2,  the  measured  H2  pressure,  and  the  confocal  parameters 
determined  for  the  fit  in  Fig.  9.  The  refractive  index  for  H  atoms  is  calculated  from  Equation  (1)  of 
Leonard56  using  calculated  oscillator  strengths,94  and  neglecting  contributions  from  the 
continuum.  This  calculation  is  sufficiently  accurate  for  our  purpose.  The  H  atom  density  is  used 
as  a  fitting  parameter.  The  resulting  tuning  profile  for  a  H  atom  number  density  equal  to  0.004 
times  the  H2  number  density  is  shown  as  the  dashed  line  in  Fig.  6.  The  fit  agrees  qualitatively 
with  the  measurement  except  very  near  Lyman-a.  Some  disagreement  between  the  measurement 
and  fit  is  not  surprising  since  the  fit  includes  many  approximations.  For  instance,  the  fit  assumes  a 
uniform  H  atom  density  over  the  full  mixing  length,  while  the  H  atoms  are  likely  to  drop 
significantly  away  from  the  focus. 

There  is  additional  reproducible  structure  in  Fig.  6  that  we  believe  is  due  to  H2  at  121.0, 
121.2,  122.5,  and  123  nm.  We  have  not  been  able  to  conclusively  assign  these  features.  Possible 
sources  of  this  structure  are  competition  from  2  +  1  (2  ArF  photons  and  one  dye  photon) 
absorption  into  H2  Rydberg  states  or  Werner  band  absorption  by  excited  states  produced  by  the 
VUV  ASE  described  in  Section  3A.  Expected  line  positions  from  the  literature  for  the  Rydberg95 
and  Werner76  transitions  are  shown  in  Fig.  6. 

We  have  tried  unsuccessfully  to  avoid  production  of  H  atoms  near  Lyman-a.  Reduction  in 
H2  pressure  and  ArF  laser  power  and  slight  detunings  of  the  ArF  laser  from  resonance  met  with 
limited  success,  but  none  eliminated  the  drop  in  power  at  Lyman-a.  Adjusting  synchronization  of 
the  ArF  and  dye  laser  to  encourage  mixing  before  significant  H  atom  production  had  no  apparent 
effect.  Although  we  could  produce  5-15  pJ  near  Lyman-a,  the  maximum  energy  we  were  able  to 
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produce  at  Lyman-a  in  H2  was  0.2  |iJ  (see  Fig.  8).  The  zero  level  for  the  measurements  in  Fig.  8 
was  determined  by  changing  the  delay  between  the  two  lasers  to  prevent  nonlinear  mixing. 

We  have  not  identified  the  mechanism  for  generation  of  the  H  atoms.  Emission  at  Balmer- 
a  (656  nm)  was  observed  in  the  focal  region.  The  intensity  of  the  Balmer-a  radiation  increases  as 
roughly  the  third  power  of  density  before  saturating  and  slowly  decreasing  above  13  mbar.  Datkos 
et  al.70  have  observed  Lyman-a  from  ArF-irradiated  H2,  which  they  attribute  to  electron 
attachment  to  high  Rydberg  states  produced  following  absorption  of  three  ArF  photons. 

Note  that  from  Fig.  8  we  observe  VUV  generation  even  directly  at  Lyman-a,  while  the  H 
atoms  concentrations  estimated  from  Fig.  6  should  absorb  all  Lyman-a.  One  possible  explanation 
is  that  H  atoms  are  predominantly  produced  at  the  laser  focus  while  some  Lyman-a  is  produced 
after  the  focal  region.  To  summarize  our  results  in  pure  H2,  we  produce  5-15  pJ  near  Lyman-a, 
with  a  1-2  order  of  magnitude  dip  on  resonance  where  peak  energies  of  about  0.2  (Jj  are  produced. 

The  lower  efficiency  for  Kr  relative  to  H2  near  Lyman-a  is  primarily  due  to  the  phase 
mismatch,  not  the  nonlinear  susceptibility  At  121.46  nm,  H2  has  a  maximum  VUV  power  at 
a  pressure  of  about  65  mbar,  while  Kr  has  a  maximum  at  about  7  mbar  in  this  region  (both  using  a 
50  cm  lens).  Thus  without  phase  matching  Kr  has  an  experimentally  measured  optimum  pressure 
that  is  about  10  times  lower  than  for  H2.  This  is  twice  the  calculated  difference  in  phase  mismatch 
between  Kr  and  H2  given  in  Fig.  2  of  Ckt  ~  -5  CH2.  The  discrepancy  is  primarily  due  to  the  phase 
mismatch  from  H  atoms  at  the  wavelength  of  121.46  nm  (see  Fig.  6).  Because  the  VUV  power 
scales  with  the  square  of  the  gas  density,  phase  mismatch  alone  should  cause  H2  to  be  a  factor  of 
about  100  times  more  efficient  than  Kr.  Because  we  find  that  H2  produces  about  100  time  more 
energy  than  Kr,  we  expect  the  %0)  for  Kr  and  H2  should  be  roughly  equivalent  near  Lyman-a. 

We  generated  high  energies  at  Lyman-a  using  a  phase  matched  mixture  of  Kr  and  Ar.  A 
mixture  of  Kr  and  H2  cannot  be  used  because  H  atoms  are  generated  even  when  the  ArF  laser  is 
tuned  to  the  Kr  resonance  (see  Fig.  7).  Kr  and  Ar  were  added  gradually  to  the  cell  to  enhance 
mixing,  using  typically  four  additions  of  each  gas  in  alternation.  A  mixture  ratio  of  3.9: 1  Ar/Kr 
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phase  matches  at  Lyman-a,  resulting  in  about  two  orders  of  magnitude  increase  in  energy  over 
pure  Kr.  By  comparison,  from  the  phase  mismatch  from  Fig.  2,  we  expect  a  ratio  of  4.7:1  Ar/Kr. 
The  generated  energy  at  a  Kr  partial  pressure  of  130  mbar  is  shown  in  Fig.  9  as  a  function  of  VUV 
wavelength.  Seven  |iJ  are  obtained  at  Lyman-a.  Prospects  for  further  increase  in  energies  and 
implications  for  the  tuning  range  are  described  elsewhere.71 

There  is  a  pronounced  asymmetry  in  the  tuning  profile  of  Fig.  9.  This  asymmetry  is  too 
large  to  be  produced  by  curvature  in  the  dispersion  of  the  Kr/Ar  mixture.  Such  an  asymmetry  can 
be  produced  by  unequal  confocal  parameters  for  the  input  beams  as  described  in  Section  2A.  We 
fit  a  phase  matched  tuning  curve  to  the  data  of  Fig.  9  using  Eqs.  (2)  to  (5)  and  the  phase  mismatch 
data  of  Fig.  2.  This  fit  is  shown  as  the  dashed  line  in  Fig.  9.  The  fit  uses  values  of  b]  =  0.6  cm 
and  b  =  4.  We  have  also  modeled  the  effects  of  an  offset  between  the  foci  of  the  input  beams51  to 
study  possible  asymmetries.  Although  an  offset  reduces  the  frequency  conversion  efficiency,  it 
produces  no  significant  asymmetry. 

4 .  VUV  MULTIPHOTON  SPECTROSCOPY 

We  have  applied  VUV  radiation  generated  with  two-photon-resonant  difference  frequency 
mixing  to  two  types  of  multiphoton  excitation  deep  in  the  VUV;  1+1  REMPI  (resonantly-enhanced 
multi-photon  ionization)  in  Xe  at  147  nm  and  two-photon  excited  fluorescence  in  Ne  at  133  nm. 

For  both  of  these  experiments,  separation  of  the  VUV  and  pump  beams  was  performed  with  a 
MgF2  Pellin  Broca  prism  rather  than  the  off-axis  lens  shown  in  Fig.  3. 

A.  1+1  REMPI  IN  XENON 

We  used  two  photons  at  147  nm  to  perform  1+1  REMPI  through  the  5p56s[3/2,l]  state  of 
xenon  at  68,045  cm*1,  as  shown  in  the  inset  in  Fig.  10.  The  147-nm  radiation  was  produced 
through  mixing  in  Kr.  Because  VUV  radiation  can  readily  ionize  many  molecular  species, 
background  ion  signals  are  a  major  consideration  for  VUV  REMPI.  With  energies  of  only  3  jjJ, 
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we  have  obtained  signal-to-noise  ratios  of  over  50  for  the  ion  signal,  a  good  indication  that  this 
approach  can  give  useful  signals.  A  1+1  REMPI  spectrum  and  an  absorption  spectrum  for  Xe  at  a 
cell  pressure  of  0.065  mbar  are  shown  in  Fig.  10.  The  ion  signal  is  collected  with  a  single 
electrode  biased  at  25  V/cm.  The  strong  resonance  absorption  in  xenon  leads  to  a  dip  in  the 
REMPI  signal  at  line  center.  At  higher  pressures,  the  ion  signal  drops  to  zero  at  line  center.  The 
power  dependence  of  the  ion  signal  from  xenon  at  a  pressure  of  0.13  mbar  is  shown  in  Fig.  11. 
The  signal  follows  the  expected  square  dependence  on  the  VUV  power. 

B .  TWO-PHOTON-EXCITED  FLUORESCENCE  IN  NEON 

We  have  also  performed  two-photon  excited  fluorescence  from  neon.  Neon  is  the  second 
most  difficult  element  to  excite  from  the  ground  state,  after  helium.  Mixing  in  H2  was  used  to 
generate  approximately  20  jllJ  at  133  nm.  Two  photons  excite  the  2p-'3p[3/2,2J  state  which  lies 
150,858.5  cm"1  above  the  ground  state.  Fluorescence  to  the  2p53s[3/2,l]  state  at  610  nm  is 
detected  with  a  photomultiplier  and  a  700  nm  short  pass  filter.  The  Ne  pressure  was  approximately 
130  mbar.  Tuning  the  VUV  demonstrates  a  linewidth  at  133  nm  of  1.8  cm’1.  This  linewidth  is 
predominantly  due  to  the  linewidths  of  the  ArF  (-1  cm'1)  and  frequency-doubled  dye  (~1.5  cm'1) 
lasers  .  Additional  details  of  these  measurements  are  presented  elsewhere.4 


5.  SUMMARY  AND  CONCLUSIONS 

We  have  shown  that  two-photon-resonant  difference  frequency  mixing  using  an  ArF  laser 
is  a  practical  method  for  production  of  high  energies  in  the  VUV.  Mixing  in  Kr  and  H2  each  have 
their  relative  advantages.  H2  has  a  smaller  phase  mismatch,  but  the  production  of  strong  VUV 
ASE  in  the  wavelength  range  of  -135-160  nm  appears  unavoidable  in  H2.  The  negative  phase 
mismatch  in  Kr  below  the  first  resonance  at  123.58  nm  allows  phase  matching  using  gas  mixtures. 
We  have  not  explored  mixing  in  HD. 
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We  have  directly  measured  over  65  |llJ  at  133  nm  when  mixing  in  pure  H2,  and  7  jjJ  at 
Lyman-a  (121.567  nm)  when  mixing  in  a  phase  matched  mixture  of  Kr  and  Ar.  Mixing  in  H2  at 
Lyman-a  is  hampered  by  the  production  of  H  atoms,  leading  to  poor  VUV  generation  and  an 
asymmetric  tuning  profile  due  to  the  phase  mismatch  of  the  H  atoms.  A  physical  basis  is  presented 
for  the  asymmetric  tuning  profile  for  the  phase  matched  mixture  of  Kr  and  Ar.  The  utility  of  this 
source  is  demonstrated  in  two  multiphoton  VUV  experiments:  1+1  REMPI  in  Xe  at  147  nm  and 
two-photon-excited  fluorescence  in  neon  at  133  nm. 
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FIGURE  CAPTIONS 


Figure  1.  Vibration  curves  for  three  cases  of  four-wave  mixing,  (a)  third  harmonic  generation,  (b) 
difference  frequency  mixing  for  b=l,  and  (c)  difference  frequency  mixing  for  b=4.  Inset 
graphs  show  the  phasematching  integral  F  as  a  function  of  bi  Ak. 

Figure  2.  Phase  mismatch  per  atom  or  molecule  for  Kr,  Ar,  and  H2  for  two-photon-resonant 

difference  frequency  mixing  with  an  ArF  excimer  laser.  Regions  of  negative  phase  mismatch 
for  Kr  are  indicated  with  dashed  lines. 

Figure  3.  Experimental  arrangement  for  vuv  generation.  Wavelengths  shown  are  for  Lyman-cx 
generation. 

Figure  4.  Vacuum  ultraviolet  amplified  spontaneous  emission  generated  on  pumping  the  E,F  1Ig+ 
(v  =  6)  <—  X  1Zg+  (v"=0)  Q(l)  transition  in  H2.  The  emission  strengths  have  not  been 
corrected  for  the  spectral  response  of  the  spectrometer  and  photomultiplier  tube.  Lyman  band 
line  positions  are  shown  at  the  top  of  the  figure. 

Figure  5.  Vacuum  ultraviolet  amplified  spontaneous  emission  on  the  (1,2)  Lyman  band. 

Figure  6.  VUV  energy  generated  in  pure  H2  in  the  region  of  Lyman-a  at  a  pressure  of  65  mbar. 
The  dashed  line  is  a  theoretical  tuning  profile  using  Eqs.  (2)-(5).  Possible  line  assignments 
for  spectral  features  are  shown  at  the  top  of  the  figure. 

Figure  7.  VUV  generation  near  Lyman-a  using  pure  Kr  and  a  mixture  of  Kr  and  H2. 

Figure  8.  Expanded  view  of  Lyman-a  generation  in  pure  H2.  The  upper  curve  has  been 
multiplied  by  a  factor  of  ten. 
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Figure  9.  VUV  generation  in  a  phasematched  mixture  of  Kr  and  Ar  at  Lyman-a  (solid  line)  and  fit 
tuning  profile  from  Eqs.  (2)-(5). 

Figure  10.  Transmission  (a)  and  1+1  REMPI  (b)  spectra  for  xenon  with  147-nm  radiation. 

Energy  level  diagram  for  1+1  REMPI  is  shown  as  inset. 

Figure  11.  Power  dependence  of  ion  signal  from  1+1  REMPI  in  xenon. 
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Figure  4.  Vacuum  ultraviolet  amplified  spontaneous  emission  generated  on  pumping  the  E,F  ^g*  (v'=  6) 
<—  X  JEg-1'  (v"=0)  Q(1)  transition  in  H2.  The  emission  strengths  have  not  been  corrected  for  the  spectral 
response  of  the  spectrometer  and  photomultiplier  tube.  Lyman  band  line  positions  are  shown  at  the  top 

of  the  figure. 
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Figure  6.  VUV  energy  generated  in  pure  H2  in  the  region  of  Lyman-a  at  a  pressure  of  65  mbar.  The 
dashed  line  is  a  theoretical  tuning  profile  using  Eqs.  (2)-(5).  Possible  line  assignments  for  spectral 

features  are  shown  at  the  top  of  the  figure. 
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Figure  8.  Expanded  view  of  Lyman-a  generation  in  pure  H2.  Upper  curve  has  been  multiplied  by  a  factor 

of  ten. 
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Figure  11.  Power  dependence  of  ion  signal  from  1+1  REMPI  in  xenon. 
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We  report  high-power  vacuum-ultraviolet  (vuv)  generation  at  the  Lyman- a  wavelength  of  121.6  nm,  using 
a  simple  experimental  system,  vuv  radiation  is  produced  through  two-photon- resonant  difference-frequency 
mixing  with  a  tunable  ArF  excimer  laser  and  a  Nd-YAG-pumped  dye  laser.  Using  phase-matched  mixtures 
of  Kr  and  Ar  at  a  total  pressure  of  650  mbar,  we  produced  7-/laJ  energies  at  Lyman-a  in  approximately  5  ns 
(1.3  kW),  as  measured  directly  with  a  pyroelectric  energy  probe.  Measurements  indicate  that  higher  powers 
are  possible  with  system  optimization.  A  tuning  range  of  0.1  nm  was  achieved  for  a  fixed  gas  mole  fraction  at  a 
total  pressure  of  650  mbar.  Qualitative  agreement  is  found  between  measured  tuning  profiles  and  theoretical 
predictions.  ©  1998  Optical  Society  of  America 
OCIS  codes:  190.0190,  140.7240,  140.2180. 


Development  of  coherent  vacuum-ultraviolet  (vuv) 
sources  continues  to  provide  new  opportunities  for 
vuv  spectroscopy  and  photochemistry.  Higher  powers 
and  broad  tunability  are  extending  techniques  such  as 
multiphoton  detection,  planar  imaging,  and  photolysis 
deeper  into  the  vuv.  Hydrogen  Lyman-a  radiation 
(121.567  nm)  is  especially  interesting  because  it  per¬ 
mits  single-photon  detection  of  ground-state  H  atoms, 
which  are  important  in  arc-jet  plumes,  tokamaks, 
combustion,  molecular  dynamics,  and  astrophysics. 
Several  researchers  have  generated  Lyman- a  with  dye 
lasers  and  phase-matched  mixtures  of  Kr  with  Xe  or  Ar 
by  use  of  frequency  tripling1"3  or  difference-frequency 
mixing.4 

Alternatively,  substantial  vuv  power  can  be  gener¬ 
ated  with  an  ArF  excimer  laser  and  a  dye  laser  by  two- 
photon-resonant  difference-frequency  mixing,  = 
2^ArF  _  ^dye  -5"7  One  can  enhance  this  process  by  tun¬ 
ing  the  ArF  laser  to  a  two-photon  resonance  in  H2, 
HD,  or  Kr,  and  tuning  is  provided  by  variation  of  the 
wavelength  of  the  dye  laser.  By  varying  the  dye-laser 
wavelength,  one  can  produce  vuv  wavelengths  of  110 
to  180  nm.  Because  of  the  importance  of  the  Lyman- 
a  wavelength,  our  objective  was  to  evaluate  the  en¬ 
ergies  available  at  Lyman- a  with  this  technique.  Kr 
has  negative  dispersion  near  Lyman-a,  allowing  phase 
matching  by  use  of  gas  mixtures. 

A  schematic  of  the  experimental  arrangement  is 
shown  in  Fig.  1.  Because  of  the  importance  of  beam 
quality  for  nonlinear  frequency  conversion,  we  modi¬ 
fied  the  ArF  laser  (Lambda  Physik  EMG-150  MSC)  to 
provide  spatial  filtering  of  the  oscillator  output  beam 
and  triple  pass  the  amplifier.8  This  laser  is  tuned  to 
two-photon  resonances  in  the  mixing  gas.  The  third 
photon  is  provided  by  a  Nd:YAG-pumped  dye  laser 
(Quanta-Ray  DCR  II  and  PDL)  using  Lambda  Physik 
LC4700  dye.  Output  energies  of  the  ArF  laser  and  the 
dye  laser  are  typically  20  and  10  mJ,  respectively,  and 
linewidths  are  of  the  order  of  1  cm-1.  The  beams  are 
combined  on  a  high-power  ArF  dielectric  mirror  and 
focused  with  a  50-cm  lens  into  the  gas  cell.  We  ver¬ 
ify  spatial  overlap  by  temporarily  directing  portions  of 


the  beams  through  a  pinhole  by  use  of  a  wedge.  Tem¬ 
poral  overlap  is  maintained  by  a  drift-control  circuit 
and  verified  by  a  fast  photodiode.  To  minimize  op¬ 
tics  transmission  losses,  we  use  a  single  MgF2  plano¬ 
convex  lens  (50-cm  focal  length  at  250  nm)  to  provide 
an  exit  window  for  the  mixing  cell,  collimate  the  vuv 
beam,  and  dispersively  separate  the  vuv  beam  from  the 
ArF  and  the  dye  beams.  An  aperture  of  7 -mm  diam¬ 
eter  spatially  selects  the  vuv  beam,  which  is  detected 
by  a  Laser  Precision  RjP735  pyroelectric  energy  meter 
and  a  KBr  photomultiplier  tube.  Two  new  vuv-grade 
MgF2  lenses  purchased  from  Optovac  and  Harshaw- 
Bicron  were  used  for  these  measurements,  and  both 
initially  yielded  >70%  transmission  at  Lyman-a.  To 
date  we  have  used  only  the  Optovac  lens  and  have  not 
encountered  the  same  degradation  that  was  encoun¬ 
tered  in  the  research  reported  in  Ref.  9,  perhaps  owing 
to  improvements  in  MgF2  growth  techniques. 

Although  the  pyroelectric  energy  probe  allows  rela¬ 
tively  straightforward  measurement  of  the  vuv  ener¬ 
gies  ,  it  is  not  without  complications .  Irradiation  of  the 
probe  with  high-energy  vuv  radiation  near  Lyman-a 
causes  the  energy-probe  output  to  decay  within  a  pe¬ 
riod  of  minutes,  becoming  less  sensitive  to  the  ArF 
and  the  dye  beams  as  well  as  to  the  vuv  beam.  The 
sensitivity  recovers  over  periods  of  10  min  to  several 
hours.  The  cause  of  this  detector-sensitivity  drop  is 


Fig.  1.  Experimental  apparatus  for  vuv  generation. 
PMT,  photomultiplier  tube. 
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not  currently  understood.  To  provide  a  stable  energy 
monitor,  we  also  use  a  solar-blind  KBr  photomultiplier 
with  a  Lyman-0?  filter  to  measure  vuv  light  scattered 
from  the  probe  at  right  angles  to  the  beam.  This  pho¬ 
tomultiplier —filter  combination  is  insensitive  to  ArF 
and  dye  light,  and  signal-to-background  ratios  of  up 
to  100  are  obtained.  We  estimate  the  accuracy  of  our 
energy  measurements  at  10%  for  direct  probe  readings 
and  30%  with  the  photomultiplier  and  the  filter. 

Without  phase  matching,  mixing  in  H2  produces  2 
orders  of  magnitude  more  energy  than  mixing  in  Kr 
near  Lyman- <2.  However,  H2  does  not  produce  vuv  ra¬ 
diation  well  at  Lyman-a  because  of  production  of  H 
atoms,  which  absorb  the  vuv.  In  addition,  resonant 
contributions  to  the  refractive  index  from  H  atoms 
produce  large  variations  in  phase  mismatch,  causing 
large  intensity  fluctuations  with  varying  vuv  wave¬ 
length.  Through  manipulation  of  experimental  pa¬ 
rameters  such  as  H2  pressure,  laser  powers,  laser  pulse 
synchronization,  and  ArF  wavelength,  we  were  not  able 
to  reduce  H  atom  production  sufficiently  to  allow  good 
conversion  at  Lyman-n.  For  this  reason  we  investi¬ 
gated  the  use  of  Kr. 

One  can  enhance  difference-frequency  mixing  in  Kr 
by  tuning  the  ArF  laser  to  two-photon  resonances  with 
the  6p[5/2, 2],  6p[3/2, 2],  and  6p[l/2, 0]  states.  Tran¬ 
sitions  to  6p[5/2,2]  and  6p[l/2,0]  lie  in  the  wings  of 
the  ArF  tuning  curve,  so  we  chose  the  6p[3/2,2]  tran¬ 
sition.  We  purged  the  ArF  beam  path  with  argon  to 
minimize  losses  owing  to  oxygen  Schumann-Runge 
absorptions  at  this  wavelength.8  Generation  of  vuv  ra¬ 
diation  in  pure  Kr  is  poor  because  of  the  relatively  large 
phase  mismatch  in  Kr.  The  pressure  dependence  of 
the  vuv  intensity  in  pure  Kr  is  shown  in  the  lower  part 
of  Fig.  2.  Peak  energies  of  ~0.2  pj  are  produced  at 
an  optimum  pressure  of  7  ±  3  mbar,  roughly  an  order 
of  magnitude  lower  than  the  optimum  pressure  in  H2. 
This  is  consistent  with  a  calculated  phase  mismatch  for 
Kr,  which  is  five  times  larger  than  for  H2  at  Lyman-n. 
Phase  mismatch  can  be  eliminated  by  the  addition  of 
a  positively  dispersive  gas  such  as  Ar  but  not  H2,  be¬ 
cause  H  atoms  are  produced  even  when  the  ArF  laser 
is  tuned  to  the  Kr  resonance. 

With  phase-matched  mixtures,  higher  partial  pres¬ 
sures  of  Kr  can  be  used.  Figure  2  shows  that  a  3.9:1 
Ar/Kr  mixture  with  260  mbar  of  Kr  yields  roughly 
2  orders  of  magnitude  of  improvement  over  pure  Kr. 
Phase-matched  vuv  generated  at  130  mbar  of  Kr  is 
shown  in  Fig.  3.  We  measured  as  much  as  7  pJ  of  en- 
ergy  with  130  mbar  of  Kr.  The  excimer  laser  has  a 
pulse  length  much  longer  than  that  of  the  dye  laser, 
so  we  use  the  dye  pulse  length  (5  ns  FWHM)  as  an  es¬ 
timate  of  the  vuv  pulse  length.  At  this  pulse  length 
we  produced  1.3  kW  of  power  at  Lyman- <2.  As  an  in¬ 
dication  of  the  brightness  of  our  source,  the  Lyman- 
a  radiation  can  be  observed  on  a  phosphor  screen 
with  the  room  lights  on.  The  highest  vuv  power  at 
Lyman- a,  that  was  previously  reported  apparently  cor¬ 
responds  to  250  W  of  power  external  to  the  mixing 
cell.4  This  previous  power  was  measured  indirectly 
through  comparison  of  dye  and  vuv  intensities  mea¬ 
sured  by  a  spectrometer  and  a  photomultiplier.  This 
measurement  approach  is  inaccurate  owing  to  a  variety 


of  systematic  errors,  such  as  comparing  beam  powers 
based  on  the  fraction  of  the  beams  sampled  by  the  spec¬ 
trometer  slit. 

Our  experiments  indicate  that  significantly  higher 
vuv  powers  are  possible  with  further  system  optimiza¬ 
tion.  Figure  2  shows  that  the  vuv  power  is  still  in¬ 
creasing  with  pressure.  Our  cell  was  not  designed  for 
operation  at  overpressure,  so  higher  total  pressures 
were  not  explored.  Further  improvement  is  also  pos¬ 
sible  by  using  larger  focal-spot  sizes  (longer-focal- 
length  lenses).  We  did  not  optimize  the  focal  length 
of  the  focusing  lens.  Because  the  generated  vuv  de¬ 
pends  on  the  path-integrated  laser  intensities,  it  varies 
with  the  incident  laser  powers,  not  with  the  peak  inten¬ 
sities  in  the  foci.10,11  Thus  larger  focal-spot  sizes  can 
reduce  saturation  effects  (which  increase  with  the  peak 
intensity)  without  sacrificing  vuv  power.  Indications 
of  saturation,  which  can  occur  as  a  result  of  dielectric 
breakdown,  multiphoton  ionization,  ground-state  de¬ 
pletion,  intensity-dependent  refractive-index  changes, 
or  other  processes,  are  present  in  our  measurements. 
Ideally,  the  vuv  power  scales  with  the  square  of  the 
density,  and  our  power  rises  more  slowly  (see  Fig.  2). 
To  gain  full  advantage  of  the  longer  conversion  path 
lengths  obtained  with  larger  focal-spot  sizes,  one  must 


Fig.  2.  Generated  vuv  as  a  function  of  Kr  partial  pressure. 
All  measurements  were  performed  at  Lyman-a. 


Wavelength  (nm) 

Fig.  3.  Absolute  vuv  energy  for  a  Kr  partial  pressure  of 
130  mbar.  7  jiJ  is  generated  at  Lyman-a. 
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Fig.  4.  Tuning  range  of  vuv  at  different  partial  pressures 
of  Kr.  Solid  curves,  experimentally  measured  tuning 
profiles;  dashed  curves,  fits  using  theoretical  profiles. 

ensure  that  the  beam  quality  is  quite  good,  which 
may  require  phase-conjugation  amplification  of  the 
ArF  laser.12  Breakdown  is  often  observed  in  the  fo¬ 
cal  region  at  higher  total  pressures.  Some  reduction 
in  breakdown  or  other  saturation  processes  may  be 
possible  by  use  of  a  buffer  gas  with  a  higher  refrac¬ 
tive  index  than  Ar,  as  the  mixture  used  here  is  only 
20%  Kr.  For  operation  at  high  pressures,  absorption 
and  (or)  dispersion  by  rare-gas  dimers  may  become 
important.1  Because  the  tuning  range  at  a  fixed  gas 
mixture  is  inversely  proportional  to  the  laser  confocal 
parameters  and  the  gas  pressure,  increased  power  may 
come  at  the  cost  of  reduced  tuning  range.  However, 
the  tuning  range  may  be  increased  by  using  a  gas  mix¬ 
ture  that  has  a  flatter  dispersion  than  Ar/Kr. 

The  spectral  tuning  ranges  for  three  Kr  partial  pres¬ 
sures  are  shown  as  solid  curves  in  Fig.  4.  As  expected, 
the  tuning  range  decreases  (approximately)  linearly 
with  increasing  pressure.  These  tuning  profiles  show 
a  strong  asymmetry,  too  strong  to  be  caused  by  a  cur¬ 
vature  in  the  phase  mismatch  as  a  function  of  wave¬ 
length.  Such  asymmetry  can  be  caused  when  ArF 
and  the  dye  beams  have  dissimilar  confocal  param¬ 
eters.  This  asymmetry  is  related  to  the  additional 
phase  change  that  occurs  when  the  beams  pass  through 
the  focus  (Guoy  effect).  When  the  confocal  parameters 
are  not  equal,  a  ripple  is  produced  on  the  phase  of  the 
driving  polarization  that  generates  the  vuv  radiation. 
Depending  on  whether  the  phase  of  the  driving  polar¬ 
ization  leads  or  lags  the  phase  of  the  vuv  wave  (i.e., 
whether  the  phase  mismatch  is  negative  or  positive), 
the  phase  ripple  can  be  compensated  for  to  a  greater  or 
lesser  extent,  leading  to  an  asymmetry. 

We  fitted  phase-matching  curves  to  the  profiles  in 
Fig.  4  based  on  the  expressions  in  Ref.  II.13  These 
fits  are  shown  as  dashed  curves  in  Fig.  4.  The  fitted 
confocal  parameters  for  the  ArF  beam  are  1.6  mm 
for  26  mbar  of  Kr  and  6  mm  for  the  two  higher 
pressures.  The  fitted  confocal  parameters  for  the  dye 


laser  are  four  times  smaller  than  the  ArF  confocal 
parameter  in  each  case.  If  we  use  the  low-pressure 
confocal  parameters  to  calculate  focal-spot  sizes,  we 
find  that  the  calculated  and  the  measured  spot  sizes 
agree  for  the  ArF  beam,  whereas  the  calculated  spot 
size  is  roughly  half  the  measured  spot  size  for  the 
dye  laser.  Our  lasers  do  not  have  diffraction-limited 
Gaussian  beams,  so  some  disagreement  in  these  values 
is  not  surprising.  The  fact  that  the  fitted  confocal 
parameters  are  not  the  same  for  all  pressures  may 
be  related  to  refractive-index  changes  induced  by  the 
lasers  or  to  other  saturation  effects.  The  dependence 
of  the  asymmetry  on  experimental  confocal  parameters 
means  that  the  direction  and  the  amount  of  asymmetry 
are  experimentally  dependent.  The  vuv  beam  profile 
also  depends  on  the  phase  mismatch  Ak  and  so  will 
vary  with  wavelength.10 

We  have  achieved  directly  measured  energies  of  7  p  J 
at  Lyman- a,  using  two-photon-resonant  difference- 
frequency  mixing  in  phase-matched  gas  mixtures  of 
Kr  and  Ar.  Significant  improvement  in  energy  is 
expected  with  further  system  improvement,  although 
this  may  come  at  the  expense  of  a  reduced  tuning  range 
for  a  fixed  gas  mixture. 
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Abstract 

Spatial  distributions  of  the  absolute  concentration  of  H,  CH,  C2,  and  C3  reactive  intermediate  species  have  been  measured  by 
laser-induced  fluorescence  (LIF)  in  the  freestream  plume  of  a  dc-arcjet.  Diamond  film  grows  from  this  reactive  mixture  on  a 
water-cooled  molybdenum  substrate  in  a  25  Torr  reactor.  The  temperature  of  the  gas  jet  is  determined  from  LIF  measurements 
of  the  rotational  distribution  of  CH  radicals  in  the  plume.  The  directed  velocity  of  the  gas  plume  is  derived  from  the  Doppler 
shift  of  LIF.  The  fraction  of  the  feedstock  hydrogen  dissociated  is  obtained  from  calorimetry.  The  atomic  hydrogen  concentration 
and  gas  temperature  are  used  to  constrain  a  model  of  the  gas  phase  chemistry  in  the  reactive  plume.  Concentration  measurements 
of  CH,  C2,  and  C3  provide  a  stringent  test  for  the  model  of  the  diamond  precursor  chemistry.  ©  1998  Elsevier  Science  S.A. 

Keywords:  Plasma  diagnostics;  Chemical  vapor  deposition;  Laser-induced  fluorescence 


L  Introduction 

The  chemistry  precursor  to  diamond  growth  in  con¬ 
ventional  diamond  reactors  scales  with  the  atomic 
hydrogen  concentration  [1-3].  The  atomic  hydrogen 
concentration  in  a  dc-arcjet  can  approach  [4]  half  of 
the  hydrogen  feedstock,  which  is  consistent  with  the 
high  quality  diamond  film  and  rapid  growth  rates 
observed  in  dc-arcjet  reactors  [5-8].  Chemical  models 
[9-12]  of  the  gas  phase  chemistry  in  the  arcjet  plume 
find  the  atomic  hydrogen  concentration  important  to 
predict  the  reactive  gas  flux  to  the  growing  diamond 
surface.  Dandy  and  Coltrin  [11]  find  the  atomic 
hydrogen  concentration  so  important  to  the  predictions 
of  their  model,  that  they  parameterize  the  results  as  a 
function  of  atomic  hydrogen  concentration.  The  atomic 
hydrogen  concentration  for  the  reactor  condition  studied 
here  has  been  determined  by  a  detailed  power  balance 
[4].  We  use  this  measured  atomic  hydrogen  concen¬ 
tration  as  an  input  parameter  for  the  model  of  the 
precursor  gas  phase  chemistry. 

The  reactive  environment  of  the  arcjet  plume  is  investi¬ 
gated  with  laser-induced  fluorescence  (LIF)  measure- 
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ments.  The  gas  temperature  is  derived  from  the 
rotational  temperature  measured  by  LIF  excitation  scans 
[13,14].  This  measured  temperature  constrains  the 
model  of  the  gas  phase  chemistry  in  the  freestream  of 
the  arcjet  plume.  The  directed  velocity  of  the  plume  is 
determined  from  the  Doppler  shift  of  LIF  between  NO 
seeded  into  the  plume  flow  and  NO  in  a  static  cell  [15]. 
The  directed  velocity  establishes  the  reaction  time  of  the 
reactive  mixture  of  gases  in  the  arcjet  freestream.  Both 
the  gas  temperature  and  the  directed  velocity  are  impor¬ 
tant  terms  in  the  calorimetry  to  determine  the  fraction 
of  the  feedstock  hydrogen  which  is  dissociated  in  the 
plume  [4], 

The  dc-arcjet  is  an  ideal  reactor  to  test  models  of  the 
chemical  mechanism  for  the  precursors  for  diamond 
growth;  the  rapid  growth  rate  in  such  reactors  suggests 
the  important  precursor  species  concentrations  are  large. 
In  addition,  the  gas  residence  time  in  a  dc-arcjet  is  well 
defined.  Common  diamond  reactors  such  as  hot  filament 
and  microwave  discharge  have  rather  small  input  feeds¬ 
tock  flow  with  a  long,  diffusion  controlled  residence 
time.  Small  quantities  of  impurities  or  minor  species  in 
the  feedstock  can  be  important  to  the  chemistry  of 
systems  with  long  residence  times  via  complex  and/or 
catalytic  chemical  reaction  schemes. 

The  feedstock  in  the  dc-arcjet  reactor  is  activated  in 
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the  discharge  and  expands  through  the  nozzle  into  a 
reactive  plume.  The  gas  plume  has  a  rapid  directed 
velocity  and  the  flight  time  to  the  boundary  layer  above 
the  substrate  is  well  defined.  The  transport  through  the 
thin  boundary  layer  is  diffusion  controlled.  In  this  paper, 
we  will  present  a  comparison  of  quantitative  LIF  meas¬ 
urements  [16]  of  CH,  C2,  and  C3  concentrations  in  the 
freestream  of  the  arcjet  plume,  with  the  predictions  of 
reactive  species  mole  fractions  from  a  detailed  chemical 
model  constrained  by  measured  gas  temperature  and 
atomic  hydrogen  concentration.  This  comparison  indi¬ 
cates  that  the  transport  in  the  arcjet  reactor  is  more 
complex  than  the  simple  picture  presented  above. 


2.  Experiment 

The  dc-arcjet  reactor  shown  in  Fig.  1  is  the  same  used 
for  previous  work  [4,13-16].  A  1.6  kW  arc  is  struck  in 
a  46:54  mixture  of  hydrogen  and  argon  at  a  pressure  of 
6  atmospheres.  The  effluent  from  this  arc  expands 
through  a  converging/diverging  nozzle  into  a  reactor 
maintained  at  25  Torr.  Methane  at  0.5%  of  the  hydrogen 
flow  is  injected  into  the  diverging  section  of  the  nozzle 
as  shown  in  Fig.  1.  For  the  velocity  experiments  0.1% 
NO  is  added  to  the  plume  via  the  methane  injector. 
Diamond  thin  film  grows  on  a  water-cooled  molyb¬ 
denum  substrate  inserted  into  the  plume  of  the  arcjet  at 
a  growth  rate  of  approximately  1  mm/min.  A  probe 
laser  beam  shown  in  Fig.  1  is  directed  across  the  arcjet 
plume,  the  subsequent  fluorescence  is  collected  normal 
to  the  plume  flow  and  the  laser  beam,  and  the  optical 
probe  volume  of  1  mm3  is  determined  by  the  intersection 
of  the  laser  beam  and  the  optical  signal  collection.  The 
arcjet  and  substrate  assembly  can  be  moved  with  respect 
to  optical  probe  volume  during  diamond  deposition  to 
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Fig.  1 .  Schematic  of  dc-arcjet  in  CVD  reactor. 


map  spatial  distributions  along  the  centerline  of  the 
plume  and  radially  across  the  cylindrically  symmetric 
plume. 

LIF  signals  from  H,  CH,  C2,  C3,  and  NO  are  excited 
with  light  from  a  variety  of  wavelength  tunable  laser 
sources,  and  the  details  of  the  excitation  and  optical 
collection  for  each  of  these  species  has  been  published 
in  detail  [14-17].  Atomic  hydrogen  is  excited  from  the 
IS  ground  state  to  the  3S  and  3D  excited  states  by  two 
photons  near  205  nm  and  the  subsequent  fluorescence 
to  the  2P  state  near  656  nm  is  observed  [17].  The  205  nm 
light  is  produced  by  frequency  tripling  a  Nd:YAG 
pumped  dye  laser.  CH  radicals  are  excited  in  the  0-0 
band  of  the  B-X  transition  near  390  nm  or  the  0-0 
band  of  the  A-X  transition  near  437  nm  with  light  from 
an  excimer  pumped  dye  laser,  and  the  subsequent  fluo¬ 
rescence  is  collected  in  the  same  band  excited  [14,16]. 
C2  radicals  are  excited  in  the  2-0  band  of  the  d-a  Swan 
system  and  fluorescence  collected  in  the  2-1  and  1-0 
bands  [16].  C3  radicals  are  excited  in  the  spectrally 
overlapped  vibrational  bands  of  the  A-X  transitions 
near  427  nm.  The  NO  LIF  is  excited  in  the  0-0  band  of 
the  A-X  transition  near  226  nm  with  frequency  doubled 
light  from  an  excimer  pumped  dye  laser  [15]. 


3.  Arcjet  plume  velocity 

A  trace  of  NO  is  seeded  into  the  arcjet  plume  and 
LIF  is  excited  by  a  laser  beam  directed  into  the  plume 
flow  through  a  small  hole  in  the  center  of  the  substrate. 
The  same  laser  beam  also  excites  LIF  from  NO  in  a 
static  cell.  The  directed  velocity  of  the  plume  produces 
a  Doppler  shift  in  the  NO  excitation  compared  with 
that  in  the  cell.  The  plume  velocity  determined  from 
these  Doppler  shift  measurements  is  shown  in  Fig.  2  as 
a  function  of  radial  distance  about  the  centerline  of  the 
plume.  The  peak  velocity  of  2.6  km/s  is  supersonic  at 
Mach  2.2  for  the  gas  mixture  of  0.26:0.27:0.47  of 


radial  position  [mm] 


Fig.  2.  Axial  velocity  vs  radial  position  at  the  nozzle  exit. 
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H:H2:Ar.  This  mixture  results  from  the  measured  [4] 
33%  dissociation  of  the  hydrogen  feedstock  with  an 
initial  mixture  of  46:54  hydrogen  and  argon. 


4.  Temperature 

The  gas  temperature  is  determined  from  LIF  measure¬ 
ments  of  the  rotational  distribution  of  CH  [14,16]. 
Because  the  CH  radicals  chemically  react  rapidly  with 
both  H  and  H2  in  the  arcjet  plume,  we  tested  the 
assumption  that  the  CH  rotational  distribution  is  equili¬ 
brated  with  the  gas  temperature.  We  compared  the  CH 
measurements  with  an  LIF  measurement  of  the  rota¬ 
tional  distribution  of  NO  seeded  into  the  flow  [16].  NO 
does  not  react  with  the  major  species  in  the  arcjet  plume 
(H,  H2,  and  Ar),  and  thus,  undergoes  hundreds  of 
collisions  to  insure  the  rotational  distribution  of  NO  is 
in  thermal  equilibrium  with  gas  temperature  in  the 
freestream  of  the  plume.  We  found  agreement  between 
the  NO  and  CH  data  and  conclude  that  the  CH  rota¬ 
tional  distribution  provides  a  good  measure  of  the  gas 
temperature  [16]. 

Fig.  3  shows  the  gas  temperature  on  the  centerline  of 
the  arcjet  plume  between  the  nozzle  exit  at  z  =  0  and  the 
substrate  at  z  —  38.2  mm.  The  gas  temperature  is  con¬ 
stant  in  the  freestream  and  there  is  nearly  a  1000  K 
temperature  increase  in  the  boundary  layer  above  the 
substrate.  This  temperature  rise  is  consistent  with  a 
shock  heating  from  the  impinging  supersonic  flow.  In 
the  shock,  both  temperature  and  pressure  will  rise 
abruptly. 


5.  Concentration  measurements 

Quantitative  concentration  of  H,  CH,  C2,  and  C3  are 
determined  in  the  reactive  arcjet  plume  from  LIF  meas- 
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Fig.  3.  Temperature  along  the  centerline  of  the  dc-arcjet  measured  by 
excitation  LIF  of  CH  B-X  (0,0)  rotational  transitions. 


urements  [14,16,17].  Quantitative  LIF  requires  knowl¬ 
edge  of  the  gas  temperature,  fluorescence  quantum  yield, 
spectroscopic  data  for  the  transition,  and  the  detection 
efficiency  for  the  optical  collection.  LIF  results  from 
optical  excitation  of  a  specific  quantum  state  of  an  atom 
or  molecule;  therefore,  gas  temperature  must  be  known 
to  correct  the  signal  for  the  Boltzmann  fraction  of  the 
population  in  the  specific  state  initially  excited.  The 
quantum  yield  is  the  probability  that  the  excited  state 
will  radiate  a  signal  photon,  including  all  the  non- 
radiative  loss  rates  for  the  excited  state.  In  the  arcjet 
plume,  LIF  from  H,  CH,  C2,  and  C3  must  compete  with 
collisional  quenching.  For  CH,  C2,  and  C3  the  quantum 
yield  is  directly  determined  by  comparison  of  the  mea¬ 
sured  fluorescence  decay  time  with  the  radiative  lifetime 
[14,16].  From  time-resolved  LIF  measurements,  the 
quantum  yield  is  obtained  for  different  spatial  positions 
in  the  plume.  For  atomic  hydrogen,  the  quantum  yield 
measured  at  very  low  reactor  pressures  compares  well 
with  collisional  quenching  rates  determined  from  a 
model  of  the  major  species  concentrations  and  literature 
values  of  the  quenching  rate  constants.  The  model 
predictions  are  then  extrapolated  to  higher  reactor  pres¬ 
sures,  where  direct  measurements  of  the  time-resolved 
fluorescence  with  nanosecond  laser  excitation  are  not 
possible  [17]. 

Different  strategies  are  used  to  determine  the  quantita¬ 
tive  scale  for  each  species  [14,16,17].  For  CH  and  C2 
there  is  a  complete  spectroscopic  data  base,  a  single 
quantum  state  can  be  excited  within  our  laser  and 
Doppler  bandwidths,  and  the  laser  excitation  and  fluo¬ 
rescence  rates  are  well  known.  The  C3  absorption 
spectrum  is  quasi-continuous  at  this  temperature  and 
spectral  region.  The  probability  for  laser  absorption  is 
obtained  from  measurements  of  temperature-dependent 
optical  absorption  [18].  Two-photon  absorption 
depends  critically  on  the  mode  statistics  of  the  laser; 
therefore,  we  chose  to  calibrate  the  quantum  yield- 
corrected  LIF  measurements  of  relative  concentration 
of  atomic  hydrogen  with  calorimetry  [4,17].  The  optical 
detection  efficiency  for  CH,  C2,  and  C3  LIF  is  calibrated 
by  Rayleigh  scattering. 

The  variation  of  the  absolute  number  density  for  H, 
CH,  C2,  and  C3  measured  along  the  centerline  in  the 
freestream  of  the  arcjet  plume  is  presented  in  Fig.  4. 
Atomic  hydrogen,  CH,  and  C2  radicals  peak  near  the 
nozzle  and  decline  in  the  freestream  by  approximately 
30%.  C3  radical  concentrations  increase  by  approxi¬ 
mately  a  factor  of  four  between  the  nozzle  and  the 
boundary  layer. 


6.  Chemical  reaction  model 

The  gas  phase  chemistry  in  the  freestream  is  modeled 
as  a  one-dimensional  (ID)  flow  with  axial  diffusion 
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Fig.  4.  Comparison  between  concentrations  of  gas-phase  species  in  the 
arcjet  determined  by  a  ID  model  (lines)  and  calibrated  UF  (symbols). 


using  the  Sandia  premix  and  Chemkin  codes  [19].  We 
have  constructed  our  own  chemical  mechanism  (E.A. 
Brinkman  and  G.P.  Smith,  pers.  comm.)  consisting  of 
Cl  and  C2  chemistry  from  the  GRI-Mech  [20]  mecha¬ 
nism  for  methane  combustion  and  estimates  of  the 
C3-C6  chemistry.  The  model  calculation  is  constrained 
by  the  measured  temperature,  gas  flow,  gas  velocity, 
and  initial  atomic  hydrogen  concentration.  The  solid 
lines  in  Fig.  4  are  predictions  of  species  concentrations 
in  the  freestream  as  a  function  of  distance  from  the 
nozzle  without  a  substrate  present.  The  methane  is 
injected  inside  the  diverging  nozzle  at  the  position  z  = 
—  8  mm.  LIF  measurements  are  possible  only  after  the 
plume  exits  the  nozzle. 

The  deviation  between  calculated  and  measured 
atomic  hydrogen  is  a  direct  measure  of  the  radial 
diffusion,  which  is  not  included  in  this  model.  The 
measured  CH  radical  concentration  is  in  remarkable 
agreement  with  the  model  prediction.  Near  the  nozzle, 
we  find  C2  and  C3  present  in  concentrations  more  than 
two  orders  of  magnitude  larger  than  predicted  by  the 
model.  This  deviation  near  the  nozzle  is  far  larger  than 
any  estimated  uncertainty,  and  strongly  suggests  that 
the  gas  re-circulation  in  the  arcjet  reactor  (see  Fig.  1 )  is 
an  important  contribution  to  the  chemical  composition 
of  the  arcjet  plume.  The  radial  distribution  of  the  C3 
forms  a  hollow  cylinder  about  the  center  of  the  plume 
[16,21].  This  radial  distribution  is  another  indicator  that 
hydrocarbon  radicals  re-circulate  from  the  shock  struc¬ 
ture  in  front  of  the  substrate  and  are  entrained  in  the 


plume  near  the  nozzle.  Such  a  fluid  flow  provides  a 
much  longer  residence  time  for  the  precursor  chemistry. 


7.  Conclusions 

LIF  measurements  of  gas  temperature  and  directed 
velocity  are  used  to  constrain  a  model  of  the  gas  phase 
chemistry  of  the  freestream  plume  in  a  dc-arcjet  reactor. 
The  model  predictions  for  H,  CH,  C2,  and  C3  are 
compared  with  absolute  radical  concentrations  mea¬ 
sured  with  LIF.  We  find  remarkable  agreement  between 
measured  and  predicted  CH  concentrations.  The  C2 
and  C3  concentrations  near  the  nozzle  indicate  that 
re-circulation  must  be  included  in  the  transport  model. 
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Spatially  resolved  optical  emission  spectroscopy  is  used  to  investigate  excited  species  in  a  dc-arcjet 
diamond  depositing  reactor.  Temperature  measurements  indicate  a  cold  plasma  with  electrons, 
excited  states,  and  gas  in  nonthermal  equilibrium.  The  H,  C,  C2,  and  CH  excited  state  number 
densities  decrease  exponentially  with  the  distance  from  the  nozzle  and  have  a  pronounced  increase 
in  the  shock  structure  above  the  substrate.  The  H  emission  increases  throughout  the  boundary  layer 
to  the  substrate  surface,  whereas  emission  from  other  species  has  a  maximum  in  the  boundary  layer 
and  then  decreases  again  towards  the  substrate.  The  reconstructed  radial  distribution  of  excited  state 
concentrations  are  Gaussian,  with  the  C  and  C2  distributions  broader  than  the  H  and  CH  ones.  The 
optical  emission  is  calibrated  with  either  Rayleigh  scattering  or  laser-induced  fluorescence  to  furnish 
absolute  number  densities.  We  find  all  the  excited  species  to  be  present  in  concentrations  two  or 
more  orders  of  magnitude  smaller  than  the  corresponding  ground  states  measured  in  the  same 
reactor  and  conditions.  We  find  that  C2(d-a)  emission  intensity  correlates  well  with  laser-induced 
fluorescence  measurements  of  C2(a)  concentration  in  the  arcjet  plume.  Ground  state  concentrations 
of  the  other  species  do  not  vary  as  their  emission  intensity  except  near  the  substrate,  where  the 
variations  of  CH(A-X),  CH(fl-X),  and  C2(d-a)  emission  intensities  are  good  monitors  of  the 
corresponding  concentration  changes.  ©  1998  American  Vacuum  Society. 
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spectroscopy  (OES)  is  another  nonintrusive  method,  with 
quite  simple  data  collection  and  straightforward  automation. 
Unfortunately  optical  emission  arises  only  from  atoms  and 
molecules  in  electronically  excited  states,  and  interpretation 
of  the  optical  emission  signal  requires  an  understanding  of 
the  link  between  excited  and  ground  state  concentrations.  In 
spite  of  these  limitations,  OES  is  widely  used  for  real-time 
monitoring  in  plasma  processing. 

The  diamond  film  quality  varies  dramatically  with  depo¬ 
sition  method  and  conditions  and  a  large  body  of  research 
has  been  undertaken  to  develop  a  simple  relationship  be¬ 
tween  OES  from  a  species,  the  ground  concentration  of  that 
species,  and  the  quality  of  deposited  diamond  film.  Klein- 
Douwel  et  al ,22  found  a  correlation  between  C2  emission  and 
the  diamond  film  quality  for  material  grown  in  a  oxyacety- 
lene  flame.  In  a  microwave  plasma  reactor,  Lang  et  al.23 
found  the  diamond  film  quality  to  correspond  to  H  atom 
actinometry24  measurements.  Gruen  et  al.25  found  OES  sig¬ 
nals  from  C2  emission  in  a  microwave  plasma  correlate  with 
the  growth  rate,  but  not  the  emission  from  either  H  or  CH.  In 
a  hot  filament  reactor,  Cui  et  al26  observed  emission  from  H 
and  CH+  relate  to  high  quality  diamond,  and  CH  emission  to 
indicate  the  presence  of  amorphous  carbon  films. 

In  CVD  dc-arcjet  reactors,  the  diamond  film  grows-on 
temperature  controlled  substrates  from  hydrogen  arcjet 
plumes  with  a  trace  (<  1  %)  of  hydrocarbon.  Although,  op¬ 
tical  emission  from  the  reactive  plumes  of  dc  arcjets  has  been 
studied,  all  of  the  published  work  focuses  on  temperature 
determination  from  the  rotational  and  vibrational  distribu¬ 
tions  of  the  excited  molecules.  Raiche  and  Jeffries27  report 
emission  from  CH(A-X),  CH(B-X),  CH(C-X),  C3(A-X), 


I.  INTRODUCTION 

Low-pressure  dc-arcjet  reactors  are  attractive  sources  of 
chemically  reactive  gas  flows  for  spray  coating  and  chemical 
vapor  deposition  (CVD).  Arcjet  reactors  have  demonstrated 
diamond  growth  with  rates  approaching  1  mm/h,1  at  least  ten 
times  faster  than  in  other  reactors.2  This  rapid  growth  rate  is 
a  consequence  of  the  large  fraction  of  dissociation  of  the 
hydrogen  feedstock  gases  in  the  reactive  gas  flow  from  the 
arcjet  discharge.3 

Gas  mixtures  containing  argon,  hydrogen,  and  methane 
are  commonly  used  to  produce  diamond  films  in  plasma  re¬ 
actors.  The  chemistry  in  most  diamond  CVD  is  dominated 
by  charge  neutral  reactions  of  ground  state  radicals  and  at¬ 
oms  formed  in  the  discharge,  and  there  are  models  predicting 
C  atoms  and  CH3  radicals4,5  as  well  as  C2  and  C2H26’7  as  the 
dominate  precursors  to  diamond  growth.  In  order  to  charac¬ 
terize  the  reactive  environment  for  diamond  CVD,  experi¬ 
ments  to  determine  the  ground  state  concentration  of  H,8-10 
c,11  C2, 12,13  C3,13  CH,13-17  CH3,11,14,15,18-20  and  C2H220  have 
been  undertaken  using  nonintrusive  optical  diagnostics  tech¬ 
niques  in  a  wide  variety  of  diamond  producing  reactors.  Al¬ 
though  this  research  has  increased  our  understanding  of  dia¬ 
mond  CVD,21  the  methodology  is  of  limited  use  for  real-time 
industrial  process  control,  as  such  careful,  quantitative,  laser- 
based  measurements  require  expensive,  high  maintenance 
equipment  with  highly  trained  technicians.  Optical  emission 
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and  C2(d-a )  in  a  dc  arcjet  (1  kw,  220  Torr),  and  do  not  find 
thermalization  between  excited  and  ground  states  of  C2. 
Reeve  and  Weimer28  observe  CH (A-X)  and  C2(d-a)  emis¬ 
sion,  and  find  quite  different  vibrational  and  rotational  tem¬ 
peratures  in  CH(A)  and  C2(d).  The  OES  study  by  Cuber- 
tafon  et  al.29  reports  a  rotational  temperature  of  CH(Z?)  near 
2700  K  along  the  plume  of  their  arcjet  (3-5  kW,  —50  Torr) 
which  is  nearly  three  times  less  than  their  C2  vibrational 
temperature  of  7500  K.  They  also  find  quite  different  vibra¬ 
tional  and  rotational  temperatures  from  the  excited  C2(d). 
Brinkman  et  al.30  find  nearly  the  same  C2(d)  rotational  and 
vibrational  temperatures  as  Cubertafon  in  their  arcjet  (1.6 
kW,  —25  Torr).  Analysis  of  the  spatial  distributions  of  the 
emission  intensities  led  them  to  suggest  that  the  excited 
C2(d)  was  formed  by  chemical  reactions  and  CH(A)  by 
electron  impact. 

In  this  work,  we  present  spectroscopic  observations  from 
a  dc-arcjet  reactor  during  diamond  CVD  using  the  same  op¬ 
erating  conditions30  where  diamond  film  grows  at  50  /zm/h 
on  a  water  cooled  molybdenum  substrate.  The  diamond  film 
is  well  faceted  with  primarily  (111)  morphology,  and  the 
Raman  spectrum  of  the  film  shows  no  trace  of  graphite  con¬ 
tent.  We  have  characterized  the  arcjet  plume  in  previous 
laser-based  and  probe  measurements.  Using  quantitative 
laser-induced  fluorescence  (LIF),  we  determined  the  number 
densities  and  spatial  distributions  of  CH,  C2,  C3 
radicals13  31,32  and  H  atoms  in  the  ground  state.10  The  gas 
temperature13  30  distribution  was  previously  determined  from 
LIF  measurements  of  rotational  distributions  of  ground  state 
molecules.  The  flow  velocity10  was  measured  from  the  Dop¬ 
pler  shift  of  LIF  from  NO  seeded  into  the  flow.  Langmuir 
probe  was  used  to  determine  the  electron  temperature  and 
ion  density.33  We  report  here  the  quantitative,  spatially  re¬ 
solved  optical  emission  of  H,  C,  C2,  and  CH.  From  these 
measurements,  we  determine  excited  state  concentration, 
correlate  the  excited  state  concentration  distribution  with  the 
respective  ground  state  distribution  and  discuss  the  chemical 
roles  and  excitation  mechanisms. 

II.  EXPERIMENTAL  METHOD 

The  CVD  reactor,  diamond  deposition,  and  details  about 
laser  measurements  have  been  described  pre¬ 
viously. 10,1 3,30,33,34  An  arc  (1.6  kW,  12  A,  140  V)  is  struck  in 
a  mixture  of  H2/Ar  (0.9:1  at  6.9  slm  total  flow  and  six  atmo¬ 
spheres  pressure);  the  effluent  from  this  arc  expands  through 
a  converging/diverging  nozzle  into  a  low-pressure  reactor 
(typically  25  Torr).  Methane  (—0.5%  of  the  H2  flow)  is  in¬ 
jected  in  the  flow  in  the  diverging  section  of  the  nozzle.  The 
expanding  gas  forms  a  luminous  plume  approximately  1  cm 
in  diameter  which  impinges  on  a  water  cooled  molybdenum 
substrate  38.2  mm  from  the  nozzle.  Diamond  grows  at  nearly 
1  yum/min  with  a  well  faceted  (1 1 1)  morphology  with  a  Ra¬ 
man  spectrum  free  from  graphitic  features.10 

The  optical  emission  from  the  arcjet  plume  is  collected 
with  a  lens  //6,  spatially  filtered  with  a  1  mm  aperture,  and 
focused  with  another  lens  identical  to  the  first  one  into  a 
monochromator  selected  to  fit  the  requirements  of  the  mea¬ 


surement.  For  rotationally  resolved  scans,  a  Heath  mono¬ 
chromator  is  used  (ffl ,  0.35  m  focal  length,  250  or  500  nm 
blaze  gratings).  Laser-based  measurements  require  a  detec¬ 
tion  bandpass  with  equal  response  for  any  rotational  level  in 
the  selected  vibrational  band;  here  we  construct  a  trapezoidal 
bandpass  with  30  nm  full  width  at  half  maximum  (FWHM) 
from  a  Bausch-Lomb  high  intensity  monochromator  (f/6 , 
0.2  m  focal  length,  500  nm  blaze  grating)  using  a  narrow 
front  slit  and  a  wide  back  slit.  The  signal  is  detected  by  a 
1P28  photomultiplier,  preamplified  (X10),  and  averaged 
with  a  boxcar  integrator  (Stanford  Research  Systems  SRS- 
250).  During  the  calibration  procedures  of  the  emission  sig¬ 
nal,  the  voltage  in  the  photomultiplier  power  supply  and  box¬ 
car  settings  were  kept  identical  to  those  used  for  the  optical 
emission  measurements.  A  SC  Technology  Plasma  Chemis¬ 
try  Monitor  401,  consisting  of  a  spectrograph  with  a  gated, 
intensified,  photodiode  array  of  512  elements,  is  used  as  an 
optical  multichannel  analyzer  to  obtain  emission  spectra  si¬ 
multaneously  in  the  range  from  200  to  1000  nm  with  a  4  nm 
bandwidth.  The  relative  spectral  response  of  this  instrument 
is  calibrated  using  standard  lamps  (Optronics  Laboratories), 
tungsten  in  the  visible  and  deuterium  in  the  ultraviolet.  The 
two  different  lamps  provide  an  overlapping  region  from  300 
to  350  nm  where  we  find  good  (<5%)  agreement  between 
the  two  independent  calibration  measurements. 

The  absolute  intensity  of  the  optical  emission  is  calibrated 
by  both  Rayleigh  scattering  and  laser-induced  fluorescence 
of  a  known  gas  concentration  and  composition.  A  XeCl  ex- 
cimer  laser  (Lambda  Physik  EMG  103)  pumping  a  dye  laser 
(Lambda  FL  2002)  is  the  source  of  vertically  polarized  laser 
light  to  excite  CH(R-Y)  or  provide  a  Rayleigh  scattering 
source.  The  14  ns  laser  pulses  have  a  spectral  bandwidth  of 
0.20  ±0.02  cm-1  at  430  nm,  measured  with  a  monitor  etalon. 
Some  measurements  around  315  nm,  were  made  with  a  fre¬ 
quency  doubled  PDL-1  dye  laser  (bandwidth  —  1  cm”1, 
pulse  length  —7  ns),  pumped  by  a  Quanta  Ray  GCR-3 
Nd:yttrium-aluminum-gamet  (YAG)  laser.  A  high  preci¬ 
sion  energy  meter  (Rj-7200,  Laser  Precision  Corp.)  monitors 
the  laser  pulse  energies  from  10  nJ  to  1  mJ. 

III.  QUANTITATIVE  OPTICAL  EMISSION 

The  optical  emission  signal  of  a  given  electronic  transi¬ 
tion  between  an  excited  state  u  and  the  ground  state  g  during 
a  time  interval  tg ,  is  equal  to 

n 

‘^OES~i^ug*  wexc*  Yfl’  ^ -€‘7]'tgi  (1) 

where  Aug  is  the  Einstein  emission  coefficient  of  the  ob¬ 
served  vibrational  band  (s-1),  nexc  is  the  steady  state  number 
density  of  excited  species  (cm-3),  and  is  the  observed 
volume.  The  remaining  factors  are  given  by  the  optics  and 
electronics  in  the  detection  of  the  fluorescence;  where  fl  is 
the  solid  angle,  e  is  the  transmission  efficiency  of  the  optics, 
and  7]  the  photoelectric  conversion.  The  emission  from  the 
total  production  of  electronically  excited  species,  nnascent,  is 
reduced  by  collisional  quenching  and  predissociation;  the 
steady  state  concentration  of  excited  species  is 
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n exc  =  ti nascent^5 »  where  3>=A/(A  +  <2  +  P)  is  the  quantum 
yield,  Q  the  quenching  rate  (s_1),  and  P  the  predissociation 
rate  (s-1).  We  note  the  nexc  is  the  excited  state  number  den¬ 
sity  which  is  measured  in  absorption. 

One  of  the  main  problems  related  to  OES  absolute  mea¬ 
surements  is  the  limited  spatial  resolution  of  optical  emis¬ 
sion.  When  a  single  lens  or  a  fiber  optic  are  utilized  to  collect 
the  emission  from  the  reactor,  the  volume  Vfl  is  defined  by 
the  lens  /  number  or  the  fiber  optic  acceptance  angle,  and  in 
both  cases  gives  poor  results.  Using  optical  spatial  filters,  we 
obtain  a  line  of  sight  measurement  of  the  emission.  However, 
only  when  the  spatial  distribution  of  emitters  is  homoge¬ 
neous  and  the  size  of  the  plasma  is  known  ( d ,  in  cm),  the 
number  density  in  the  line  of  sight  volume  is 

H*xc=d  Jo  n^x)dx  =  A^d  ’  (2) 

where  Nexc  is  the  number  of  emitters,  Vfl= And,  Afl=  7rr2  is 
the  collection  cross  section  given  by  the  spatial  filter,  d  is  the 
plasma  length.  The  plasma  size  and  emitter  distribution  are 
usually  unknown,  and  we  can  only  obtain  integrated  column 
densities  (emitters/cm2).  This  is  the  same  problem  that  oc¬ 
curs  for  absorption  measurements,  and  it  is  solved  for  geom¬ 
etry  with  cylindrical  symmetry  through  Abel’s  inversion. 
The  arcjet  plume  is  cylindrically  symmetric  and  the  spatial 
distribution  of  the  different  species  can  be  reconstructed 
from  the  set  of  measurements  taken  along  the  perpendicular 
axes  to  the  line  of  sight.  There  are  many  proposed  methods 
for  the  computation  of  Abel’s  inversion.  In  this  work,  we 
have  followed  the  method  proposed  by  Buie  et  al., 35  which 
fits  the  data  to  a  polynomial  with  even  powers.  The  param¬ 
eters  are  the  input  to  an  analytically  derived  form  of  the  Abel 
inversion.  We  will  find  below  that  all  the  species  have 
Gaussian  radial  distributions  in  the  arcjet  plume.  The  peak 
number  densities  are  fit  by  Eq.  (2)  with  the  parameter  d 
equal  to  the  FWHM  of  the  reconstructed  Gaussian  distribu¬ 
tion. 

The  instrumental  factors  (ft, 6, 77)  in  Eq.  (1)  can  be  deter¬ 
mined  by  other  optical  measurements  which  include  the 
same  set  of  parameters.  For  example,  Rayleigh  scattering  of 
a  gas  with  the  same  optical  collection  given  by36 
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where  n  is  the  number  density  of  the  gas,  v  is  the  frequency 
of  the  scattered  photon  (cm-1),  (daldfl)  is  the  averaged 
cross  section,  EL  is  the  laser  energy  (J) ,  V  is  the  collection 
volume,  and  A  L  is  the  laser  cross  section,  set  by  apertures 
approximately  equal  to  Afl.  From  the  plot  of  Rayleigh  scat¬ 
tering  signal  versus  the  product  of  laser  energy  and  pressure, 
we  determine  the  experimental  factor  £lerjVlAL. 

The  optical  emission  can  also  be  calibrated  by  comparison 
to  the  signal  from  a  known  excited  state  number  density.  LIF 
can  be  used  if  the  ground  state  number  density  of  one  of  the 
plasma  species  is  known  as  well.  With  low  laser  spectral 
irradiance  to  insure  linear  excitation  and  full  time  integration 


of  the  fluorescence  signal,  the  relationship  between  the  de¬ 
tected  LIF  signal  SLIF  and  the  molecular  number  density  n0 
is  given  by  ’ 
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fB ,  the  Boltzmann  factor,  is  the  fraction  of  molecules  in  the 
probed  quantum  state  and  is  calculated  from  the  temperature 
determined  by  LIF  excitation  scans.  B  is  the  absorption  co¬ 
efficient  for  the  rotational  transition  excited  (cm2  J-1  s-1), 
El  is  the  laser  energy  ( J),  T  is  the  lineshape  overlap  integral, 
Av  is  the  laser  bandwidth  (cm-1),  the  ratio  reff/r0  is  the 
fluorescence  quantum  yield  with  reff  as  the  effective  lifetime 
including  collisional  quenching,  and  r0  is  the  radiative  life¬ 
time.  Ffl  is  the  fraction  of  the  fluorescence  collected  in  the 
spectral  bandwidth  of  the  detector. 

The  LIF  signal  is  produced  only  in  the  laser  probed  vol¬ 
ume  and  can  be  related  to  the  OES  volume.  Assuming  the 
LIF  measurement  excites  a  rotational  line  in  the  same  elec¬ 
tronic  transition  observed  in  optical  emission 
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The  analogous  expression  for  Rayleigh  calibration 

^Ray  =  _n_  477 -Tp 

oes  «exc  h  2rdhc-v-Ffl 

at  a  given  laser  energy  EL ,  and  with  integration  time  tg  long 
enough  to  integrate  the  whole  laser-induced  fluorescence  sig¬ 
nal  or  the  Rayleigh  signal,  we  can  estimate  the  steady  state 
number  density  of  the  corresponding  excited  state  (nexc). 


IV.  RESULTS  AND  DISCUSSION 

The  wavelength  resolved  optical  emission  spectra  from 
the  arcjet  plume  varies  dramatically  as  the  feedstock  gas  is 
changed.  When  the  feedstock  is  pure  argon,  we  see  strong 
emission  features  from  the  Ar  and  weak  lines  from  Ar+. 
With  a  mixture  of  argon  and  hydrogen,  we  see  a  few  of  the 
strongest  Ar  lines,  but  the  dominant  emission  is  from  the 
Balmer  series  of  atomic  hydrogen  (Fig.  1);  we  are  not  able  to 
detect  the  vacuum  ultraviolet  Lyman  a  to  compare  its  emis¬ 
sion  intensity.  Adding  just  0.5%  methane  to  the  flow  again 
alters  the  emission  from  the  plasma  plume,  we  find  signifi¬ 
cant  atomic  emission  from  H  atoms  in  the  Balmer  series,  C 
atoms,  and  very  weak  Ar  I  lines.  There  is  strong  emission 
from  CH(A-X),  CH (B-X),  CH(C-X),  and  C2  (d-o)  bands. 
The  remaining  structure  in  the  spectra  is  assigned  to  mol¬ 
ecules  as  well.  Emission  from  two  other  C2  electronic  sys¬ 
tems  is  assigned:  the  C2(D-X)  Mulliken  bands  at  232  nm 
and  the  C2(C-A)  Deslandres-D’Azambuja  bands  near  360 
nm. 

The  molecular  emission  becomes  more  dominant  down¬ 
stream  from  the  expansion  nozzle;  Fig.  2  shows  the  spectrum 
at  z  =  20  mm  which  is  roughly  halfway  between  the  nozzle 
and  the  substrate.  The  optical  emission  spectrum  also  varies 
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Fig.  1.  Upper  panel:  Optical  emission  spectrum  (4  nm  FWHM),  taken  near 
the  nozzle  in  a  dc  arcjet  with  25  Torr  reactor  pressure,  and  55%  Ar  and  45% 
H2  feedstock  at  1.6  kW.  Lower  panel:  Electron  temperature  from  atomic 
hydrogen  Balmer  series. 


Fig.  2.  Optical  emission  spectra  (4  nm  FWHM),  taken  20  mm  downstream 
from  the  arcjet  nozzle  with  the  standard  conditions.  Upper  panel.  Center  of 
plasma  plume.  Lower  panel.  At  a  radius  of  5.5  mm  from  the  center. 


with  radial  position  in  the  plume;  the  upper  panel  shows  the 
spectrum  from  the  center  of  the  plume  which  is  dominated 
by  CH  emission  with  only  a  small  underlying  broadband 
visible  background.  Figure  2  lower  panel  is  taken  at  a  radius 
of  5.5  mm  where  the  total  emission  intensity  scale  has  been 
multiplied  by  20.  Here  we  see  the  dominant  structured  emis¬ 
sion  is  from  the  C2  Swan  bands  ( d-a ).  In  addition,  the 
broadband  visible  background  contains  most  of  the  emission 
intensity.  The  C3  radical  ground  state  has  its  maximum  in¬ 
tensity  at  large  radius  in  the  cylindrically  symmetric  plume,32 
similarly  at  large  radius  near  the  edge  of  the  plume,  we  iden¬ 
tify  C3  emission  as  the  405  nm  feature  of  the  C3(A-X)  (000)- 
(000)  band,  and  we  infer  C3  as  the  source  of  the  less  struc¬ 
tured  background  emission  between  330  and  500  nm.  In 
addition  to  C3  emission,  the  visible  background  likely  in¬ 
cludes  contributions  from  other  hydrocarbon  fragment  like 
C2H37  or  carbon  clusters  C4  and  C5  by  analogy  to  observa¬ 
tions  in  neon  matrices.38  This  speculation  is  further  rein¬ 
forced  by  ab  initio  calculations;38,39  of  both  radicals  which 
find  electronic  bands  in  the  visible  region,  near  500  nm. 
However,  C4  and  C5  have  not  been  spectroscopically  identi¬ 


fied  in  gas  phase  yet.  In  addition,  we  observe  light  from 
OH(A-X)  and  NH(A-X)  at  308  and  336  nm,  respectively; 
these  molecules  are  present  in  the  plume  either  from  impu¬ 
rities  in  the  feedstock  gases  or  minor  leaks  in  the  vacuum 
system.  The  only  ion  observed  is  emission  assigned  to 
CH+(A-X)  as  it  can  be  seen  in  Fig.  3.  Details  of  the  most 
important  species  and  transitions  identified  can  be  found  in 
the  Table  I. 

A.  Temperatures 

The  energy  distribution  in  the  plasma  is  an  important 
characterization  parameter.  We  determine  excited  state  tem¬ 
peratures  from  the  rotational  (or  vibrational)  distribution  of 
molecular  optical  emission,  the  ground  state  (gas)  tempera¬ 
ture  from  LIF  measurements  of  the  ground  state  rotational 
distribution,  and  electron  temperature  from  the  atomic  hy¬ 
drogen  excited  state  distribution.  Atomic  hydrogen  emission 
lines  are  a  good  measure  of  the  electron  temperature  when 
electron  impact  is  the  excitation  mechanism  and  the  excita¬ 
tion  energies  are  all  in  the  tail  of  the  electron  energy  distri¬ 
bution.  Assuming  the  various  atomic  hydrogen  excited  states 


J.  Vac.  Sci.  Technol.  A,  Vol.  16,  No.  2,  Mar/Apr  1998 


401 


Luque  et  a!.:  Excited  state  density  distributions  of  H,  C,  C2,  and  CH 


401 


stock  mixture  of  only  Ar  and  H  we  obtain  a  temperature  of 
1.25±0.25eV  near  the  nozzle  (Fig.  1,  bottom),  in  good 
agreement  with  earlier  Langmuir  probe  measurements  of 
electron  temperature.33  Although  some  slight  systematic  de¬ 
viations  from  a  linear  Boltzmann  plot  are  observed  in  Fig.  1, 
these  are  minor  and  can  be  attributed  to  either  quenching 
effects  or  nonequilibrium  features  of  the  plasma.41  The  addi¬ 
tion  of  methane  lowers  the  temperature  to  —  1  eV 
(-11  600  K). 

We  examine  the  rotational  distribution  in  the  excited 
CH(Z?)  and  CH(C)  and  compare  the  results  to  our  previous 
results30  for  excited  CH(A)  and  C2(d).  The  emission  is  ana¬ 
lyzed  by  spectral  simulation  with  LIFB  ASE  42  including  both 
predissociation  and  collisional  quenching.  These  corrections 
are  important  to  determine  the  nascent  population  distribu¬ 
tions,  especially  in  the  case  of  the  CH(C).  From  time  re¬ 
solved  LIF,  we  measure  the  quenching  of  the  CH(Z?)31  to  be 
3850  3900  3950  4000  4050  4100  1  jjs~1  Torr"'1,  and  found  no  rotational  dependence  of  this 

Wavelength  (A)  collisional  quenching  rate.  There  is  no  data  for  the  CH (C) 

state,  but  there  is  evidence  of  its  quenching  with  H2  that  is 

Fig.  3.  Optical  emission  spectrum  from  the  plasma  plume  on  the  centerline  faster  than  the  quenching  of  the  CH  A  and  B  States,  and 

with  z=20  mm  (0.1  mm  FWHM)  showing  emission  from  CH(fl),  CH +(A),  we  estimate  a  quenching  rate  of  -2  /xs  Torr  \  Predisso- 

C2<C),  and  atomic  C  and  H  in  the  lower  panel.  Upper  panel  shows  Simula-  ciation  is  Strong  in  CH(£,  u'=0,  N'>  15),  CH  (B  u'=l, 

tion  of  CH(fl)  with  a  rotational  temperature  of  2500  K,  and  vibrational  of  >6)  45  an(J  fov  aJ]  level§  of  CH(C),46’47  varying  with  TO- 

il0°  K  tational  and  vibrational  quantum  numbers. 

Emission  from  CH(Z?)  is  prominent  in  the  wavelength 
have  similar  excitation  cross  sections  and  collisional  quench-  dispersed  optical  emission  in  the  region  384-410  nm  shown 

ing  rates,  the  atomic  hydrogen  emission  intensity  from  each  in  the  lower  panel  of  Fig.  3.  This  emission  is  collected  20 

level  is  proportional  to  the  number  density  of  electrons  at  mm  downstream  of  the  nozzle  from  the  center  of  the  plume, 

that  energy.  The  lower  panel  in  Fig.  1  shows  the  Boltzmann  The  upper  panel  is  a  simulation  of  the  CH (B)  (0,0)  and  (1,1) 

plot  of  the  atomic  hydrogen  Balmer  series.  The  relationship  bands  with  a  rotational  temperature  of  2500  K.  We  find  the 

between  the  fine  intensities  /,  and  temperature  T  is  excited  CH(£,  v’=0)  state  is  well  fit  by  rotational  tempera¬ 

tures  between  2300  and  2700  K.  These  temperatures  are  sub- 
(7)  stantially  lower  to  the  ones  measured30  in  CH(A)  and  C2(d), 
and  in  good  agreement  with  measurements  of  CH(Z?)  of  Cu- 
where  is  the  degeneracy  of  the  state  z,  Aiy-  the  emission  bertafon  et  al29  in  a  similar  reactor.  In  addition,  the  rota- 

coefficient,  and  Et  the  energy.40  When  we  calculate  an  emis-  tional  distributions  of  CH(tf )  produced  by  electron  impact 

sion  temperature  for  the  Ha,  H^,  Hy,  and  H*  with  a  feed-  on  CH4  and  C2H2  near  the  appearance  potential  (-14  eV) 


Table  I.  Summary  of  the  identified  species,  electronic  transitions,  and  their  dominant  band  or  line  positions. 
Also,  a  comparison  of  the  ground  and  excited  states  number  densities  in  the  middle  of  the  arcjet  plume  at  20 
mm  from  the  nozzle. 


Species 

Transition 

Line 

positions 

(nm) 

Emission 

coefficient 

Cs-') 

Excited  state 
number  density 
(cm-3) 

Ground  state 
number  density 
(cm-3) 

H 

3n  — 2n 

656 

4.4X  107  (Ref.  40) 

107 

3  X 1016  (Ref.  10) 

4n  —  2n 

486 

8.4X  106 

107 

C 

lP°-lD 

193 

2.4  X  10s  (Ref.  40) 

5X106 

not  measured 

lP°-xS 

247 

3.4X  107 

CH 

A-X(2 A-2n)  Au  =  0 

431 

1.8X  106  (Ref.  47) 

1.2X  109 

b-xC  x--2n>  Aw=  — 1,0 

387 

3  X  I06  (Ref.  47) 

3X108 

3.5  X1012  (Ref.  13) 

c-x(2s+-2n)  Av=o 

314 

1 X 107  (Ref.  64) 

6X107 

CH+ 

A-X('n-'2+)  Aw  =  1,0 

396,  423 

2.5  X 106  (Ref.  65) 

-3X107 

not  measured 

c2 

rf-a(3nr3ny)  Aw =1,0-1 

465,  480,  515 

1 X 107  (Ref.  13) 

2X108 

3  X 1010  (Ref.  13) 

D-X('2*-'X*)  Au=0 

231 

5.5  X 107  (Ref.  66) 

-5X105 

c-A('n,-1nj  Au=i,o,-i 

360,  385,  410 

3.2X  107  (Ref.  66) 

-2X106 

C3 

A-x('nu-'2;) 

405 

5X106  (Ref.  13) 

<108 

-3X1012  (Ref.  13) 

Ei 

=-lf+c • 
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Wavelength  (A) 


Fig.  4.  Optical  emission  from  the  plasma  plume  on  the  centerline  with  z 
=  20  mm  (0.2  nm  FWHM),  comparing  emission  from  CH(C)  in  the  upper 
trace  with  a  simulation  at  5500  K  in  the  lower  trace.  Note  the  experimental 
spectrum  also  shows  emission  from  OH  (A  -X)  near  308  nm,  and  its  simu¬ 
lation  is  set  to  7’=  4000  K. 


have  been  reported  in  the  literature,  with  temperatures  in 
v '  =  0  of  —2600  and  —2300  K,  respectively,48,49  reinforcing 
the  hypothesis  that  CH(A)  and  CH (B)  are  products  of  elec¬ 
tron  impact  dissociation. 

Figure  4  compares  the  wavelength  resolved  emission  for 
CH(C)  with  our  simulation  using  a  rotational  and  vibrational 
temperature  of  5500  K.  All  of  the  CH  electronic  states  have 
different  rotational  temperatures,  5500  K  for  CH(C),  2500  K 
for  CH (B),  and  3500  K  for  CH(A).  Figure  5  shows  the  LIF 
excitation  spectrum  of  CH(C-X)  taken  in  the  same  position 
in  the  plume  which  compares  well  with  the  simulation  at  a 
temperature  of  2200  K  in  good  agreement  with  previous 
ground  state  LIF  measurements13,30,31  with  the  same  reactor 
conditions.  These  earlier  experiments  found  the  LIF  rota¬ 
tional  distribution  of  ground  state  CH  to  be  a  good  measure 
of  the  gas  temperature.  Thus,  we  find  the  dc-arcjet  plume  is  a 
nonequilibrium  plasma  with  different  thermal  distributions 
for  the  electrons,  molecular  excited  states,  and  the  bulk  gas. 

B.  Excited  state  number  densities 

The  number  densities  are  proportional  to  the  ratio  be¬ 
tween  band  areas  and  emission  coefficients  in  spectrally  cor¬ 
rected  scans  like  the  ones  in  Fig.  2.  The  optical  detection  is 
calibrated  to  determine  quantitative  excited  state  number 
densities  from  the  optical  emission.  We  use  Rayleigh  scat¬ 
tering  to  calibrate  the  optical  emission  signal  collection.  The 
scattering  was  produced  with  a  laser  wavelength  of  388  nm, 
in  the  same  spectral  region  where  the  CH(Z?-X)  emission 
occurs.  The  calibration  of  the  emission  signal  is  performed 
20  mm  downstream  from  the  nozzle,  with  the  plasma  jet  off 


Fig.  5.  Laser-induced  fluorescence  spectrum  exciting  CH(C)-CH(X)  in 
the  plasma  plume  on  the  centerline  with  z  —  20  mm  in  the  lower  trace  (0.012 
nm  FWHM).  The  lower  trace  is  a  simulation  spectrum  at  2200  K. 


and  the  chamber  filled  with  argon  to  a  pressure  between  10 
and  100  Torr.  The  number  density  of  CH  X  is  known  from 
earlier  measurements,31  and  LIF  on  /?i(10)  CH(2?-X)  (0,0) 
was  used  as  a  secondary  calibration  standard.  We  found  both 
methods  in  agreement  within  15%,  but  the  LIF  calibration  is 
less  reliable  because  it  is  an  indirect  measurement,  which 
requires  calibration  as  well.  For  the  excited  state  number 
density  measurements  reported  here,  the  errors  are  domi¬ 
nated  by  the  mismatch  of  the  calibrated  detection  volume 
and  the  actual  emission  volume,  and  we  expect  overall  25%- 
50%  uncertainty. 

Once  the  ratio  between  optical  emission  and  calibration 
signal  is  obtained,  the  number  density  of  CH(i?)  is  corrected 
by  the  radial  profile  to  account  for  line  of  sight  integration  of 
the  emission  signal.  Table  I  shows  the  number  densities  of 
excited  state  species  after  normalizing  to  the  CH  B  value,  the 
emission  coefficients,  and  ground  state  number  densities 
whenever  they  are  available.  The  excited  state  number  den¬ 
sities  estimated  are  in  the  105-109  mol/cm3  range,  and  are 
all  several  orders  of  magnitude  less  than  the  respective  | 
ground  state  number  density.  Those  values  are  reasonable  in 
a  system  where  the  chemistry  is  dominated  by  ground  state 
species. 

The  highest  ratio  of  excited  to  ground  state  concentration 
is  found  for  C2(rf)/C2(a)— 0.01.  Even  higher  ratios,  with 
population  inversion  between  C2(d)  and  C2(a)  have  been 
measured  by  absorption  spectroscopy  in  an  arc  plasma  dur¬ 
ing  the  deposition  of  amorphous  hydrogenated  carbon 
coatings.67  Excited  C2(d)  appears  readily  in  these  systems, 
but  it  is  not  accompanied  by  a  large  amount  of  C2  ground 
state;  for  example,  Gruen  et  al  ,25  observe  bright  emission 
during  hydrogen  free  diamond  growth  in  a  microwave 
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Fig.  6.  Upper  panel:  Radial  distribution  of  excited  state  number  density  at 
z-  20  mm  in  the  freestream  of  the  plasma  plume.  Note  the  C2(d )  structure 
is  wider  than  that  for  CH(£).  Lower  panel:  Radial  distribution  from  Ref.  13 
of  ground  states  from  LIF  measurements  in  the  same  arcjet  conditions. 

plasma  reactor  and  they  infer  C2  is  an  important  reactant  for 
diamond  growth.  In  the  dc  arcjet,  a  strong  chemilumines¬ 
cence  does  not  correlate  with  a  large  amount  of  C2,  and  the 
total  number  density  of  this  radical  is  too  small  to  support  the 
observed  diamond  growth  rate. 

1 .  Spatial  distributions 

Radial  distributions  of  the  excited  state  concentrations  of 
H,  C,  CH  (A),  CH  (£),  CH  (C),  C2(rf),  and  C2(D)  states 
have  been  obtained  by  the  Abel  transform,  and  measure¬ 
ments  of  CH(A)  and  C2(d)  at  z  =  20  mm  from  the  nozzle  are 
shown  in  the  upper  panel  of  Fig.  6.  We  find  the  recovered 
distributions  to  have  a  Gaussian-like  shape,  and  the  C2  dis¬ 
tribution  is  wider  than  for  CH(A)  (7  mm  vs  5  mm  FWHM), 
which  is  consistent  with  observed  differences  between  their 
respective  ground  states.  Nonetheless,  the  excited  states  are 
observed  in  a  region  of  the  plume  significantly  narrower  than 
for  the  ground  states,13  where  we  find  the  radial  distributions 
with  a  FWHM  greater  than  10  mm  (Fig.  6,  lower  panel).  We 
find  two  distinct  radial  distributions  with  excited  H,  CH(A), 
CH (B),  and  CH (C)  have  narrow  distributions,  and  C, 
C2(d),  and  C2(D)  have  wide  distributions. 

Figure  7  is  a  contour  map  of  the  cylindrically  symmetric 
excited  state  distribution  for  CH(R);  the  top  and  center  of 
the  figure  is  the  nozzle  exit  and  the  substrate  is  at  38.2  mm. 


Radial  distance  (mm) 

Fig.  7.  Contour  plot  of  the  spatially  resolved  concentrations  of  CH(fl)  be¬ 
tween  the  nozzle  and  the  substrate;  note  the  increased  emission  in  the 
boundary  layer  above  the  substrate.  Number  densities  are  in  units  of 
108  molecules/cm3. 


Along  the  centerline  in  the  plume,  we  see  the  number  density 
decreases  noticeably,  before  rising  again  just  above  the  sub¬ 
strate.  However,  the  radial  distribution  does  not  change  be¬ 
tween  the  nozzle  exit  and  the  shock  structure  just  above  the 
substrate.  Figure  8  shows  the  distribution  of  excited  C2(d), 
where  the  vertical  gradient  is  not  so  large,  and  again  the 
radial  distribution  does  not  change  between  the  nozzle  and 
the  substrate.  Emissions  from  H,  C,  CH(A),  and  CH(C)  also 
have  maxima  on  the  centerline  of  the  plume.  However,  the 
emission  from  C3  and  larger  hydrocarbons  does  not  peak  on 
the  centerline  as  seen  in  the  radial  variation  of  the  broadband 
visible  background  in  Fig.  2. 

There  is  a  large  variation  of  the  optical  emission  intensity 
as  a  function  of  distance  from  the  nozzle.  The  upper  panel  of 
Fig.  9  plots  the  excited  state  number  density  along  the  cen¬ 
terline  of  the  plume  versus  distance  from  the  nozzle  for 
atomic  carbon  (lP),  CH(A),  atomic  hydrogen  (n  =  3), 
C2(D),  and  C2(d);  CH(tf),  and  CH(C)  are  omitted  because 
they  have  a  similar  behavior  to  CH(A).  The  lower  panel  of 
Fig.  9  contrasts  the  excited  state  concentrations  on  the  cen¬ 
terline  of  the  arcjet  plume  with  the  ground  state  concentra¬ 
tions  measured  earlier  by  laser-induced  fluorescence.13,52 
Note  the  excited  state  number  densities  are  several  orders  of 
magnitude  smaller  than  the  respective  ground  states.  Both 
atomic  hydrogen  and  CH  radical  ground  state  concentrations 
are  nearly  constant  with  distance  from  the  nozzle  while  their 
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Radial  distance  (mm) 

Fig.  8.  Contour  plot  of  the  spatially  resolved  concentrations  for  C2(d). 
Number  densities  are  in  units  of  108  molecules/cm3. 


excited  state  concentrations  decline  rapidly.  Whereas  the 
C2(d)  and  C2(a)  state  concentrations  have  quite  similar 
variation  with  distance  from  the  nozzle. 

The  excited  state  number  density  for  all  of  the  species 
decreases  with  distance  from  the  nozzle  until  the  shock 
heated  boundary  layer  just  above  the  substrate  surface.  The 
plume  velocity  is  greater  than  the  sound  speed,  therefore  a 
shock  is  created  when  this  supersonic  plume  impinges  on  the 
substrate.  The  shock,  called  a  normal  shock  because  the  flow 
velocity  is  normal  to  the  substrate,  creates  a  pressure  and 
temperature  rise  in  the  boundary  layer.  The  excited  state 
number  density  for  all  the  species  except  atomic  hydrogen 
rises  and  shows  a  peak  between  the  shock  and  the  substrate. 
Excited  atomic  hydrogen  also  rises  in  the  shock  heated 
boundary  layer;  however,  instead  of  peaking  in  the  boundary 
layer  and  decreasing  near  the  surface,  the  excited  H  atoms 
continue  to  increase  with  a  maximum  at  the  surface. 

In  the  freestream,  the  region  of  the  plume  flow  before  the 
shock,  the  excited  state  number  density  exponentially  de¬ 
creases,  as  evident  from  the  linear  semilogarithmic  plots  in 
Fig.  9;  all  of  the  excited  state  concentrations  decrease  by 
more  than  a  factor  of  100  except  the  factor  of  8  for  the 
decrease  of  C2(d).  The  variation  in  excited  state  number 
density  is  much  more  dramatic  than  the  decline  of  less  than  a 
factor  of  2  observed  for  the  ground  state  CH31  and  atomic 
hydrogen.52  The  data  in  Fig.  9  show  that  each  of  the  different 
excited  species  has  a  different  exponential  decay  with  dis¬ 


Distance  from  nozzle  (mm) 

Fig.  9.  Centerline  number  density  of  CH(/4),  C2{d ),  C 2(D),  C,  and  3 n 
hydrogen  atoms  as  a  function  of  distance  from  the  nozzle  with  a  substrate  at 
38.2  mm  in  the  upper  panel  contrasted  with  LIF  measurements  of  CH(X), 
C '?(<*),  and  ground  state  hydrogen  atoms  from  Refs.  13  and  52  in  the  lower 
panel. 

tance,  or  in  the  semilogarithmic  plots  a  different  slope  /?: 

n»c(d)  =  nexc(0)-e~fid.  (8) 

We  can  use  the  rapid  decrease  in  optical  emission  with  dis¬ 
tance  from  the  nozzle  and  the  variation  of  the  /3  to  unravel 
excitation  mechanism  of  the  optical  emission. 

The  radiating  species  are  excited  in  the  plume  and  emit 
their  optical  emission  without  moving  very  far  from  position 
of  excitation.  The  distance  from  the  nozzle  in  Fig.  9  is  re¬ 
lated  to  time  via  the  plume  velocity.  For  the  conditions  re¬ 
ported  here,  we  have  measured10  the  plume  velocity  to  be 
—  2.5  km/s;  therefore,  the  transit  time  between  the  exit  of  the 
nozzle  and  the  boundary  layer  above  the  substrate  is 
— 15  jj,s.  This  time  is  very  long  compared  to  the  excited  state 
lifetimes;  we  measured  collisionally  shortened  effective  life¬ 
times  of  CH (A,B)  and  C2(d)  to  be  30-45  ns.  Therefore, 
even  at  the  supersonic  plume  velocity,  the  excited  species  are 
produced  and  deactivated  within  a  distance  of  only  0.1  mm. 
There  are  four  excitation  mechanisms  which  might  produce 
excited  atoms  and  molecules  in  the  arcjet  plume:  radiative 
transfer  from  the  discharge  inside  the  nozzle,  thermal  exci- 
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tation,  electron  impact,  and  chemiluminescent  reactions.  We 
can  quickly  eliminate  radiative  transfer  and  thermal  excita¬ 
tion  as  important  sources  of  optical  emission  from  the  arcjet 
plume. 

Radiative  transfer  from  optical  emission  in  the  discharge 
region  inside  the  nozzle  can  produce  an  exponential  decay  of 
excited  state  population  with  distance  from  the  nozzle  only  if 
the  lower-state  densities  in  the  optical  transition  are  large 
enough  for  the  plume  to  be  optically  thick.  Previous  quanti¬ 
tative  LIF  measurements13  show  that  CH(X)  and  C2(a)  have 
number  densities  of  1012  and  1010  cm-3,  respectively.  Thus, 
the  CH  and  C2  transitions  can  not  be  optically  thick.  Ground 
state  atomic  hydrogen  has  a  density  of  greater  than 
1016  cm-3  and  the  plume  is  optically  thick  for  Lyman  tran¬ 
sitions  which  terminate  on  the  ground  state.  However,  we 
find  the  number  density  of  H(n  =  3)  to  be  in  the  range  of 
106- 109  cm-3  and  Balmer  transitions  are  optically  thin  in 
the  arcjet  plume.  Atomic  carbon  is  similarly  optically  thin. 
Therefore,  the  H,  C,  CH,  and  C2  excitation  mechanisms  in 
the  arcjet  plume  are  not  radiative  transfer  from  the  discharge 
inside  the  nozzle. 

Thermal  excitation  is  not  important  in  the  plume  either. 
The  temperature  is  nearly  constant  between  the  nozzle  and 
the  boundary  layer;30,31  therefore,  thermal  excitation  would 
not  produce  the  exponential  decline  in  excited  state  popula¬ 
tion  with  distance  from  the  nozzle.  In  addition,  thermal  equi¬ 
librium  calculations  with  temperature  of  2200  K  predict  a 
ratio  of  10~7- 10“5  for  CH  (A  or  B)  or  C2(d)  and  the  respec¬ 
tive  ground  states,  much  smaller  than  the  ratio  of 
10_2-10"3  observed. 

Chemiluminescent  chemical  reactions  are  an  important 
excitation  mechanism  in  many  reacting  flows;  for  example, 
chemiluminescent  reactions  are  the  primary  production  path¬ 
way  for  the  formation  of  excited  CH(A),  CH(Z?),  and  C2(d) 
in  flames.  However,  the  number  density  and  the  spatial  dis¬ 
tribution  of  most  of  the  excited  species  in  the  arcjet  plume 
are  inconsistent  with  a  chemiluminescent  reaction  mecha¬ 
nism. 

Most  of  the  combustion  reactions  which  produce  emission 
from  CH  require  oxygen  containing  reactants.50  In  the  arcjet 
plume,  oxygen  is  only  present  as  an  impurity;  thus,  these 
chemiluminescent  reactions  are  unlikely  to  produce  the 
bright  emission  from  CH(A),  CH(£),  and  CH(C)  observed 
from  the  arcjet  plume.  The  possible  chemiluminescent  reac¬ 
tion  excitation  mechanisms  are  also  constrained  by  the  large 
decrease  in  optical  emission  intensity  with  distance  from  the 
nozzle.  Fast  chemical  reaction  rate  coefficients  are 
1011  cm3/s,  and  the  plume  feedstock  consists  of  only  hydro¬ 
gen,  argon,  and  less  than  1%  hydrocarbon.  Therefore,  the 
chemical  composition  of  the  freestream  is  in  nearly  frozen, 
and  we  do  not  expect  large  gradients  in  concentrations  along 
the  centerline  of  the  freestream  of  the  plume.5’51  Our  mea¬ 
surements  of  the  ground  state  concentrations  of  the  reactive 
intermediates  H,52  CH,  C2,  and  C313  find  only  modest  con¬ 
centration  variations  between  the  nozzle  and  the  boundary 
layer  (Fig.  9,  lower  panel);  H,  CH,  and  C2  decrease  less  than 
a  factor  of  2  and  C3  increases  approximately  a  factor  of  2. 


Except  for  C2(d),  all  of  the  excited  state  concentrations  in 
Fig.  9  decrease  at  least  a  factor  of  100  between  the  nozzle 
exit  and  the  boundary  layer.  Thus,  chemiluminescent  reac¬ 
tion  can  be  the  dominant  excitation  pathway  only  for  the  C2 
swan  band  ( d-a )  emission.  The  most  likely  reaction  to  ex¬ 
plain  the  chemiluminescence  of  C2,  as  suggested  by 
Brinkman30  and  Wakisaka,53  is  the  formation  of  C2(d)  by 
three  body  recombination: 

C+C+M-+C2(J)  +  M. 

This  is  consistent  with  the  predictions,5’51  that  in  the  arcjet 
plume  with  its  large  molefraction  of  atomic  hydrogen, 
atomic  carbon  is  the  most  abundant  carbon  species. 

Electron  impact  is  an  important  excitation  mechanism  for 
much  of  the  optical  emission  in  the  plume.  The  electron 
number  density  is  the  only  plasma  parameter  with  a  signifi¬ 
cant  variation  along  the  vertical  axis,  with 
— 1012  electron/cm3  at  12  mm  downstream  from  the  exit 
plane  of  the  nozzle  which  decreases  by  an  order  of  magni¬ 
tude  at  32  mm  downstream.33  Measurements  and  models  of 
electron  densities  in  dc  arcjets  show  the  electron  number 
density  is  much  higher  inside  the  nozzle  than  in  the  plume.41 
Thus,  the  electron  density  declines  rapidly  with  distance 
from  the  nozzle  exit  similar  to  the  decrease  in  optical  emis¬ 
sion,  and  the  high  electron  density  in  the  nozzle  region  con¬ 
tributes  to  enhance  the  dissociation  of  methane. 

Electron  impact  dissociation  of  methane  can  have  excited 
CH  products  when  the  electron  energy  is  greater  than  the 
appearance  potential  of  approximately  14  eV 

*-(~14  eV) 4- CH4— >CH(A ,  B ,  C  or  X)  +  H2+H+e~. 

At  higher  electron  energies,  excited  H,  C,  and  CH+  species 
as  products  are  possible.49,54"57  Hydrocarbon  fragments  can 
also  be  dissociated  or  excited  by  direct  electron  impact.58  For 
example,  the  CH  radical  can  be  directly  excited  by  low  en¬ 
ergy  electrons 

*-(~3_4  eV)  +  CH— CH (A,B,C)  +  e~. 

The  most  abundant  two-carbon  hydrocarbon  in  the  arcjet 
plume  was  measured  and  predicted  to  be  C2H2.59  Its  electron 
impact  dissociation  has  been  studied48,56,57,60  and  the  elec¬ 
tronically  exited  fragments  detected,  in  addition  to  the  ones 
observed  from  methane,  are  the  C2(d ),  C2(C),  and 
C2(D). 56,60  We  find  emission  from  all  of  these  radicals  in  the 
dc-arcjet  plume.  Collisions  of  metastable  argon  atoms  could 
contribute  to  the  excitation  mechanism.  However,  argon 
metastables  are  very  rapidly  deactivated  by  molecular 
hydrogen,61  and  we  estimate  a  lifetime  in  the  plume  of  less 
than  0.5  /j&.  Even  at  the  fast  directed  plume  flow  of  2.5  km/s, 
the  argon  metastables  are  collisionally  deactivated  in  the 
nozzle  or  in  the  first  few  mm  of  the  freestream  plume. 

Assuming  the  number  density  of  electrons  in  the  plume 
decreases  with  distance  with  a  constant  electron  temperature, 
we  should  observe  the  same  gradient  in  the  excited  state 
emission  for  all  the  species  produced  by  electronic  impact. 
However,  we  find  a  different  slope  (3  for  each  excited  spe¬ 
cies,  and  Fig.  10  shows  the  correlation  of  /3  with  the  excita¬ 
tion  energy.  We  see  a  remarkable  linear  correlation  between 
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Fig.  10.  Variation  of  the  exponential  factor  p  vs  excitation  energy  of  the 
corresponding  radical  or  atom. 


P  and  the  excitation  energy,  with  only  C2(d )  falling  off  this 
line.  We  find  a  similar  correlation  between  the  appearance 
potential  for  electron  impact  dissociation  to  produce  excited 
CH(A),  CH(Z?),  and  CH(C)  from  methane,  excited  C2(D) 
from  C2H2,  and  excited  H  from  H2.  Thus,  we  could  explain 
the  spatial  distribution  of  optical  emission  for  all  the  species 
except  C2( d)  with  an  electron  impact  excitation  mechanism 
if  the  high  energy  tail  of  the  distribution  depletes  with  dis¬ 
tance  from  the  nozzle  in  addition  to  the  drop  in  electron 
density.  Such  a  decrement  is  consistent  with  the  Langmuir 
probe  observation33  that  the  high  energy  component  to  the 
electron  distribution  was  greatly  reduced  by  the  addition  of 
hydrocarbon  molecules  into  the  plume. 


2.  Effect  of  methane  addition 

The  amount  of  methane  in  the  gas  mixture  is  one  of  the 
most  critical  parameters  controlling  the  growth  rate  and  qual¬ 
ity  of  the  deposited  diamond  films.  Methane  addition  also 
reduces  the  electron  temperature  and  number  density  in  the 
plume  of  the  dc-arcjet  plasma.  We  found  feedstock  mixtures 
of  Ar  and  H2  without  added  methane  had  bimodal  electron 
energy  distributions  characterized  by  temperatures  ~2  and 
~  8  eV,  respectively,33  when  methane  is  added,  the  high  en¬ 
ergy  component  disappears.  This  observation  is  common  to 
other  plasma  reactors  with  H2  and  CH4  mixtures.62  The  ad¬ 
dition  of  methane,  a  polyatomic  collision  partner  with  many 
low-energy  vibrational  modes  and  overtones,  provides  effi¬ 
cient  relaxation  of  the  high  energy  electrons  in  the  plasma 
plume.63  Thus,  as  seen  in  Fig.  11,  the  addition  of  methane 
significantly  reduces  the  H  atom  optical  emission,  consistent 
with  an  electron  impact  excitation  mechanism. 

Figure  1 1  shows  excited  state  densities  deduced  from  the 
optical  emission  near  the  nozzle,  in  the  arcjet  plume,  and  in 
the  shock  heated  boundary  layer  as  a  function  of  methane 
addition.  Near  the  nozzle,  the  addition  of  methane  produces  a 
strong  decrease  in  the  H  atom  concentration,  and  the  rise  of 
emission  from  carbon  containing  species  from  electronic  im¬ 
pact  dissociation  of  hydrocarbons;  above  1%  of  methane  ad¬ 
dition,  this  increase  slows.  In  the  arcjet  plume,  the  atomic 
hydrogen  emission  decreases  with  methane  addition;  the 
emission  increases  for  the  C  species  up  to  0.5%  of  added 
methane,  but  beyond  this  value,  emission  from  all  the  C 
species  except  C2(d)  decreases.  The  C 2(d)  behavior  is  con¬ 
sistent  with  the  hypothesis  of  a  chemiluminescent  production 
mechanism  from  atomic  carbon  recombination.  The  rise  and 
fall  of  the  other  species  is  related  to  the  depletion  of  electron 


CH4/H2  Flow  ratio(%)  CH4/H2  Flow  ratio  (%)  CH4/H2  Flow  ratio  (%) 


Fig.  11.  Variation  of  excited  state  number  density  vs  methane  addition  near  the  nozzle  (z  =  0),  in  the  freestream  of  the  plume  (z= 20  mm),  and  in  the  shock  ] 
heated  boundary  layer  (z  =  5  7 . 5  mm).  | 
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energy  and  number  density,  similar  to  the  reduction  in  emis¬ 
sion  with  distance  seen  in  Fig.  9  for  constant  methane  addi¬ 
tion.  Finally,  in  the  shock  heated  boundary  layer,  where  the 
temperature  and  pressure  are  higher  than  in  the  other  regions, 
there  are  more  paths  for  the  radical  production.  Excitation  of 
H,  C,  and  CH(C)  by  electron  impact  is  still  likely,  but 
CH(A),  CH(Z?),  and  C2(d)  are  proportional  to  the  methane 
flow  and  more  likely  to  be  produced  by  chemiluminescent 
reactions.  At  this  point,  we  do  not  know  which  reactions  are 
responsible  for  CH(A)  and  CH(#)  production,  nor  why  the 
CH(C)  is  not  chemically  produced. 

V.  SUMMARY  AND  CONCLUSIONS 

In  the  arcjet  plume,  the  predominant  production  mecha¬ 
nism  of  electronically  excited  H*,  C*,  CH(A),  CH(i?), 
CH(C),  CH+(A),  C2(C),  and  C2(D)  is  electronic  impact 
excitation,  whereas  C2(d)  is  the  product  of  chemilumines¬ 
cent  reactions.  The  excited  state  populations  are  several  or¬ 
ders  of  magnitude  lower  than  the  ground  state,  with  peak 
concentrations  near  the  nozzle,  where  the  electron  density  is 
higher.  The  spatial  variation  of  excited  state  number  density 
is  quite  different  than  the  spatial  distribution  of  the  ground 
state  atoms  and  radicals.  The  number  density  of  excited 
states  is  sufficiently  low  that  excited  state  chemistry  is  un¬ 
likely  to  play  any  important  role  in  the  overall  chemistry  of 
diamond  CVD. 
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Atomic  hydrogen  concentration  in  a  diamond  depositing  dc  arcjet 
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The  fraction  of  hydrogen  dissociated  in  the  plume  of  a  dc  arcjet  used  for  diamond  deposition  is 
determined  by  calorimetry  to  be  0.33  ±0.12.  A  dc  arc  is  struck  in  a  mixture  of  argon  and  hydrogen 
at  90  psi  and  the  effluent  is  expanded  through  a  converging/diverging  nozzle  into  a  reactor 
maintained  at  25  Torr.  Methane  (<  1%)  is  added  to  the  luminous  gas  plume  in  the  diverging  nozzle. 
This  supersonic  jet  impinges  on  a  water  cooled  molybdenum  substrate,  and  diamond  thin  film  grows 
from  the  reactive  mixture.  The  electrical  power  input  of  the  arcjet  (1.6  kW)  is  balanced  by  the  power 
losses  due  to  cooling  of  the  nozzle,  enthalpy  change  in  the  gas,  ionization  of  the  gas,  dissociation 
of  H2,  and  the  directed  velocity  of  the  gas  phase.  The  gas  temperature  is  determined  by  linear 
laser-induced  fluorescence  (LIF)  measurements  of  several  rotational  lines  of  NO  seeded  to  the  gas 
plume.  The  velocity  of  the  gas  plume  is  obtained  via  the  Doppler  shift  between  LIF  signals 
measured  simultaneously  in  a  stationary  reference  cell  and  in  the  arcjet.  The  atomic  hydrogen 
fraction  increases  as  a  function  of  the  feedstock  H2/Ar  ratio  and  with  increasing  input  power  to  the 
arcjet.  ©  1997  American  Institute  of  Physics.  [S0021-8979(97)04112-l] 


I.  INTRODUCTION 

Thin  films  of  diamond  deposited  by  chemical  vapor 
deposition  (CVD)  are  used  in  a  variety  of  industrial  applica¬ 
tions  such  as  thermal  management,  wear  resistant  tool  coat¬ 
ings,  and  optical  windows.  Over  the  past  ten  years,  a  great 
deal  of  research  has  been  undertaken  to  understand  the  basic 
chemistry  and  physics  of  thin  film  diamond,  and  results  of 
this  research  were  recently  summarized  by  Goodwin  and 
Butler.1  CVD  diamond  is  grown  in  a  wide  variety  of  reactor 
designs  with  a  plethora  of  feedstock  gas  recipes;  however, 
nearly  all  these  techniques  can  be  reduced  to  the  growth  of 
CVD  diamond  from  a  small  amount  of  hydrocarbon  added  to 
a  hydrogen  plasma. 

A  mechanistic  view  of  CVD  diamond  growth  has 
evolved  from  the  insight  of  many  researchers.1-7  Atomic  hy¬ 
drogen  plays  a  crucial  role  in  this  mechanism;  gas-surface 
abstraction  reactions  of  atomic  hydrogen  produce  active  radi¬ 
cal  sites  on  the  growing  diamond  surface,  carbon  containing 
radicals  bind  to  these  sites,  and  recombination  reactions  ter¬ 
minate  the  surface  carbon  to  maintain  sp 3  bonding.  Using 
this  mechanistic  view  of  diamond  growth,  Goodwin8  and 
Bachmann  et  al?  developed  scaling  laws  for  growth  rates 
and  material  quality  that  depend  strongly  on  the  amount  of 
atomic  hydrogen. 

High  quality  diamond  film  can  be  grown  in  a  dc-arcjet 
reactor  with  deposition  rates  of  up  to  100  /zm/h,  among  the 
fastest  of  any  growth  rates  reported.10-12  In  an  arcjet  reactor, 
a  dc  arc  is  struck  in  the  feedstock  gas,  this  effluent  is  ex¬ 
panded  through  a  converging/diverging  nozzle  into  the  reac¬ 
tor,  and  diamond  grows  from  the  stagnation  point  flow  of 
reactive  gas  onto  a  substrate.  These  reactors  are  attractive 
candidates13-16  for  testing  the  models  of  diamond  CVD  be¬ 
cause  of  the  well-defined  gas  residence  time  in  the  reactive 
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plume  and  rapid  growth  rate.  The  gas  phase  chemical  com¬ 
position  of  a  dc-arcjet  plume  is  strongly  dependent  on  the 
atomic  hydrogen  concentration.  Coltrin  and  Dandy15  param¬ 
eterize  their  model  results  as  a  function  of  the  hydrogen  dis¬ 
sociation  for  the  very  wide  range  between  0.05  and  0.95.  To 
further  understand  the  CVD  process  and  to  test  the  existing 
models,  it  is  therefore  necessary  to  determine  the  absolute 
concentration  of  atomic  hydrogen  in  dc-arcjet  reactors. 

Atomic  hydrogen  has  been  measured  in  hot  filament  re¬ 
actors  using  nonintrusive,  in  situ  optical  detection  techniques 
including  laser-induced  fluorescence,17  resonance  enhanced 
multiphoton  ionization,18  and  third-harmonic  generation.19  It 
has  also  been  measured  by  mass  spectroscopy  of  the  molecu¬ 
lar  beam  sampled  through  a  hole  in  the  substrate.20,21  Indirect 
probe  measurements  have  determined  atomic  hydrogen  con¬ 
tent  in  hot  filament  reactors  by  measuring  the  heat  input  to  a 
probe  from  the  heterogeneous  recombination  of  atomic  hy¬ 
drogen  to  form  H2. 22,23  However,  measurements24-26  of 
atomic  hydrogen  in  arcjet  plumes  have  been  limited  to  very 
low  pressure  reactors  (few  mTorr)  to  study  arcjet  thrusters 
and  none  of  these  measurements  was  performed  during  dia¬ 
mond  deposition. 

We  report  the  determination  of  the  atomic  hydrogen 
mole  fraction  in  the  reactive  plume  of  a  dc-arcjet  reactor 
using  a  detailed  power  balance.  The  electrical  input  power  is 
balanced  with  the  power  losses  due  to  arcjet-nozzle  cooling, 
radiation,  ionization,  enthalpy  change,  and  directed  velocity 
of  the  gas  phase  species.  From  this  balance  the  hydrogen 
dissociation  rate  can  be  determined. 

II.  EXPERIMENT 

In  our  CVD  reactor,  schematically  shown  in  the  inset 
in  Fig.  1,  a  0.7  MPa  (90  psi)  mixture  of  H2/Ar  (0.8:1, 
7  slm  total  flow)  is  activated  by  a  dc  discharge  and  expands 
through  a  converging/diverging  water-cooled  nozzle  into  a 
reactor  maintained  at  constant  pressure  (nominally  3.3  kPa  or 
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FIG.  1.  Experimental  setup  for  gas  temperature  measurements  (case  a)  and  Doppler  velocimetry  (case  b)  in  a  dc  arcjet  using  laser  induced  fluorescence. 


25  Torr).  Methane  (—0.5%  of  the  H2  flow)  is  injected  into 
the  flow  in  the  diverging  section  of  the  nozzle.  The  lumi¬ 
nous,  chemically  reactive  gas  mixture  forms  a  supersonic 
stagnation  point  flow  onto  a  water-cooled  molybdenum  sub¬ 
strate  (38  mm  from  the  nozzle  exit),  where  diamond  grows  at 
an  average  growth  rate  of  approximately  50  /xm/h.  The  ex¬ 
tremely  high  quality  diamond  deposited  using  our  standard 
conditions  is  demonstrated  in  Figs.  2  and  3.  A  scanning  elec¬ 
tron  micrograph  (Fig.  2)  of  the  surface  morphology  shows 
the  (111)  texture  of  a  free-standing  500-/xm-thick  diamond 
plate,  and  Fig.  3  shows  a  Raman  spectrum  of  the  same  dia¬ 
mond  sample.  This  spectrum  was  taken  using  a  cw  argon  ion 


60  pm 


FIG.  2.  Scanning  electron  micrograph  of  the  CVD  diamond  deposited  in  the 
arcjet. 


laser  with  an  unfocused  300  pm  spot  size.  It  shows  a  char¬ 
acteristic  diamond  peak  at  1332  cm-1,  which  is  more  than  20 
times  the  nonresonant  background  luminescence.  There  are 
no  signs  of  graphitic  impurities.  The  2.2  cm-1  full  width  at 
half-maximum  (FWHM)  of  the  diamond  Raman  line  is  in¬ 
distinguishable  from  the  linewidth  observed  in  this  Raman 
spectrometer  from  natural  type  Ila  diamond. 

A.  Calorimetry 

We  determined  the  amount  of  dissociation  of  feedstock 
hydrogen  using  a  detailed  power  balance: 


FIG.  3.  Raman  spectrum  of  the  CVD  diamond  in  Fig.  2. 
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^ heating^”  ^accelerating)* 


(1) 


The  voltage  drop  in  the  arc  and  the  arc  current  were  mea¬ 
sured  to  determine  the  input  power.  This  is  balanced  with  the 
power  losses  of  the  nozzle  water  cooling,  radiation  from  ac¬ 
tivated  gases,  gas  ionization,  the  temperature  rise  in  the  gas, 
and  the  power  to  accelerate  the  gases  and  produce  the  di¬ 
rected  velocity  of  the  plume. 

Values  for  all  these  terms  except  gas  temperature  and 
gas  velocity  can  be  monitored  routinely.  At  our  standard 
conditions  the  arc  current  is  1 1.4  A  and  the  voltage  drop  is 
140  V  for  a  power  input  of  1600  W.  The  converging/ 
diverging  nozzle  is  water  cooled,  and  the  flow  and  tempera¬ 
ture  rise  of  this  coolant  are  monitored;  at  standard  conditions 
485  W  are  lost  to  the  cooling  water.  The  plume  of  the  arcjet 
has  atomic  line  and  molecular  band  emission  leading  to  a 
loss  term  due  to  radiation.  We  integrate  this  emission  and 
estimate  a  blackbody  radiation  loss  with  a  generous  upper 
bound  of  50  W.  Langmuir  probe  measurements27  show  that 
the  ion  density  in  the  plume  is  less  than  1013  ions/cm3,  re¬ 
sulting  in  a  negligible  power  loss  (1  W)  to  gas  ionization.  In 
the  high  pressure  arc,  the  ion  density  is  much  larger  than  this 
value.  However,  most  of  the  ions  and  electrons  apparently 
recombine  on  the  tungsten  surface  of  the  converging/ 
diverging  nozzle,  producing  part  of  the  large  heat  load  to  the 
nozzle.  Very  similar  ion  densities  were  measured  in  a  radia¬ 
tion  cooled  arcjet  thruster.26 


B.  Gas  temperature  and  velocity  measurements 

We  use  laser-induced  fluorescence  (LIF)  of  NO  seeded 
into  the  gas  plume  (0.1 1%  of  total  flow)  to  determine  the  gas 
temperature  and  velocity.  LIF  is  a  well  established  probe  of 
gas  temperature  and  velocity  in  combustion  systems.28,29  NO 
is  relatively  unreactive  in  this  environment  and  therefore  is 
not  expected  to  change  the  gas  temperature  or  the  flow  dy¬ 
namics.  LIF  measurements  of  CH  in  the  arcjet  plume  during 
diamond  deposition  and  NO  seeded  into  the  flow  give  iden¬ 
tical  gas  temperatures.30  Even  though  the  Langmuir  probe 
measurements27  show  an  increase  in  the  ion  density  when 
NO  is  added  to  the  flow,  this  is  a  negligible  perturbation  to 
the  power  balance.  Having  7  slm  of  gas  flowing  through  a  10 
mm-diameter  into  a  25  Torr  reactor,  we  expect  the  NO  mol¬ 
ecules  to  have  reached  the  mean  velocity  of  the  other  species 
by  the  time  they  pass  the  nozzle  exit,  where  the  velocity  is 
measured. 

A  tunable,  pulsed  excimer-pumped  dye  laser  (10  ns 
pulse  width,  —0.3  cm-1  bandwidth,  frequency  doubled  with 
a  BBO  crystal)  produces  the  light  to  excite  the  NO  A22  + 
-X2n  transition  near  226  nm.  The  laser  beam,  with  ener¬ 
gies  around  1  /xJ  to  avoid  saturation  of  the  spectroscopic 
transitions,  is  constricted  with  a  1-mm-diam  aperture  and 
aligned  through  the  arcjet  plume  perpendicular  to  the  di¬ 
rected  velocity.  The  fluorescence  is  collected  perpendicular 
to  both  plume  and  laser  beam.  The  collected  light  is  spatially 
filtered  with  an  aperture  to  1  mm  resolution  and  is  detected 
broadband  using  a  solar  blind  photomultiplier  (Hamamatsu 
R166H)  with  an  UG-5  colored  glass  filter  (Schott).  The  spec¬ 


FIG.  4.  Excitation  LIF  of  seeded  NO  A-X  (0,0)  for  gas  temperature  mea¬ 
surements  in  a  dc  arcjet. 


tral  bias  introduced  by  the  filter  was  less  than  3%.  A  boxcar 
is  used  to  time  integrate  the  fluorescence.  A  section  of  the 
excitation  spectrum  of  the  (0,0)  band  of  the  NO  A-X  tran¬ 
sition  is  shown  in  Fig.  4.  The  gas  temperature  is  determined 
from  the  rotational  temperature  derived  by  Boltzmann  analy¬ 
sis  of  such  spectra.  We  find  the  gas  temperature  to  be  2200  K 
at  the  exit  of  the  nozzle  in  the  middle  of  the  arcjet  plume 
with  a  statistical  uncertainty  of  100  K.  These  results  are  in 
good  agreement  with  rotational  temperature  measurements 
on  CH  B-X  previously  reported.31  Storm  and  Cappeili  found 
similar  gas  temperatures  in  1  kW  radiatively  cooled  arcjets.32 

The  directed  gas  velocity  of  the  arcjet  plume  is  deter¬ 
mined  from  Doppler  shift  measurements  of  LIF  of  NO 
seeded  into  the  gas  flow.  The  laser  beam  is  aligned  through  a 
3-mm-diam  hole  in  the  substrate  into  the  counter- 
propagating  gas  plume  as  shown  in  Fig.  1.  Before  the  laser 
beam  enters  the  CVD  reactor,  it  passes  through  a  reference 
cell  filled  with  NO  (—1  Torr).  A  second  detection  system 
collects  the  NO  LIF  signal  from  the  stationary  reference.  The 
LIF  signal  from  both  the  reference  and  the  arcjet  are  shown 
in  Fig.  5  for  excitation  of  the  Rx(  17.5)  of  the  NO  A-X  (0,0) 
transition.  Both  spectra  are  fit  to  a  Gaussian  shape  using  a 
Lervenberg-Marquadt  least-squares  algorithm.  The  Doppler 
shift  of  the  two  line  centers  gives  the  directed  velocity  v 
=  cAX/\,  where  c  is  the  speed  of  light,  AX  represents  the 


FIG.  5.  LIF  of  seeded  NO  in  the  reference  cell  (open  circles)  and  in  the 
arcjet  (closed  circles)  fitted  to  a  Gaussian  (line)  for  Doppler  velocimetry. 
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TABLE  I.  Power  balance  to  determine  the  fraction  of  hydrogen  dissocia¬ 
tion. 


Power  into  the  arcjet 

Power  supply: 

1600±32  W 

Power  losses 

Cooling  water 

485+55  W 

Accelerating  the  gases 

410±82  W 

Heating  the  gases 

330+50  W 

Radiative  loss  (estimated) 

<50  W 

Ionization 

<1  w 

Power  to  dissociate 

325+118  W 

-s-  Energy  to  dissociate  H2 

4.35  X105  J/moi 

Rate  of  atomic  hydrogen 

7.3X10-4  mol/s 

=>  Fraction  of  dissociation 

0.33±0.12 

Doppler  shift,  and  X  is  the  excitation  wavelength.  We  find 
the  directed  velocity  of  the  plume  to  be  v=  2600 ±  100  m/s  in 
the  center  of  the  plume  at  the  exit  plane  of  the  nozzle.  The 
bandwidth  of  the  laser  dominates  the  linewidth  in  Fig.  5  for 
the  reference  cell  data  and  is  a  significant  contribution  to  the 
linewidth  observed  in  the  arcjet  plume.  The  single  scan  mea¬ 
surement  of  both  the  Doppler  shifted  signal  and  the  reference 
signal  avoids  any  calibration  uncertainty  in  wavelength. 
Great  care  was  taken  to  avoid  problems  in  the  linearity  of  the 
laser  step  size  due  to  backlash  in  the  steppermotor  when 
scanning  back  and  forth.  Since  the  directed  velocity  is 
greater  than  the  speed  of  sound  in  this  system 
(^sound=1*2  km/s),  the  free  stream  velocity  is  not  perturbed 
by  the  hole  in  the  substrate.  Of  course,  any  measurements 
near  the  substrate  would  be  influenced. 

C.  Determination  of  atomic  hydrogen 

We  can  now  use  Eq.  (1)  to  determine  the  fraction  of  the 
feedstock  hydrogen  that  is  dissociated  in  the  arcjet  plume. 
We  define  the  dissociation  fraction,  X=0.5(nH//iH2),  where 
is  the  number  of  moles  for  atomic  hydrogen  and  is 
the  initial  number  of  moles  of  molecular  hydrogen  in  the 
feedstock  gas.  If  all  hydrogen  molecules  are  dissociated,  the 
fraction  of  dissociation  is  X=  1.  The  hydrogen  heat  of  for¬ 
mation  is  43.5  kJ/mole;  we  determine  X  from 

f*in~  ^cooling-  ^radiation ~~  ^ionization 

-(C^-Ar+^2-39.95X10-3)nAr, 

-  (c“2-  A7-=  \v2  •  2.016X  10_3)(  1  -  X>h2  , 

-(C”-A7’+|v2-1.008X10_3)2X-nH,  (2) 

where  Cp  is  the  specific  heat  capacity,  v  the  directed  gas 
velocity,  A  T  the  difference  between  room  temperature  and 
gas  temperature,  and  n  the  gas  flow  is  in  mol/s. 

We  find  the  fraction  of  hydrogen  dissociated  is  0.33 
±0.12  at  our  standard  diamond  deposition  conditions;  the 
power  balance  data  are  summarized  in  Table  I.  The  value  of 
0.33  is  a  lower  limit  because  we  assume  the  maximum  tem¬ 
perature  and  velocity  to  be  constant  over  the  entire  exit  plane 
of  the  nozzle.  If  we  calculate  the  fraction  of  dissociation  for 
the  average  temperatures  and  velocity  at  the  exit  plane  of  the 
nozzle,  we  find  a  dissociation  fraction  of  0.51,  in  good  agree- 


FIG.  6.  Power,  velocity,  and  fraction  of  atomic  hydrogen  vs  hydrogen  feed¬ 
stock  gas. 


ment  with  the  results  of  Probst  et  al2A  who  find  a  dissocia¬ 
tion  fraction  of  —0.50  for  a  1  kW  arcjet  thruster  operating  on 
pure  hydrogen  feedstock  in  a  reactor  at  pressures  of  a  few 
mTorr. 

As  the  ratio  of  hydrogen  to  argon  in  the  feedstock  is 
varied,  we  find  a  maximum  in  the  hydrogen  dissociation.  The 
variation  of  arcjet  power,  gas  velocity,  and  hydrogen  disso¬ 
ciation  is  shown  in  Fig.  6  as  a  function  of  the  fraction  of 
hydrogen  in  the  feedstock.  The  gas  temperature  remains  con¬ 
stant  within  5%  for  this  range  of  reactor  conditions.  The 
velocity  increases  steadily  with  increasing  molecular  hydro¬ 
gen  feedstock.  As  the  fraction  of  hydrogen  in  the  feedstock 
increases,  the  arc  voltage  increases  for  a  given  current, 
thereby  increasing  the  operating  power.  This  is  a  typical 
current-voltage  characteristic  of  a  supersonic  dc  arcjet  due 
to  the  negative  impedance  of  the  arc  column.33  At  50%  hy¬ 
drogen  feedstock,  the  dissociation  fraction  in  Fig.  6  shows  a 
maximum.  A  similar  relationship  was  measured30  for  the  CH 
concentration  as  a  function  of  the  H2/Ar  ratio;  the  CH  con¬ 
centration  increased  with  increased  hydrogen.  This  correla¬ 
tion  between  H  atoms  and  CH  molecules  agrees  with  the 
model  predictions13'15  that  H  and  CH  are  in  partial  equilib¬ 
rium  in  the  free  stream  of  the  arcjet  plume. 

The  method  presented  is  adequate  for  measuring  large 
dissociation  fractions  of  hydrogen  in  plasma  reactors  and  has 
some  very  important  advantages  over  quantitative  spectro¬ 
scopic  techniques  in  the  case  of  supersonic  gas  jets.  A  prob¬ 
lem  related  with  compressible  flows  is  the  uncertainty  of  the 
local  pressure  at  the  nozzle  exit,  which  usually  does  not  cor¬ 
respond  to  the  background  pressure.  As  a  consequence, 
knowing  the  degree  of  dissociation  requires  the  number  den¬ 
sity  determination  of  both  H  and  H2,  and  the  error  associated 
with  the  two  independent  concentration  measurements  is  at 
least  of  the  same  order  as  in  the  power  balance  method.24 

III.  SUMMARY 

A  detailed  power  balance  determines  the  dissociation 
fraction  of  hydrogen  in  the  plume  of  a  diamond  depositing 
dc  arcjet.  The  gas  temperature  at  the  nozzle  exit  is  measured 
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by  LIE  on  seeded  NO  to  be  2200  K,  and  the  directed  gas 
velocity  is  determined  to  be  2600  m/s  by  Doppler  velocim- 
etry.  We  find  hydrogen  dissociation  fraction  to  be  0.33  and 
0.5  if  peak  or  average  velocities  and  temperatures  are  cho¬ 
sen. 
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